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The Neutron Absorbing Isotopes in Gadolinium and Samarium 


R. E. Lapp,! J. R. VANHorn,? anp A. J. DempsTER 
Argonne National Laboratory, Chicago, Illinois 
(Received February 20, 1947) 


The large absorption for slow neutrons in samarium was found to be caused by the isotope at 
mass 149, and in gadolinium to the isotopes at masses 155 and 157, the cross section for 157 
being about three times that for 155. Mass spectra were made with the double focusing mass 
spectrograph, using the spark with nickel tubes packed with oxide as electrodes. Mass spectra of 
normal samples were compared photometrically with those of samples in which this isotope 
abundance had been altered by long exposure to slow neutrons. In the altered gadolinium the 
abundance of mass 157 was less than 155 instead of greater, and 156 became greater than 160 
rather than less. In samarium the abundance of mass 149 was reduced to be nearly as weak as 
148, and mass 150 was correspondingly increased. 





T was shown in 1935 by Dunning, Pegram, 
Fink, and Mitchell* that neutron absorption 
in samarium and gadolinium is extremely large. 
It is of interest to find which isotopes in these 
elements are responsible for the large absorption. 
The isotopes have been identified by measuring 
the.reduction in their intensity, produced by 
long exposure to neutrons in the Clinton Pile. 
The samples were prepared and exposed in the 
Clinton Pile by Dr. L. Borst in the form of very 
thin layers of oxide, 1.6 milligrams per sq. cm 
for the gadolinium and 1.9 milligrams per sq. 
cm for the samarium. With this thickness the 
inside of a sample receives almost the same 
bombardment as the outside layers. 
The isotopes were photographed with a double 
focusing _ mass spectrograph,‘ using a vacuum 
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spark and an electrode containing about 0.5 
milligram of the sample. The abundance of the 
isotopes was deduced from photometer traces 
made with a Moll recording microphotometer. 
The intensity standards were the normal isotope 
abundances as discussed below. 

For normal gadolinium the abundances meas- 
ured by Wahl® give a normal photographic 
density curve in terms of isotope abundance as 
shown in Fig. 1. The curve is made from five 
traces of different exposures on the same plate. 
The points shown on the curve are plotted with 
the abscissae marked on the five scales below. 
In the normal material the mass at 157 is more 
abundant than the one at 155, and the mass at 
160 is more abundant than the one at 156. 

With the exposed gadolinium, the abundances 
were radically altered, and it was necessary to 
take the new abundance ratios shown in Fig. 2 


* Wahl, Soc. Sci. Fenn. 11, 1 (1941), According to a 

paper by J. Mattauch and H. Ewald, Zeits. f. Physik 122, 
eI (1944), the standard intensity marks were made by 
x-rays and a rotating sector. 
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A sample of samarium oxide was exposed in 
the same way. Here there is the difficulty that 
the normal abundances are not accurately known 


-°* Normal Gadolinium 





Normal Samarium 











Log Abundence 


Fic. 1. Photographic density curve made from micro- 
‘mvonne = mae traces of five mass spectra on one plate. The 
ogarithmic isotope abundance ratios are given by the 
five scales. 














Exposed Gadolinium 

















Log Abundence 


Fic. 3. Density-abundance curve drawn from micro- 
photometer traces of five mass spectra of normal samarium. 
The isotope abundance scales are shown on the five 
lines below. 








Exposed Samarium 























Log Abundance 
Fic. 2. Density curve plotted from microphotometer 
traces of five mass spectra of exposed gadolinium. The 


isotope abundance scales are indicated on the five lines 
below. One normal scale is shown for comparison. 


to make the observed densities fit on the pho- 
tographic density curve. The points shown on - aed 
the curve are plotted with the abscissae marked 
on the five altered abundance scales. The normal Saanyoeny eee 
scale is shown for comparison. The isotopes at wou? 
masses 157 and 155 are reduced, and those at 
158 and 156 are increased. The normal and altered 
abundances in percents are given in Table I. 

We may conclude that the absorbing isotopes 
are at masses 155 and 157 with the absorbing 

Fic. 4. Density curve of five mass spectra of samarium 


. section for 157 approximately 3.5 times exposed to neutrons. The o— isotope abundance scales 
as large as that for 155. are given below with one f®rmal scale for comparison. 
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as Dr. Aston’s estimates made from incompletely 
resolved spectra are the only ones now available.*® 
New abundances were deduced from an indirect 
comparison with the gadolinium, in which the 
densities from standard light intensities served 
as a reference in both cases. These new normal 
abundances are used in plotting the photographic 
density curve in Fig. 3, made from five spectra 
with different exposure times on one plate. 

The abundances of the isotopes in the exposed 
sample were altered by a reduction in the isotope 
at mass 149 and an increase in the mass at 150. 
This is shown in Fig. 4 where the new abundances 
shown on the lines at the bottom for the five 
spectra are those required to make the observed 
densities fall on one photographic density curve. 
The normal abundances of Fig. 3 are shown for 
comparison in the top scale. Although the 





‘F. W. Aston, Proc. Roy. Soc. A146, 46 (1934). 
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TABLE I. 








Isotope 15S 156 157 158 160 


Gd (normal) 15.6 20.6 16.42 23.45 20.87 
Gd (altered) 13.6 22.3 9.86 30.60 20.62 


Change —-2.0 +17 —-6.6 +7.1 











changes are less than in the case of gadolinium, 
the samarium isotope at mass 149 has become 
less abundant, from 12.8 to 10.0 percent, and 
the isotope at mass 150 has increased from 5.0 to 
7.1 percent. We may conclude that the large 
absorption in samarium is due primarily to the 
isotope at mass 149. 

The authors are indebted to Dr. L. Borst for 
preparing the thin deposits and exposing them 
in the Clinton Pile. 

These observations were carried out while the 
authors were associated with the Metallurgical 
Laboratory at the University of Chicago. 
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Note on the Nuclear Electric Quadrupole Spectrum of a Homonuclear Diatomic 
Molecule in a Magnetic Field* 


H. M. Foiey 


Columbia University, New York, New York 
(Received February 20, 1947) 


It is shown that for very large values of the rotational quantum number J the energy of 
interaction of the electric quadrupole moments of a homonuclear diatomic molecule may be 
expressed as the sum of the energies of states of the single nucleus problem. This approxima- 


tion becomes less accurate for any J value as the applied magnetic field is decreased. 


HE effect of the nuclear electric quadrupole moment in the magnetic spectrum of a hetero- 
nuclear molecule has been given by Feld and Lamb.' In this work the molecule was treated 


as if the two nuclei were entirely independent. The calculations were actually carried out for the 
case of a single nucleus interacting with the molecular field. 

In the case of a homonuclear molecule certain apparent complications arise due to the exchange 
degeneracy and to the fact that both nuclei are coupled to the molecular rotational angular mo- 
mentum through the quadrupole interaction. The Hamiltonian expression for two nuclei with 


quadrupole moments in a magnetic field is 


e*qiQi 





~ 2I(2I —1)1,(2hh—1) 


{3(ieJ)?+3(hieJ) -1’J?} + 
* & [3(leeJ)?+$ (hee J) — he? J*} + wogiHels + wogeH ele + uogs HJ, 


* Publication assisted by the Ernest Kempton Adams Fund for Physical Research of Columbia University, New York. 
'B. T. Feld and W. E. Lamb, Phys. Rev. 67, 15 (1945), hereinafter referred to as FL, 


e*g2Qe 
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in which Q; and Q2 are the quadrupole moments, I, and I; the spin vectors of the nuclei, and a 
and gq: are directly related to the electric field gradients at the positions of the nuclei.! It is assumed 
that the energy of interaction between the magnetic moments of the nuclei, and of the molecular 
rotational magnetic field with the nuclear magnetic moments, is small compared to the quadrupole 
energy. 

For the discussion of the energy levels of a homonuclear molecule it is convenient to employ a 
“strong-field” representation, in which the rotational angular momentum, the spins of the individya| 
nuclei, and the components of each of these quantities in the direction of the magnetic field are 
diagonal matrices. In the absence of the quadrupole interaction all states with the same value of m, 
and of Msg, the component of total nuclear spin in the magnetic field direction, are degenerate. 

If the wave function representing a nucleus with z component of spin given by 7 is a;, then the 
part of the total nuclear wave function which depends on the nuclear spin coordinates may be 
written a;(1)-ox(2) where 1+k= Ms. It is useful to alter the representation by forming the sym. 
metric and antisymmetric combinations of these functions 


1 
lT'a* =! (1) +a, (2) +a (1)a;(2) } » Te=ai(1)-a,(2). 


The magnetic terms in the Hamiltonian are diagonal in this representation with energy values 
g:M sH,+ gsm sH,. 


The matrix elements of the quadrupole interaction for an individual nucleus have been given by 
Kellogg, Rabi, Ramsey, and Zacharias.? Non-vanishing matrix elements exist between states with 
the following relationships of quantum numbers 























Ams= 0 Am;= 0 
Ams=+1 Am;=-—1 
Ams=-—1 Am;=-+1 
Ams=+2 Am,;=—2 
Ams=—2 Am ;=+2 








In the case of a homonuclear diatomic molecule there are no matrix elements of the two quadrupole 
terms between the states defined in the I representation having a common value of Ms. Thus if 
the quadrupole interaction is considered to be a small perturbation (strong-field condition), the I 
representation is the correct set of wave functions in first-order approximation. The diagonal matrix 
element of the quadrupole operator for a single nucleus is given by 


Fim = 4[3#?— h(i +1) 13m?—J(J+1)]. (2) 
For two identical nuclei in the T representation the diagonal element is 
tm km 
Fimg+ Femy =[3i2+38?— I, (; +1) 3m s?—J(J+1)]. (3) 


The total nuclear spin is not diagonal in the ['y functions, that is, the wave functions in terms of 
which the quadrupole energy is a diagonal matrix consist of combinations of functions corresponding 
to different values of total nuclear spin. The first-order quadrupole energy consists simply of the 
sum of the energies of two individual nuclei. It will be shown below that if the rotational angular 
momentum is very large the eigenvalues of the magnetic and quadrupole energies may always be 
expressed as the sum of two eigenvalues ‘of the individual nucleus problem.’ 

The Hamiltonian of Eq. (1) may be written as 


i = Hy (I,JH)+H2(I2JH). (4) 
2J. M. B. Kellogg, I. I. Rabi, N. F. Ramsey, and J. R. Zacharias, Phys. Rev. 57, 677 (1940). 
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(The rotational moment term is not important for the following discussion.) The matrix elements 
for an individual nucleus in the strong-field representation are 





i my ‘ 
VajRmyrimy:, my +i-i+ VM joi, ; (5) 





2 — . : i_my 
in which Vm; is the magnetic interaction energy (diagonal), and Vo;Rmy-dmy-, my +i-tis the quadrupole 
matrix element between nuclear spin states 7 and j and rotational states my and my. The energy 
matrix of Eq. (4) in the I’ representation is then 











tkm J i_my i k my k 
Halmy: = VonRm sdk, lémy, my+i-n+ Vm,6; nok, Vaud, nRmybmy:, my+k—-i+ Vmidx, 161, ne (6) 






In this matrix non-zero elements appear only between states of equal values of M, the zs component 
of total angular momentum. Within each sub-matrix for any value of M there are no matrix elements 
between states of the same Ms. In such a sub-matrix the value of the rotational magnetic quantum 
number has at most 4J;+1 values distributed in the interval M—2I,<m;<M+2],. If J is very 
large compared to J;, the variation of the matrix elements Rnj, with m,, with m;—m~y fixed, 
will be small in this interval. (This will not_be the case for the small fraction of states in which my 


or J—m, is comparable with J;.) If this variation be neglected Rx4, depends only on 
















My—My=n—-t Or mMy—my=l—k 








and the matrix element Eq. (5) may be written V;*. The homonuclear matrix element of Eq. (6) 
becomes 






Hi; = Vn "bk, 1+ Vi" dy, n. 7 (7) 






In the matrix element V,‘ the rotational matrix element R%4, may be set equal to its average value 
for which m, is equal to M. 

We assume that the energy matrix of the individual nucleus problem is brought into diagonal 
form by a transformation matrix A. Thus W=AVA~ in which W is a diagonal matrix. The direct 
product matrix 







Q@=AXA or ay™=a;/"a,' (8) 







is to be applied to the complete energy matrix of Eq. (7) 





He =E DL anHaa-) py 


nl ik 





=D X asmay'{ Vn'd:, + Vi*din} (a-") p*(a") 


nl ik 


= W514, 5m, pt Wo'dm, 951, 0° 















Thus the transformation matrix @ diagonalizes the energy matrix, and the diagonal values are 
just the sums of eigenvalues of the individual nucleus problem. 

It is clear that any system consisting of two (or more) parts, the energy matrix of which can be 
put into the form of Eq. (7), will show this same additivity of the energies of the individual parts. 
Thus a heteronuclear molecule of large rotational angular momentum, and with quadrupole inter- 
actions, will yield eigenvalues expressible as the sums of the energies of the individual nuclei. The 
inclusion of an interaction between the nuclei will of course destroy the form of the energy matrix 


of Eq. (7). 
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A weak oscillating magnetic field H’ introduces a perturbation of the form 
KH’ = wogrH’ 1, + wogrH’ el, 


leading to transitions between the levels given by Eq. (9). From the properties of the matrix @ it 
may be shown that the possible transitions are restricted to changes in the value of only one of the 
terms of Eq. (9), and the selection rules for these transitions are identical with those in the spectrum 
of a single nucleus. Subject to the approximation that the rotational angular momentum is large, 
we have the general rule: the magnetic nuclear spectrum of a homonuclear diatomic molecule, 
including electric quadrupole interaction, is identical with that of a single nucleus with the same 
spin. This rule is not in agreement with the recent results of Feld.* 

The physical meaning of Eqs. (3) and (6) may be expressed quite simply. In a strong magnetic 
field the two spin vectors and the rotational angular momentum are decoupled and precess inde- 
pendently in the field. Equation (3) expresses the fact that the energies of the precessing spin vectors 
are additive. As the field is reduced to zero, the spin vectors are each coupled to the rotational 
angular momentum through the quadrupole interaction, and the three vectors form a resultant 
total angular momentum. If the rotational vector J is very large the spin vectors may be considered 
to be precessing about J. In this case the rotational vector has the role of a strong applied field 
and the energies are additive (Eq. (6)). 

The effect of a finite value of rotational angular momentum is to make important the variation 
of the rotational matrix elements with m, within a group of states of fixed M. In this case the com- 
plete matrix of Eq. (6) may not be diagonalized by a transformation of the form @=A XA with 
eigenvalues expressed simply as sums of eigenvalues of the individual nucleus problem. The strong 
field first-order solution of Eq. (3), however, shows the additive property even for small values of J, 
and it is clear from the foregoing discussion that the deviations from additivity become more im- 
portant in higher approximations. To investigate the maximum effect we proceed directly to the 
zero-field case. The following discussion will deal only with the quadrupole interaction in zero field. 
In this case we may use a representation in which F, mr, I, J, I;, Iz are diagonal. F is the quantum 
number of total angular momentum, and m, is its component on a fixed space axis. J is the total 
nuclear spin. The matrix elements of the quadrupole interaction are given for reference in the 
Appendix. Non-zero matrix elements exist between states differing in total spin by AJ=0+2. 
For any value of ip = F—J the states of total spin values J>i%) must be diagonalized. The combina- 
tions of total spin wave functions in which the quadrupole perturbation is diagonal for any value 
of io, are the same as the total spin combinations in the strong field case for an identical value of Ms. 

For a single nucleus, in the zero-field case, Feld and Lamb gave an expression (FL III 7) for the 


energy which is accurate for large J values 
Eq=eq0[3i?—L(hi+1))/4h(2h-1), in=F-J. (10) 


These authors also gave the exact solution of the single nucleus problem (FL III 5). The eigenvalue 
solution in the homonuclear problem described above may be expressed as the sum of two expressions 
of the type of Eq. (10) for large J values, in accordance with the general rule (Eq. (9)). For smaller 
J values, ahd particularly for large values of nuclear spin, deviations from simple additivity will 
occur. Investigation of several special cases indicates that these deviations are of the same magnitude 
as, though not identical with, the differences between the exact single nucleus solution (FL III 5) and 
the values given by Eq. (10). For nuclear spins as large as J=5 and with J=20 these deviations 
are important, and the simple formula of Eq. (10) will not apply accurately in either the homonuclear 


or heteronuclear case. 
The author wishes to thank Professor I. I. Rabi for several interesting discussions of this problem. 


3B. T. Feld, Phys. Rev. 70, 112A (1946). 
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APPENDIX 


The matrix elements of the quadrupole interaction in Eq. (1) in the weak-field representation in 
which F, mr, I, I:, Iz, J are diagonal are given here. The matrix elements of (I,*J) are* 


1 1,(1+1) — L222 +1) +1(I+1) 
(JIF |hieJ | JTF) => 2II+1) { F(F+1)—I([+1) -J(J+)}, 
I tF)=- COE eee 
(JIF|LeJ|JI-1) => 41?(2I —1)(2I+1) 


x {(F+I—J)(F+J—I+1)(F+I+J+1)U+J—F)}! 








=(JI-—1F|lyeJ|JIF), 
1((0+1-ht+h)(0+14+h-h)(htht+l+2)(h+h-D}\! 
ererere ns 4(I-+1)*2T+1)(2I+3) | 
x {(FEI+1—J)(F+J—D(F+I+I+2)I+14J-F)}} 
=(JI+1F|IyeJ|JIF). 





The off diagonal elements for (I:*J) are minus the corresponding elements in (I,J). 
In the homonuclear case J;=J2, and the matrix elements of the quadrupole terms of Eq. (1) are 
obtained by matrix multiplication : 


e*gQ 

= F(F+1) —I1([+1) —J(J +1)? 
(JIF|Hg|JIF) 2707 —) hh —1) a(F(F+1) —1(1+1) —J(J+1)] 
3 (2h+1)?-F 


8 4/°-1 








(F+I—J)(F+J—I+1)(F+I+J +1) U+J—F) 





3 (24+1)?—U +1)? 
8 (22+1)(2I+3) 





(F+I+1—J)(F+J—D)(F+I+J+2)I+1+J—F) 


+30 F(F+1) —1(+1) - I+) 1-2) IT+}, 





(JIF|Hg|JI+2F)= 


3e%Q | tt ~J)(F+J-1) 


2I(2J—1)I(2,—1) | 4(27+1)(27+3) 
(2h,+1)?— (+2)? 


‘ 
I+J+2)I+J+1—-F F+I+2-J 
Ee | 4(aI+3)Qi+s) 





x (F+J—I—1)(F+I+J+3)I+2+J—F) 
=(JI+2F|He|JIF), 
3e2gQ (24,+1)?—I2 
Te] 4(4I?—1) 
4/ (21,+1)?—(—-1)? 
4(2I —3)(2I—1) 





(JIF|Hg|JI—2F) = (F+1—J)(F+J—I+1)(F+1+J+1) 





4 
x(I+J-P)| (FHI-J-1)(F+J-N(F+I+DU+I-F-1)| ; 





*E. U. Condon and G. H. Shortley, The Theory of Atomic Spectra (Cambridge University Press, 1935), p. 71. 
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Diffraction of Neutrons by a Single Crystal 


W. H. Zinn 
Argonne Laboratory, University of Chicago, Chicago, Illinois 
(Received February 26, 1947) 


The intensity of neutrons emitted by the nuclear chain reactor has been shown to be suffi- 
ciently great to permit the observation of the Bragg reflection of neutrons by a single crystal. 
A single crystal monochromator has been constructed and used for the investigation of reso- 
nances. The well-known resonance in Cd has been measured and the Breit-Wigner constants 


determined. 





INTRODUCTION 


NUMBER of experiments designed to 

demonstrate the diffraction of slow neu- 
trons by crystals have been made in the years 
following the discovery of the neutron. Early 
methods of treating the interaction between the 
neutron and the nucleus, developed by Fermi! 
and extended by Wick,’ anticipated difference 
in scattering in crystalline and disordered sub- 
stances. Experimental verification was attempted 
by Preiswerk and von Halban* who measured 
the temperature variation in the distribution 
of neutrons transmitted by a polycrystalline 
medium, and Mitchell and Powers,‘ as well as 
Preiswerk,' presented evidence for a coherent 
component of a neutron beam scattered from an 
array of crystals. 

More recently, further observations of neutron 
diffraction have been made in the measurements 
of slow neutron scattering by Whittaker and 
Beyer® and in the work of Anderson, Fermi, and 
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CALCITE CRYSTAL 





Fic. 1. First experimental arrangement. 


1E. Fermi, Ricerca Scient. 7, pt. 2, 13 (1936). 
2 Wick, Physik. Zeits. 38, 403 (1937). 
a 6 Halban and Preiswerk, Comptes rendus 203, 73 
‘ Mitchell and Powers, Phys. Rev. 50, 486 (1936). 
5 Preiswerk, Helv. Phys. Acta 10, 400 (1937). 
* Whittaker and “e Phys. Rev. 55, 1101 (1939); 
Beyer and Whittaker, Phys. Rev. 57, 976 (1937). 


L. Marshall’ and of Rainwater and Havens? who 
observed variations with the energy of the 
neutron beam in the scattering of crystalline 
substances. The high neutron flux generated by 
the chain reacting pile has permitted, as re. 
ported here, a series of experiments which show 
the existence of an intense diffracted beam of 
neutrons from a single crystal, and which demon- 
strate that this beam is of sufficient intensity for 
the study of the energy dependence of neutron 
reactions by methods similar to those of x-ray 
spectrometry. The work of this paper has been 
reported briefly® and similar results have been 
given by Borst e¢ al.'° 


DESCRIPTION OF APPARATUS 


Neutrons diffusing from the main part of the 


chain reacting pile are permitted to enter a 


large block of graphite which acts as a slowing 
down column, bringing the average neutron 
energy to the equilibrium value for graphite at 
room temperature. Since the de Broglie wave- 
length of these neutrons is of the same order of 
magnitude as that of x-rays, it was possible to 
attempt the diffraction experiments with the 
customary crystals, using a spectrometer similar 
to that used with x-rays. As is shown schemati- 
cally in Fig. 1, a collimated beam of neutrons 
was obtained by introducing two thick cadmium 
slits, 5 meters apart, in the neutron beam from 
the graphite thermal column of the heavy water 
pile. Cadmium slits sufficed here since the neu- 


7 Anderson, Fermi, and L. Marshall, Phys. Rev. 70, 
102 (1946). 

® Havens and Rainwater, Phys. Rev. 70, 154 (1946); 
Rainwater and Havens, Phys. Rev. 70, 136 (1946). 

9W. H. Zinn, Phys. Rev. 70, 102A (1946). 

1°L. B. Borst, A. J. Ulrich, C. L. Osborne, and B. 
Hasbrouck, Phys. Rev. 70, 108A (1946); Phys. Rev. 70, 
557 (1946). 
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SINGLE CRYSTAL DIFFRACTION OF NEUTRONS 








PILE SHIELD 











Fic. 2. The crystal spectrometer. 


tron radiation from a thermal column is singu- 
larly free from fast neutrons. The collimated 
beam fell upon the (1, 0,0) planes of a large 
single calcite crystal, and the diffracted beam 
was detected by means of a BF; proportional 
counter filled, for greater detection efficiency, 
with gas enriched in the B"® isotope. Only early 
experiments were performed at the thermal 
column; subsequently the apparatus was set up 
at a side hole of the reactor where the high 
energy neutron flux is much greater than at the 
thermal column. 

The final form of the spectrometer (Fig. 2) is 
identical in principle with the usual x-ray instru- 
ment; however, in order to support the heavier 
collimators and detectors which are necessary for 
neutron diffraction experiments, the unit is made 
very massive. To improve the resolution the arm 
is long, having a length of 80 inches and a width 
of 8 inches. A 6-inch diameter crystal table is 
provided for mounting the crystals. The two 
graduated circles indicating the position of the 
arm and crystal table are 12 inches in diameter 
and are calibrated in degrees, while greater 
precision is attainable from readings on a cylin- 
drical scale calibrated in quarter minutes and 
placed at the drive wheel. The table and arm 
are driven by means of large knurled wheels 





through a train of gears, and in use can be 
adjusted independently or coupled by a friction 
clutch. The table and arm are driven in the two 
to one angular ratio necessary for reflection 
experiments. 

In order to measure with greatest efficiency 
the flux of diffracted neutrons, the beam was 
passed axially through the proportional counter. 
In this manner it was possible to detect nearly 
100 percent of all the neutrons of the low 
energies. The counter was 5 cm in diameter, 
60 cm long, and filled to a pressure of 40 cm of 
Hg with boron trifluoride. This gas was enriched 
about five times in the B® isotope and thus 
transmitted but 5 percent of neutrons of kT 
energy. The required degree of collimation was 


.obtained by inserting into the shield of the 


reactor a steel block containing a long channel. 
This channel was approximately 4” wide, 1” 
tall, and 8’ long. The channel was directed at a 
part of the reactor containing heavy water and 
uranium. The beam collimated by this means 
contains neutrons rangingin energy from thermal 
to those of the fission spectrum itself. The in- 
tensity obtainable is a sharp function of the 
degree of collimation desired and the spectral 
region to be investigated. At the maximum of 
the thermal spectrum, with quite excellent reso- 
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Fic. 3. Rocking curve (1, 0, 0) LiF crystal. 


lution, counting rates of 5000 counts per minute 
can be obtained. At the suggestion of W. H. 
Zachariasen, roughening the crystal surfaces was 
tried and gives an improvement in intensity of 
about a factor of two. At energies in the so-called 
resonance region considerably smaller counting 
rates are obtained. A detector more sensitive to 
neutrons in this region is very much to be 
desired. 


RESULTS 


A typical rocking curve is shown in Fig. 3. 
No attempt has been made to correct these data 
for resolution or detector sensitivity. The de- 
tector, boron, follows a 1/v law for absorption, 
and as a result the sensitivity of the counter is 
not uniform over the entire range of velocities. 
Therefore, the rocking.curve is not exactly sym- 
metrical, showing relatively higher values on the 
low energy side of the diffraction peak. 

Sufficient intensity of diffracted neutrons ex- 
isted throughout the range of energies emitted 
from the thermal column to make measurement 
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of its neutron spectrum possible by this method 
Intensities at low energies, since they are cea 
taminated by higher order components and are 
more efficiently detected, show considerably 
higher values than the real distribution in this 
region. Figure 4 shows the distribution of ney. 
trons from the thermal column as obtained with 
a calcite crystal. Figure 5 shows the spectrum of 
neutrons directly from the reactor. These ney. 
trons have not diffused through any considerable 
column of moderator material. Again no cor. 
rection to the data has been made for detector 
sensitivity or for resolution. Neither of the cor. 
rections is sizable in the low energy region 
covered in these curves. Conversion from glane- 
ing angle to energy was made by the relation 


ev = (2.22 X 10-*) /sin?}@, 


where @ is the glancing angle. Measurements 
were made in the energy interval.between 0.004 
ev and 0.30 ev, which is in the range of glancing 
angles between 50° and 2°. Both spectra show 
a strong, approximately Maxwellian component, 
distorted by high order contributions on the 
long wave-length side of the maximum. The 
spectrum of neutrons from the wall of the 
reactor, having passed through a smaller amount 
of moderator, has an additional component at 
small glancing angles representing neutrons not 
yet in equilibrium with the moderator. Attempts 
to fit the experimental data theoretically have 
been made by Goldberger and Seitz" for various 
moderator temperatures and are included on the 
figures. 

» To demonstrate that the diffracted beam con- 
sists primarily of neutrons in a very small energy 
band, a measurement of the transmission of a 
Pyrex plate was made over the range 0.018 to 
0.5 ev. This plate, previously standardized by 
means of a mechanical velocity selector, was 
shown to have a total cross section, because of 
its boron content, very close to 1/v in the region 
investigated. With the calcite crystal, measure- 
ments were made of the transmission of the 
plate from which effective neutron velocities 
were obtained by comparison with the measure- 
ments made on the mechanical velocity selector. 


Goldberger and Seitz, Phys. Rev. 70, 116 (1946). . 
| —— discussion of their method will be published 
shortly. 
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4 , ; ; ; aoe Reciprocals of these velocities as a function of 
the sine of the glancing angle are plotted ag 
j , a ; c 
Py 4 points in Fig. 6. The dotted curve gives the 
| fis reciprocal velocity as calculated from the glance. 
ff ing angle and the crystal spacing. Similar meas. 
ab va . urements show that the straight line behavior - * 
Pal shown extends to 15 volts which was the upper 
| P ; limit of the instrument at the time these meas. 
+ Ps y urements were made. 
“ a“ The total cross section of cadmium has been 
af - ' measured by the method of the modulated 
FI F | cyclotron beam” and a single strong level ob. ‘ 
P served. The results obtained for cadmium with 
ash F 4 ; the crystal monochromator are shown in Fig, 7, 
at a : The constants of the resonance are similar to 
v those of Rainwater and Havens* but show a 
mS v4 7 greater maximum than those of Baker and 
Fs Bacher. The measurements of Fig. 7 extend to 
o4- PA 4 
iz Pe 15 electron volts, and no further resonance 
2. « ——+-_—+-—+ structure is found. A series of six foils of cadmium . 
. vn Nn 
4 40% om sec of thickness varying from 0.000142 x 10” to : 
0.11510" atoms per square centimeter were 0 
Fic. 6. Boron absorption as a function of neutron energy. ; possi 
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Fic. 7. Total cross section of cadmium. Solid curve is Breit-Wigner fit to experimental points. 
¢ in units of 10-* cm? per atom. 


2 C, P. Baker and R. F. Bacher, Phys. Rev. 59, 332 (1941). N 
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between 0.50 and 0.75 to reduce effects of high 
energy beam components. The experimental re- , 
sults are best fitted with a curve of the Breit- 
Wigner function having the constants: Ey) = 0.180 
ev, ao= 7800 X10 cm?/atom, !' =0.122 ev, and 
g,=6.0X 10-* cm?/atom. 

Acknowledgment is hereby made for the in- 
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valuable advice and skill in construction of the 
instrument contributed by T. J. O’Donnell, 
D. DiCostanzo, and J. Getzholtz of the labora- 
tory work shop. 

The initial results reported here were obtained 
in July-September, 1944. This work was carried 
out under the auspices of the Manhattan District. 
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Measurement of Neutron Cross Sections with a Crystal Spectrometer*t+ 
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By use of the monoenergetic beams of neutrons diffracted 
from the (100) planes of lithium fluoride, survey measure- 
ments of the total neutron cross section of a series of ele- 
ments have been made. The strong diffracted beams were 
obtained by placing a single crystal in the high flux of 
neutrons emitted by the heavy water moderated pile at the 
Argonne Laboratory. Observations on this group of ele- 
ments, which includes several strongly absorbing rare 
earths, and elements commonly used as neutron filters and 
detectors, have been made in the range of energy between 
0.04 and 65.0 electron volts. The method permits energy 
dependent cross-section measurement by activation or by 


I, INTRODUCTION 


ARLY experiments of Szilard,! of Amaldi 

and Fermi,? and of Tillman and Moon® on 
the energy variation of neutron absorption 
cross section, indicated that not all elements 
showed a simple dependence of cross section upon 
the energy of the reacting neutron. Several 
techniques have subsequently been applied in an 
effort to study the nature of this function, be- 
cause of its relation to the energy levels of the 
compound nucleus, in the limited region of low 


*Some of these results have been re previously : 
Sturm and Turkel, Phys. Rev. 70; 103 (A) (1946). 

This document is based on work ormed under 
Contract No. W-31-log-eng-38 for the Manhattan Project 
and the information covered therein will appear in Division 
IV of the Manhattan Project Technical Series as part of 
the contribution of the Argonne National Laboratory. 

** Now at the University of Wisconsin. 

1 Szilard, Nature 136, 150 (1935). 

* Amaldi and Fermi, Ricerca Scient. 1, 310 (1936). 

* Tillman and Moon, Proc. Roy. Soc. A153, 476 (1936) ; 
Nature 135, 904 (1935). 





transmission. However, all values reported were obtained 
by calculation from observed transmission data because the 
range of available intensities in the diffracted beam was 
best measured with a boron counter. Correction methods 
for resolution and order effects have been developed which 
are valid for low neutron energies. The following absorption 
levels have been observed: Rh, 1.28 ev; Au, 5.4 ev; Ir, 
0.635, 1.35, and 6.0 ev; Gd, 0.031 ev; Sm, 0.096, 10.0, and 
33.0 ev; Eu, 0.465, 3.3, 9.2, and 22.0 ev. Resonance levels at 
energies less than zero were found in Dy at —1.01 ev and 
in Eu at —0.011 ev. 


neutron energy. By employing filters and de- 
tectors of substances containing single strong 
resonance levels'~ it is possible to obtain an idea 
of the cross section of an element at the energy of 
the absorbing or detecting resonance. The obvi- 
ous limitations of this method include the 
uncertainty in the knowledge of the position of 
the absorbing resonance level and the random 
nature of the distribution of the resulting data on 
the energy scale. 

A method for determining the region of se- 
lective absorption for slow neutrons, based on the 
assumption that the capture cross section for 
boron and lithium is inversely proportional to the 
neutron velocity, has been suggested by Weeks, 
Livingston, and Bethe and by Frisch and 


‘ Bjerge and Wescott, Proc. Roy. Soc. A150, 709 (1935) ; 
Pontecorvo, Ricerca Scient. 6-7, 145 (1935); Bethe, Rev. 
Mod. Phys. 9, 113 et seq. (1937); Hanstein, Phys. Rev. 59, 
489 (1941); Horvath and Salant, _— Rev. 59, 154 
(1941); Feeny, Lapointe, and Rasetti, Phys. Rev. 61, 469 
(1942); Yalow and Goldhaber, Phys. Rev. 68, 99 (1945). 
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Placzek. According to this method, the quanti- 
tative measure of the absorption of neutrons in 
boron is taken as a measure of the energy of the 
captured neutron. Using this procedure, Gold- 
smith and Rasetti® have been able to locate with 
fair accuracy the position of the predominant 
levels in a series of strong absorbers. 

A third indirect technique for observing the 
energy dependence of capture cross section has 
been employed by Powers, Fink, and Pegram’ 
and others. By varying the temperature of their 
neutron moderator they were able to demonstrate 
an increase in cross section with decreasing 
neutron velocity for all elements tested. Since the 
method depends on evaluation of average cross 
section over a Maxwellian distribution of neutron 
energies having differing peak values, it cannot be 
expected that this type of measurement should 
do more than exhibit small variations in cross 
section due to the 1/v absorption in this region. 

The earliest direct and continuous measure- 
ments of cross section as a function of neutron 
energy over a small range of energies was made 
possible by the rotating shutter velocity selector 
of Dunning, Pegram, Fink, Mitchell, and Segré.® 
The range of this type of apparatus is limited 
to neutrons having energies less than 0.30 ev 
because of the nature of the cadmium cross 
section, yet with a sufficiently strong neutron 
source instruments of this type can be used for 
quite adequate measurements in this energy 
range. 

With the modulation of the cyclotron beam by 
Alvarez’ it was possible to observe the energy 
dependence of neutron cross section in still an- 
other manner. This technique, improved and ex- 
tended by many observers,’ has most recently 
been applied to the problem by Rainwater and 


5 Weeks, cavinapaa, and Bethe, Phys. Rev. 49, 471 


(1936); Frisch and Placzek, Nature 137, 357 (1936). 

* Goldsmith and Rasetti, Phys. Rev. 50, 328 (1936). 

7 Powers, Fink, and P m, Phys. Rev. 49, 650 (A) 
(1936); Libby and Long, Phys. Rev. 52, 592 (1937). 

Dunning, Pegram, Fink, Mitchell, and é, Phys. 
Rev. 48, 704 (1935). Fink, Dunning, Pegram, and Mitchell, 
Phys. Rev. 49, 103 (1936); Rasetti, é, Fink, Dunning, 
and Pegram, Phys. Rev. 49, 104 (1936); Fink, Phys. Rev. 
50, ‘738 (1936). 

® Alvarez, Phys. Rev. 54, 235 (A) (1938); Alvarez, Phys. 
Rev. 54, 609 (1938). 

10 Fertel, Gibbs, Moon, Thompson, and bf seeming — may 
Proc. Roy. Soc. 175, 316 (1940); Baker and Bacher, Phys. 
Rev. 59, 332 (1941); Manley, Haworth, and Luebke, 
Phys. Rev. 61, 152 (1942); Bacher, Baker, and McDaniel, 
Phys. Rev. 69, 443 (1946). 


Havens" who calculated the energy of the 
reacting neutron by measuring its time of flight 
over a fixed distance and evaluated many crogs. 
section functions in this manner for energies up 
to 1000 ev. ‘ 

The present method, utilizing the wave prop- 
erties of the neutron, brings techniques analogous 
to those of x-ray diffraction into practice for 
these measurements. Using the high neutron flux 
made possible by the heavy water moderated 
uranium chain reactor, Zinn” has been able to 
demonstrate the diffraction of neutrons from a 
single crystal and to show the application of the 
phenomenon to cross section measurements. The 
crystal spectrometer for neutrons, unique in that 
it produces monoenergetic beams of low energy 
neutrons, provides a relatively simple tool for 
continuous cross-section measurement by activa- 
tion or transmission for energies up to 65 ev with 
present collimators. 


II. DESCRIPTION OF APPARATUS 


The instrument used for the measurements was 
the heavy duty analog of an optical spectrometer 
described by Zinn,” permanently mounted in the 
concrete floor of the pile room in line with port 
number 4 of the Argonne heavy water pile. In its 
final, accurately centered and leveled position, 
tests of the angular reproducibility of the arm 
showed that its extreme end reproduced within 
0.0005 inch repeated settings made on the dial 
plates of the instrument. For measurement of 
spectra and of cross section the table and arm 
could be coupled together by means of a friction 
clutch to drive the two elements in the two-to- 
one angular ratio required by these experiments. 

A 61-cm long boron trifluoride counter filled to 
a pressure of 40 cm of mercury, operated in the 
proportional region and mounted with its axis 
parallel to the direction of the diffracted beam, 
served as the neutron detector in the apparatus. 
Use of enriched boron, consisting of the single 
isotope B”, increased the detection efficiency for 
high energy neutrons approximately by a factor 

1! Rainwater and Havens, Phys. Rev. 70, 136 (1946) 
and Havens and Rainwater, Phys. Rev. 70, 154 (1946); 
Rainwater, Havens, Wu, and Dunning, Phys. Rev. 71, 65 
(1947); Havens, Wu, Rainwater, and Meaker, Phys. Rev. 
71, 165 (1947); Wu, Rainwater, and Havens, Phys. Rev. 


71, 174 (1947). : ’ 
2 Zinn, Phys. Rev. 70, 102 (A) (1946); Zinn, preceding 


paper. 
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Fic. 1. Schematic diagram of apparatus for total neutron 
cross-section measurements. 


of 5. Operated at 1500 volts above ground with 
electron collection on the central wire, the effi- 
ciency of detection of low energy neutrons ap- 
proached 100 percent. One stage of preamplifica- 
tion, together with a cathode follower was 
mounted on the boron-cadmium shield of the 
counter, from which point the signal was carried 
by means of an 18-foot cable to an amplifier and 
scaler on an electronic instrument panel near the 
apparatus. The normal neutron background 
counting rate with the pile running was about 15 
counts per minute, a small fraction of the 6400 
counts per minute of Bragg reflected neutrons 
from the (100) planes of LiF in the region of the 
Maxwellian maximum. Counting rates in the 
energy region above the Maxwellian portion were 
rably lower, since here the number of 
s in the spectrum as well as the detection 
y decreased. 

shown in Fig. 1, collimation of the pri- 
ary neutron beam from the tank wall was 
achieved by inserting a 9-foot channel into the 
shielding of the reactor. This channel, a steel- 
bounded orifice with a 1.272.54-cm aperture, 
was surrounded by a mixture of paraffin, iron 


_ shot, and boric oxide to confine the path of 


neutrons and gamma-rays. The resultant pri- 
mary neutron beam had a divergence of 32 
minutes, and the diffracted beam a base width of 
114 minutes. A cadmium diaphragm, consisting 
of a series of 12 long cadmium metal sheets 
spaced 3 mm apart, mounted on the spectrometer 
arm before the detector, served as a means of 
varying the divergence of the diffracted beam 
and the resolving power of the instrument. For 
the cross-section measurements, the half-width of 
the Bragg reflected beam of neutrons was 
standardized at 22 minutes. 


MEASUREMENT OF NEUTRON CROSS SECTIONS 


Crystals and Absorber Mounts 


Since the diffracted neutron beam for this 
series of experiments was obtained by reflection 
from the flat surface of a crystal rather than by 
transmission, the crystals employed were gener- 
ally mounted in a long array in order to increase 
the strength of the diffracted beam at small 
glancing angles. Aside from minor variations 
necessitated by the differing shapes of available 
crystals, mounts consisted simply of a metal 
frame supported vertically on the spectrometer 
table. With this standardized type of mounting, 
transfer from crystal to crystal was simplified by 
use of a guide bar on the crystal table, since the 
planes of the crystal were aligned parallel to the 
metal frame. The mounts were usually made of 
aluminum with no metal in the path of the pri- 
mary beam to minimize the incoherent scattering 
from the holder itself. 

It was possible to measure the transmission, 
and consequently the total cross section of an 
element, by introducing it into the primary 
beam of neutrons, or by placing it in the Bragg 
reflected monoenergetic beam. Because the dif- 
fracted beam was considerably less intense, 
samples were mounted in the path of this beam to 
avoid irradiating them excessively. In this posi- 
tion also, the number of neutrons scattered by 
the sample was lower than if it were in the pri- 
mary beam. 

Except for elements available in foil form, it 
was necessary to add a second substance, usually 
sulfur or a solvent, as a support for the sample 
while its transmission was being measured. This, 
however, made it necessary to consider the effect 
of the binder or solvent on the transmission 
measured. To do this, a second sample containing 
the binder alone was prepared, and the number of 
neutrons passed through this sample was meas- 
ured. The ratio of the counting rate with the 
bound sample in place to that with the binder 
alone in the beam was taken as an uncorrected 
measure of the transmission of the sample. A 
sliding holder with several apertures was used to 
support the two samples successively in the beam. 
Three apertures were provided in the holder to 
make possible measurement of two samples with 
identical supporting materials concurrently. 
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Monitors and Electronic Circuits 


Two monitors were regularly employed with 
the spectrometer in the measurement of neutron 
cross section. The first, a monitor of the neutron 
intensity in the primary beam, was a fission 
chamber mounted near the reactor tank wall, in 
a removable channel adjacent to that used as a 
primary source for the spectrometer. The position 
of the monitor with respect to the tank wall could 
be varied so that varying pile intensities yielded 
the same reading on the monitor. Although the 
pile was automatically controlled, and did not in 
general fluctuate by more than 3 percent, the 
monitor was useful in averaging any short term 
fluctuations and in taking into account any long 
term drift of the pile. The monitor was an alumi- 
num fission counter containing a thin foil of 
uranium oxide, U;Os, enriched some ten times in 
the isotope 235. This monitor, operated at 1000 
volts above ground, was used as a measure of 
integrated incident neutron flux in timing the 
measurements. 

A second monitor measured the variations in 
the intensity of neutron background which were 
caused by changes in neutron leakage from other 
points of the pile. This monitor, an unshielded 
boron trifluoride counter operated in the pro- 
portional region, was identical with the main 
detector described above, except that it was 
filled with normal boron trifluoride gas to a 
pressure of 40 cm of mercury. Since this monitor 
was 100 times more sensitive to change in back- 
ground neutron intensity than the main de- 
tector, it increased the reliability of cross-section 
measurements by decreasing the uncertainty in 
background variation. In general the background 
neutron intensity in the pile room was quite 
small, and this monitor was used to detect only 
great variations in background. However, rather 
than correct the data for these background 
variations, the reading which was indicated by 
this monitor to be disturbed was generally re- 
peated under more favorable conditions. 

Each of the three counters, the main detector, 
the fission monitor, and the background counter, 
had associated with it a preamplifier, cathode 
follower, amplifier, scaler, and high voltage 
supply. The preamplifiers and cathode followers 
were mounted as near as possible to the counter, 


J. STURM 


while all high voltage supplies, amplifiers, and 
scalers were mounted together on a control panel 
on a work bench near the instrument. The de. 
tector and background monitor signals were 
recorded on two scales of 64, while the fission 
signal was recorded on a scale of 512. 

Between the amplifier and the scaler of the 
main detector circuits an additional limiting cir. 
cuit has been provided,* the function of which 
was to eliminate large negative or positive pulses, 
and any small pulses with sign opposite to that of 
the neutron pulse. This circuit has proved usefyl 
in eliminating pulses caused by an electromag. 
netic disturbance in the vicinity of the spectrome- 
ter. For facility, and to assure simultaneous 
stopping and starting, all of the circuits were 
operated from a single switch. An oscilloscope 
was provided to watch the signal at various 
stages of amplification and scaling, and regular 
checks were made to assure the perfect function- 
ing of the instrument. During measurements of 
cross section the oscilloscope constantly regis- 
tered the output of the main detector. 


Ill. ALIGNMENT 


The condition for a maximum in the diffrac- 
tion pattern of a space lattice is given by the 
familiar Bragg law, 


n\ = 2d siné@, 


where d is the interplanar distance, m a 
order and wave-length respectively, and 
glancing angle. Expressed in terms of energy, 
relation becomes 

ev = kn?/sin?6 (2) 


where k is a constant for each interplanar spacing. 
The problem of alignment, then, was one of 


determining the glancing angle 0, and thus in- ’ 


volved centering the instrument and crystal 
table with respect to the primary beam, aligning 
the detector, and finally aligning the crystal. 


Alignment of the Table, Detector, and Crystal 


Two methods ‘have been used to center the 
crystal spectrometer and its table with respect 
to the center of the primary beam. In the 
first, a long, snug-fitting steel bar was inserted 


* This circuit was designed by Mr. Thomas Brill of this 
laboratory. 
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into the primary beam collimator and extended 
over the spectrometer table. With this method, 
by inversion of the bar to cancel out any lack of 

allelism, the center line of the slot could be 
located within } mm. 

The second method used for preliminary de- 
termination of the center of the beam involved 
the use of photographic film. In this method a 
sheet of orthochromatic film was placed between 
two layers of cadmium metal, and this sandwich 
exposed to the beam of neutrons. A small hole 
drilled through the cadmium sheets was used to 
orient the film with respect to the crystal table, 
and the subsequent blackening of the emulsion 
indicated the shape and position of the neutron 
beam. The clear spot in the darkened area, repre- 
senting the hole in the cadmium layers, served 
then to relate the position of the crystal table to 
that of the beam. This method, although less 
accurate than the foregoing because of the asso- 
ciated flux of gamma-rays which caused about 15 
percent of the blackening, was useful because it 
showed both shape and rough position of the 
primary beam at any point on the arm. 

After the detector had been carefully aligned 
parallel to the arm of the spectrometer and the 
table centered with respect to the beam, it was 
necessary to locate the angular position of the 
arm at which the arm was parallel with and 
centered in the primary beam. In order to center 
the detector with respect to the primary beam a 
curve of neutron intensity as a function of 
angular position of the detector was plotted. 
Since the primary beam consisted of a very in- 
tense flux of neutrons and gamma-rays, it was 
necessary to reduce its intensity so that this 
measurement could be made. In order to reduce 
the neutron intensity, a thin slab of paraffin was 
introduced directly into the beam, and, to elimi- 
nate entirely the effect of the gamma-rays, the 
gain of the amplifier was changed so that only the 
higher neutron pulses were counted. The zero 
position of the detector was then determined to 
be the point of intersection of two lines which 
represent the slopes of the two sides of the curve 
taken in this manner. This position could easily 

be determined within 2 minutes of arc, and thus 
contributed an uncertainty of but 1 minute of arc 
in the glacing angle @. 
Having determined the position at which the 
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detector was in the center of and parallel with the 
direct beam, the detector was turned off this zero 
position by a known amount, usually 40 degrees, 
and the cadmium diaphragm and crystal care- 
fully placed in position on the instrument. As the 
crystal approached an angle of 20 degrees with 
the primary beam, the diffracted maximum could 
be detected. Leaving the arm in its fixed position, 
the crystal was rotated for several degrees on 
either side of this maximum, and a plot of the 
measurement of the intensity of diffracted neu- 
trons as function of glancing angle made. A 
lack of symmetry about the maximum of the 
observed diffraction curve arose from the fact 
that it was taken with a boron detector, since 
boron follows the well-known 1/v law for absorp- 
tion of neutrons and is not uniformly sensitive 
over the entire spectrum range. These data, 
corrected for the 1/v effect, determined the peak 
of the diffraction maximum within two minutes 
of arc. 

The uncertainty in the glancing angle @, caused 
by uncertainties in centering of the crystal table, 
could be taken as negligible because the width 
of the beam and the width of the counter were 
very nearly alike and because of the precision 
with which this position was ascertained. Be- 
cause the index of refraction of the crystal used 
for neutrons is not equal to unity, as has been 
shown for many substances, one might expect 
this to introduce an additional error in the 
measurement of the energy of the diffracted 
beam. However, since the value of the index of 
refraction is very likely of the order of 110~-*, 
this effect too was completely negligible. Thus 
the total errors were introduced in aligning the 
detector, and in aligning the crystal, and resulted 
in an uncertainty of three minutes of the glancing 
angle throughout the range of the instrument. 


IV. CRYSTALS 


Although many crystals, including mica, cop- 
per, calcium fluoride, calcite, and lithium fluo- 


TABLE I. Relative reflectivity from several crystals. 











Maximum Relative 

Crystal kX108 energy (ev) reflectivity 
LiF (100) 5.08 65.0 6.4 
Calcite (100) 2.22 28.0 2.5 
Mica (001) 0.20 2.5 9.2 
Cu (100) 1.56 20.0 ~5.0 
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ride have been used and tested for crystal spec- 
trometry, all of the measurements of cross section 
as a function of energy were made with the 
lithium fluoride crystals. The spectrum from the 
side hole of the heavy water pile has been 
measured with several crystals, and the relative 
peak intensities observed with each of these 
crystals have been tabulated in column 4 of 
Table I. 

As in the case of x-rays, the reflectivity of a 
crystal for neutrons can be increased by increas- 
ing the mosaic structure of the surface. A single 
crystal of LiF(100) has been carried through 
several steps, each involving a change in the 
surface of the crystal. The relative reflectivity 
from these surfaces is shown in Table II. In 
each case the rocking curve was 22+1 minutes 
in half-width, and since this value was deter- 
mined by the collimating geometry, no quanti- 
tative conclusion with regard the change in 
mosaic structure could be made from these data. 
However, obviously the roughened crystal, with 
its higher reflectivity, is best for use in spec- 
trometer measurements with this type of collima- 
tion. Since the range of energies covered by the 
spectrometer is a function of the interplanar 
spacing of the crystal used, a consideration which 
must be weighed in selecting a crystal for spec- 
trometer use is the value of d. As is obvious from 
Eqs. (1) and (2), the highest energy which can 
be reached by any crystal as a function of its 
lattice spacing varies as the square of the recip- 
rocal of the lattice spacing. Column 3 of Table I 
shows the maximum energy which can be reached 
using these crystals, assuming that the collima- 
tion is such as to permit the utilization of dif- 
fracted beams whose incident glancing angle is 
30 min. 


Calcite and Lithium Fluoride 


All early measurements on the thermal col- 
umn, and some of the first measurements at the 
side port of the Argonne heavy water pile, were 
made with a calcite crystal. This crystal, how- 
ever, because of its lattice perfection, showed a 
lower reflectivity than many of the other crystals 
tried. Since, in addition, it is quite a hard crystal 
and increased mosaic structure cannot easily be 
introduced, this crystal was used only for pre- 
liminary work. 
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TABLE II. Variation of reflectivity with crystal] 
surface. LiF (100). 








Surface 





Smooth cleavage face 
Surface ground 
Surface polished 


i 
Reflectivity 
as 
1.00 
1.65 
34 
Surface roughened with coarse abrasive 2.36 








The simple cubic crystal, lithium fluoride, was 
the crystal selected for most of the work reported 
in this paper. Because of the small size of the 
unit cell, and because of the consequent small 
interplanar spacing in this crystal, the range of 
energies available for this series of neutron spec- 
trometer studies extended up to 65 ev. These 
crystals, quite easily available, had the additional 
advantage that their mosaic structure could be 
increased simply by the various techniques dis- 
cussed above. 

The 100 planes of this crystal were used for all 
of the measurements of cross sections reported 
here. The unit cell dimension for LiF is 4.01 ang. 
stroms, and the interplanar spacing is 2.005 ang- 
stroms. With this crystal, the range of useful 
energy lay between 0.04 ev and 65.0 ev, a factor 
of 1.6X10* in energy. The low energy limit was 
set by choosing that energy characterized by the 
presence of less than 5 percent second and higher 
order components in the Bragg-reflected beam, 
as determined by a series of experiments of the 
type to be described; in the case of the high 
energy extreme, the limit was fixed by the 
smallest useful glancing angle which could be 
achieved with reasonable collimation, in this 
case 30 minutes. The LiF crystals used, a bank 
of small individual crystals mounted in an alu- 
minum holder to form a long array, consisted of 
5 small crystals whose individual dimensions 
were 3.7 X2.7 X.2 cm. The planes of each of the 
individual crystals were aligned within 1 minute 
of arc. 


V. SPECTRUM OF THE SOURCE 


Two ports of the Argonne heavy water chain 
reacting pile were used for all of the measure- 
ments. The first, the thermal column, interposed 
a thickness of carbon between the source and 
the spectrometer rendering the greater part of 
the neutrons thermal. Their distribution, roughly 
Maxwellian, was such as to provide a useful 
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source of neutrons in the energy range below 
0.2 ev. The second port, a side hole of the same 
reactor, interposed little moderator between the 
source and the spectrometer, and while a strong 
Maxwellian component was present in this beam 
because of the moderator inside the reactor, this 
beam contained, in addition, a fairly strong com- 
ponent of higher energy neutrons which rendered 
the extension of the range of the instrument to 


65.0 ev possible. 
The Maxwellian Region 


Attempts to explain the results observed in 
the measurements of spectra from the thermal 
column and side hole of the reactor, have been 
made by Seitz and Goldberger.** The form of the 
equation relating the crystal reflectivity R with 
the glancing angle 0, for detection by a 1/v 
absorber as developed by them is 


R,=Clan|nK;, csc*0 exp [—a,? csc?6]. (3) 


In this equation m is the order of the reflection, 
C is an arbitrary constant, and 


K,=22/d 


 dn=LstLo,'(b) } exp [2ri(hixj+hay ;+hse;)n] 


h(2x/d,)? 
= ———_ 7? 
8mkoT 


a,” 


in which d, is the interplanar distance, m is mass 
of the neutron, ¢,’(b) is the scattering cross sec- 
tion of the j’th atom in the unit cell, o is Boltz- 
mann’s constant, and 7 is the average neutron 
temperature. 

The fit of this expression for T=400°K to the 
experimental points in the Maxwellian region is 
shown plotted as a solid line in the curve of 
Fig. 2. The distribution at higher energies, when 
corrected for detector efficiency, resolution, and 
amount of beam intercepted, corresponds to the 


‘expected 1/E function. 


VI. ORDER EFFECTS 
A study of the intensities of neutron diffrac- 


tion in the various orders by a series of crystals 


rode L. Goldberger and F. Seitz, Phys. Rev. 71, 294 
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ELECTRON VOLTS 


Fic. 2. Spectrum of reactor as measured with LiF (100) 
crystals. The solid line represents the theoretical Max- 
wellian distribution for T=400°K. 


has recently been made by Fermi and Marshall.” 
Their results show that the relative amounts of 
first and higher order Bragg reflection varies 
greatly from crystal to crystal and from plane 
to plane, depending on the scattering cross sec- 
tions and the phase of neutron scattering by the 
target nuclei. Thus the choice of a crystal for 
cross-section measurements determines the useful 
range over which data can be taken with dif- 
fracted beams uncontaminated by higher order 
components. 

From Eq. (2) we see that for a constant 
glancing angle, 0, the energies of the first and 
second order have a ratio 1 to 4 in energy. A 
series of experiments were performed to evaluate 
approximately the fractional composition of the 
beam quantitatively to determine the useful 
range of the spectrometer for measurements of 
cross section. 


18 Fermi and L. Marshall, Phys. Rev. 70, 103 (A) (1946). 
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Experimental Determination of the Order Com- 
position of the Diffracted Beam 


Since total cross sections were calculated from 
measurement of transmission, the effect of a 
beam composed of several energy components 
will alter the observed transmission of a sample 
as a function of both quantity of flux having 
each energy and of the cross section of the sub- 
stance at each energy. A series of measurements 
was made on the samarium sample discussed 
below, the cross section of which is shown in 
Fig. 13. This sample has a strong resonance 
level at 0.1 ev, and detailed measurement of its 
cross section shows two artificial peaks at 
0.025 ev and at 0.011 ev, corresponding to ab- 
sorption of large fractions of the second and 
third order components at these energies by the 
primary resonance at 0.1 ev. 

Although the position of the two peaks in 
cross section at low energy is exactly that ex- 
pected for absorption of second and third-order 
components from the diffracted beam (0.025 x4 
=0.10 and 0.011 X9=0.099 ev), a preliminary 
experiment was performed to show that the two 
observed peaks at low energy are not real reso- 
nance levels. A very large part of the 0.1-ev com- 
ponent of the primary neutron beam at 0.025 ev 
was removed by filtering this beam through a 
second samarium solution, the transmission of 
which, for the high order component was about 
0.10. A measurement of the transmission of the 
first solution was made with and without the 
filter in place. Without the filter the transmission 
was 0.441, with the filter the transmission meas- 
ured 0.814, indicating a lower absorption because 
the strongly absorbed 0.1-ev component was 
minimized in this filtered beam. 

In order to get a quantitative estimate of the 
amount of second-order component in the beam 
diffracted by LiF(100) at 0.025 ev, a measure- 
ment of a standardized thin film of boron was 
made at this energy. The boron layer was a 
plate of Pyrex, the transmission’ of which had 
been previously measured in a mechanical veloc- 
ity selector in which order effects were absent. 
Measured in this way, the transmissions of this 
. plate for the first-order (0.025 ev) component 

- was 0.441 and for the second-order (0.1 ev) com- 
ponent was 0.663. A measurement on the crystal 
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spectrometer yielded a transmission of 0,490, 
intermediate between the values for the first and 
second orders alone. 

It is clear that in this measurement with a 
beam of several energy components impinging 
on an absorber of energy-dependent cross section, 
if A, is the fractional order component, and T, 
is the transmission of the absorber for the energy 
of the order n, then 


AoT9=AiT\+A2T24As373::: (4) 


where 7) is the observed transmission and A, 
is the total flux and is equal to one. Since there 
are but two orders of significance at this angle, 
substitution of the transmission values given 
above and combining with the knowledge that 


Ayp=Ai+A2=1, (5) 


showed that, for first order equal to 0.025 ev, 
the fraction of first and second-order components 
was: A,;=0.78; A2=0.22. 

A second evaluation of the composition of the 
diffracted beam at the same energy was made 
by altering the composition of the beam as in 
the case above, with a filter having known ab- 
sorption for the two energy components.‘ The 
filter for this measurement was a second and 
thicker samarium sample, having very low trans- 
mission for the second-order (0.10 ev) energy 
component, thus lessening the fractional absorp- 
tion of this strongly absorbable component at 
0.025 ev by the sample. A schematic representa- 
tion of the transmission of this series of absorbers 
is shown in Fig. 3. The primary beam of intensity 
Ao is considered to be divided into its first and 
second-order components A; and Az, and the 
lessening of their intensities is shown schemat- 
ically in the figure. The observed transmission 
of the samarium absorber in the filtered beam is 
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given by 


TiAi1Tpt+T2A2T ye 
To= (6) 
AiTpt+A2T pr 





where 7,=the transmission of the absorber for 
neutrons of the energy of the order n, 7;, = trans- 
mission of the filter for neutrons of the energy of 
the order n, A,=the intensity of the order n. 
Values of 7, and 7;, were obtained by calcula- 
tion from the quantity of absorber in each cell 
and from the cross section shown in the Breit- 
Wigner curve plotted in Fig. 13. The values 
thus obtained were: 7;,=0.374, Ty2.=0.090, 
T)=0.593, 7, =0.613 and: 7,=0.280. Again com- 
bining these by use of Eq. (5) we calculate that 
A,=0.78 and A2=0.22. 

As a third method of evaluation of the ratio 
of order components at 0.025 ev, a filter measure- 
ment, similar to the last one was made. In this 
case the filter used was the same standardized 
Pyrex plate of the earlier measurement, and the 
principle the same as that of the previous experi- 
ment. Values of 7, were calculated from Fig. 13 
and those of 7;, from the standardization men- 
tioned earlier. The values used in Eq. (6) were: 
Ty, =0.441, Ty2=0.663, T1=0.436, T.=0.294, 
and 7)=0.400. The calculated fractions of first 
and second order in this case were: A,=0.79 
and A:=0.21. The average value of the first- 
order component of the reflected beam, 78 per- 
cent, agrees quite well with the theoretical value 
at this energy calculated from Eq. (3). 

The theoretical order intensity curve calcu- 
lated from Eq. (3) which was used for all cross- 
section corrections is shown in Fig. 4. 


Reduction of the Effect of Higher 
Order Components 


In an effort to remove the effect of high order 
components in the diffracted beam at large 
glancing angles several techniques were used. 
Since the total cross sections for neutron reac- 
tions in this energy region in general show a 
decreasing cross section with increasing neu- 
tron velocity, cross sections in the region of 
hundredths of electron volts were measured using 
samples with transmissions somewhat greater 
than 0.5. In this manner, since the transmission 
of the higher energy components of the reflected 
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beam is generally much higher than that of the 
first order, the effect of the second-order neutrons 
on the observed cross section is minimized. 
Other means of reducing the effects of higher 
order components involved reduction of the in- 
tensity of these components in the primary beam, 
This was accomplished by use of filters with 
resonance levels at the energy of the second and 
higher order components, or by using the pri- 
mary neutron beam of lower mean energy and 
minimal 1/E component which is emitted from 
the thermal column. Attempts were made also 
to remove neutrons of the energy of the second, 
fourth, etc., orders, by introducing a second 
(mica) crystal in advance of the lithium fluoride 
in the primary beam at the angle necessary to 
diffract these components from the beam. These 
efforts were unsuccessful, since, aside from reflec- 
tivity considerations, the imperfection of even 
a mica crystal is not as great as was the geomet- 
rically determined divergence of the primary 
beam. These methods of beam purification were 
considered necessary only for the measurements 
made at energies less than 0.04 ev, below which 
energy the diffracted beam contained less than 
95 percent first-order neutrons. The small correc- 
tion necessary at energies lower than 0.08 ev 
was applied to the observed data at the end of 
the cross-section measurement by use of rela- 
tions (4) and (5). Since no order correction was 
made to data taken at energies above 0.08 ev, 
this correction was applied only to the few points 
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Fic. 4. Fractional order composition of diffracted beam as 
a function of glancing angle. 
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ELECTRON VOLTS 


Fic. 5. Resolution of the neutron spectrometer calculated 
from the width of the rocking curve. 


measured at very low energies, and had little 
effect on the final results. 


VII. RESOLUTION 


The angular divergence of the diffracted beam, 
as obtained from measurements of the crystal 
rocking curve, can be taken as a first approxima- 
tion and upper limit to the energy spread of 
that beam. The crystal rocking curve, as experi- 
mentally measured for glancing angles between 
45° and 33°, showed uniform shape and half- 
width; and although the shape of the rocking 
curve was not exactly triangular because of the 
1/v detector used, the shape of the energy dis- 
tribution function was assumed to be an ideal 
triangle for calculation of the resoliition function 
and for resolution correction. The base of these 
triangles indicated an energy spread correspond- 
ing to thé base of the rocking curve, 44 minutes, 
for energies up to 1.4 ev, and the corresponding 
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values of E/5E vary between 28 at 0.01 ey and 
somewhat less than 1 at 65.0 ev. (See Fig. 5.) 

The spectrometer resolution as a function of 
neutron energy varies as 1/d, the extent of the 
useful energy range varies as 1/d*. Thus the 
cleavage (100) face of LiF because of its small ¢, 
provided a wide range of energies of sufficiently 
well-resolved monoenergetic neutrons for crogg. 
section measurements. In practice, the resolving 
power of the instrument was determined by the 
minimum width of collimating slits which pro- 
vided ample diffracted intensity at high energies, 
This width, 1.25 cm, could have been made 
smaller at low energies where the intensity was 
higher, but this was not considered necessary, 
as the resolution was more than ample with the 
wide slit in this region, and under these condi- 
tions the intensity was considerably increased, 

The base width of the rocking curve was not 
constant for energies greater than 1.4 ev with 
LiF(100), since for glancing angles less than 
3° 28’ and energies greater than 1.4 ev, the 18-cm 
bank of crystals was not sufficiently long to 
intercept the entire primary beam. As a con- 
sequence, the Bragg-reflected beam at these small 
glancing angles is less divergent and the width 
of the rocking curve smaller, varying from 44 
base width for a 3° 28’ glancing angle to 19 at 
30’ glancing angle. This resulted in an improve- 
ment in resolution and a loss of intensity in the 
high energy region. The triangles shown on the 
curves of cross section indicate the idealized 
shape of the final intensity-energy distribution 
of the diffracted neutron beam at the energy of 
the apex. 


Resolution Correction 


It was possible to correct the transmission 
measurements in the energy region below 0.2 ev 
for the slight distortion in the nature of this 
curve caused by the finite resolution.* To do 
this, we note (see Fig. 6) that the measured 
transmission, ®, is given by is 


E'+y 
f f(|E—E’|)6(E) N(E)dE 
On (1) 





E’ty 
f f(|E—E’|)N(E)dE 
E’—¥ 


4M. Goldberger, unpublished work. 
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where f(|£—£’ |) is the resolution function, ¢(£) 


_ jg the true transmission at energy £, and N(£) 


is the incident flux distribution over the energy 
interval. Assuming that N(Z) in the region of 
interest is very nearly constant, we can substitute 
in the simplified expression the triangular reso- 
lution function (Fig. 6) 


f(\E-E’|)=(8/y((IE-EP+7} (8) 


we obtain 


voce) = f ” B6(EME- f 


E’-—y E’ 


EB’+y7 


Eo(E)dE 


E'’+y7 
+ f (E'+7)6(E)dE 
E 


a f (E’—y)o(E)dE. 
E’—y 


Expanding the right side of the equation about 
7=0, 


2 4 


Y Y 
6(E’) =9(E’) +,-¢"(E')+—o"(E’).._ (9) 
12 360 


To express ¢, the true transmission, in terms of 
operations on ®, the measured transmission, we 
neglect in the first approximation the term in- 
volving the fourth derivative. Solving the sim- 
plified last equation for ¢, we find 


« (—1)*b?(Z) 
¢(Z)= > 


K=0 atk 





in which (2k) refers to the 2kth derivative, and 
a=[12/y*]!, whence to this degree of approxima- 
tion, 

o(E) = #(E) —(7?/12)6"(E). (10) 


By perturbation theory, taking into account the 
term in the fourth derivative, 


$(E) = (E) — (7°/12)®”(E) 
— (7*/360)8""(E), (11) 


the final form of the relation, valid where y is 
small relative to the half-width of the resonance 
level. 


Vill. MEASUREMENT OF TOTAL NEUTRON 
CROSS SECTIONS 


Cross sections can be measured by transmis- 
sion or by activation in the monoenergetic dif- 
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fracted neutron beam of the crystal spectrometer. 
Because, however, the intensity of the Bragg 
beam does not exceed 30,000 neutrons per min- 
ute. Activation measurements are thus limited 
by the low neutron flux to elements with large 
cross sections and short periods, and that only in 
the energy range of high flux near the peak of the 
Maxwellian distribution. Results reported here 
were made by measuring transmissions of a small 
sample of the absorbing substance by means of 
a boron counter. 

Total cross sections are calculated from the 
measured transmission, 7; from the relation 


T=e"", (12) 


in which ¢o is the total neutron cross section in 
square centimeters per atom and N is the surface 
density of atoms in the sample. Solving for o, 
we have 


ao =(1/N) In (1/T) (13) 


which is the form of the relation used for the 
calculation of ¢. The value of N was obtained by 
simple assay of the foil, powder, or solution con- 
taining the element of interest, and T was ob- 
tained from the measurement of the intensity 
of the diffracted beam with and without the 
sample in place. In general, the thickness of the 
sample was made such that its transmission was 
somewhat over 0.5 for reasons indicated above. 
It was not necessary to apply a correction to the 
transmission data for beam hardening because 
of the fact that the diffracted beam consisted 
almost entirely of neutrons of a single energy. 

It was, however, necessary to correct both the 
numerator and denominator of the transmission 
measurement for a background of neutrons in- 
coherently scattered by the crystal. Since, in 





















































addition, the absorbing sample was interposed 
between the crystal and the detector, the back- 
ground varied when the sample was placed in the 
diffracted beam. Because of this variation, to- 
gether with variations in background caused by 
the differing positions taken by the crystal and 
detector, it was found necessary to measure this 
background counting rate for both the directly 
diffracted beam and the diffracted beam with 
sample interposed at each energy. The back- 
ground measurements were taken by turning the 
crystal one degree off the Bragg maximum, an 
angle at which only the incoherently scattered 
neutrons reached the detector. The transmission 
of the sample was then calculated from the four 
measurements according to the relation 


T=(Ri—n)/(R2—12) (14) 


in which R; is the counting rate of the beam 
transmitted through the sample, R, the counting 
rate of the open diffracted beam, and 7 and rz 
the background rates with and without the 


sample interposed. 
Samples 


Samples for cross-section measurement were 
prepared in one of four ways. For rhodium, gold, 
and iridium varying thicknesses of metal foils 
of these substances were used. Part of a pre- 
liminary boron curve was taken with pellets of 
amorphous boron powder pressed into the form 
of cylindrical disks of varying thicknesses and 
one inch in diameter. A third method, involving 
the use of a sulfur binder in the pressed pellets, 
was used for samples with high total cross sec- 
tions (part of the iridium curve), where it was 
impossible to get samples sufficiently thin other- 
wise. Transmission measurements in these cases 
involved the use of very fine powders of the 
substance, and necessitated the use of a second 
blank pellet of sulfur of weight equal to that of 
the sulfur binder in the first disk. 

The method most frequently used was a solu- 
tion method, which was used particularly for 
the rare earths gadolinium, samarium, europium, 
and dysprosium, where cross sections were 
usually high and samples were small. For this 
purpose two very precisely ground quartz cells, 
1 inch in diameter and 1 centimeter thick were 
employed. The nitrate of the rare earth was 
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prepared, a proper amount dissolved in a known 
amount of heavy water and this solution added 
to one cell. Heavy water was used as the solvent 
in order to avoid the loss in beam intensity dye 
to the strong scattering of hydrogen-containing 
solvents. To the heavy water in the second cel] 
was added a sufficient amount of DNO; to equal 
the effect of the nitrate ion in the rare earth cel], 
Measurements were made then as in the case of 
the sulfur bound pellets. 


Evaluation of Approximate Theoretical Constants 
for the Resonance Levels 


Although the experimental data are not suff- 
ciently precise to make very definite assignments 
to the constants in the Breit-Wigner relation, and 
thus cannot be considered either a test of the 
theory or the last word on the exact shape of the 
cross-section curves, it was considered of interest 
to evaluate the constants of the theoretical ex. 
pression which best fitted the experimental 
curves. The form of the Breit-Wigner one-level 
equation for the energy dependence of the cross 
section is 


ooEo! 
~14+[(E-E,)/3r}? 


in which Ep is the energy of the resonance level, 
I is the width of the level at half-maximum, » 
is the cross section at the resonance maximum. 
The values of oo, Ey) and I can be assigned from 
the observed data, and are included in the discus- 
sion of cross sections below where their values can 
be considered to approximate the correct values 
because of the quality of the experimental data. 

In general, the observed data taken at neutron 
energies less than one electron volt, when cor- 
rected for the effect of higher orders and resolu- 
tion as discussed above, were considered suffi- 
ciently precise to justify evaluation of these 
constants. Data at higher energies are seriously 
disturbed by lack of resolution, and although the 
energy of the peak is probably accurate, the 
measured widths of the levels are in general 
much too broad, while the peaks are for the 
same reasons very much flattened. 


cE} 





IX. RESULTS 


The total, absorption plus scattering cross sec- 
tions for neutron reaction .have been measured 
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of a series of well-known strong absorbers, de- 
tectors, and rare earths using the method dis- 


cussed. All the cross sections reported have been. 


calculated from measurements of the transmis- 
sion of a sample of the substance over the range 
of energy indicated. Although isolated measure- 
ments of the activation cross sections have been 
made on the same instrument, notably in the 
case of indium, dysprosium, and europium, these 
results were not nearly as complete nor as accu- 
rate. 

In the cross-section curves shown in Figs. 7 
to 16 corrections to the observed data have been 
made for the variations in background caused 
by change in the amount of scattering as dis- 
cussed above. In addition, these data are cor- 
rected for resolution and for order composition 
in the ranges of low energy where the methods 
of correction discussed are valid. An assay of 
the chemical composition of each sample was 
made, and corrections for the energy dependent 
effects of highly absorbing impurities upon the 
measured transmission at each energy were made. 
In the range of energies higher than 1.0 ev, only 
background and impurity corrections were ap- 
plied and the cross-section values in the maxima 


Fic. 7. Total cross section of rhodium. The cross section, ¢, is given in units of 10~* cm*/atom (barns). 
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are lower limits only. The effect of the poorer 
resolution at these energies also widened these 
peaks, and the half-widths taken from the curves 
are maximum values. None of the cross sections 
have been reduced to isotopic cross sections or 
cross sections of the absorbing isotope, but all 
have been reported as cross sections per average 
atom. Data taken at energies greater than 30 
ev, where the detector was very near to the 
primary beam, can be considered to give only 
an indication of the nature of the cross-section 
function in this region. 


Rhodium 


Attempts to evaluate the energy-dependent 
variations in the cross section of rhodium have 
been made using the boron method* and by 
varying the mean energy of the reaction neutron 
source by changing the moderator temperature.’ 


‘Further, self and resonance absorption tech- 





niques have been used to obtain the width and 
relative position of the resonance level.'*-'* By 


( oon H. Goldsmith and F. Rasetti, Phys. Rev. 50, 328 
1936). 

16H. Feeney, C. Lapointe, and F. Rasetti, Phys. Rev. 
61, 469 (1942). 
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Fic. 8. Total cross section of gold, ¢ in units of 10“ cm?/atom. 


use of a series of six metal foils of the element 
augmented by an assortment of sulfur-bound 
pellets of the finely powdered element, it has 
been possible to measure this cross section in 
the continuous energy range from 0.04 ev to 
65 ev with the crystal spectrometer (Fig 7). 

It has been possible to place the resonance 
peak, the existence of which was predicted by 
the self-absorption and activation measurements, 
at 1.28 ev. However, because of the insufficient 
resolution of the spectrometer in this region, the 
observed maximum in the cross section repre- 
sents only a lower limit, and since, for the same 
reason, the measured width is not correct, no 
attempt was made to evaluate theoretical con- 
stants. At the high energy end of the curve 
evidence is indicated for some increase in cross 
section at energy greater than 30 ev. 


17J. H. Manley, H. H. Goldsmith, and J. Schwinger, 
Phys. Rev. 55, 39 (1939). 

18 J. Hornbostel, H. H. Goldsmith, and J. H. Manley, 
Phys. Rev. 58, 18 (1940). 

18 See also spectrometer measurements of L. B. Borst, 
A. J. Ulrich, C. L. Osborn, and B. Hasbrouck, Phys. Rev. 
70, 557 (1946); Phys. Rev. 70, 108 (A) (1946). 





Combining this result for the energy of the 
level with a cyclotron measurement of the posi- 
tion of indium level (1.44 ev) recently made by 
Havens and Rainwater," we find the value of 0.16 
ev for the energy difference between the peaks. 
This is in very good agreement with the results 
of Hornbostel, Goldsmith, and Manley who by 
mutual absorption methods found | Zora —Eora| 
=0.15 ev. 


Gold 


Two gold plates, which alone and in combina- 
tion provided three thicknesses of absorber for 
the various cross-section regions, were used for 
the measurements. The surface densities of these 


samples, 1.674 g/cm’, 2.898 g/cm* and 4.572. 


g/cm?, made it possible to measure the cross 
section of gold in the region between 0.04 and 
8.0 ev. The variation (Fig. 8) below 2.5 ev is 
very nearly as 1/v, and the cross section at kT 
was found to be 99.3 10-** cm?/atom, in close 
agreement with other measurements of the cross 
section at this energy. 
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?'W. Stupp, Ann. d. Physik 43, 630 (1943). 


Earlier measurements of the position of the 
resonance level by absorption methods'® and 
by modulated cyclotron beam™ placed it at 
2.6 ev and 4.8 ev, respectively. Crystal data 
indicate the position of the level to be at 5.4 ev. 
The more detailed recent studies using the 
time-of-flight technique, however, indicate this 
value to be somewhat high. The resolving 
power of the crystal, very nearly the same 
as that of Rainwater and Havens’ earlier 
apparatus in this region, results in agreement 
with regard to the apparent width and height of 
the level as measured by the two methods. 
(T'spparont = 2.48 ev, cyclotron, = 2.50 ev, crystal ; 
Soapparent = 198 X 10-* cm*/atom, cyclotron, =211 
x 10-* cm*/atom, crystal.) 
Iridium 

Early measurements which Amaldi and Fermi 
made of the absorption of cadmium-filtered neu- 
trons showed strong A group neutron absorption 
by iridium (~5 ev), and later evidence from 


boron absorption data of resonance absorption 
in iridium was obtained by Goldsmith and 
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Rasetti,* and the energy of this level was placed 
by them at 1.16 ev. 

Spectrometer measurements (Fig. 9) of the 
total cross section of iridium, at energies up to 
8.0 ev, show three absorption resonances, one 
at 0.635 ev, a second at 1.35 ev, and the third 
at 6.0 ev. Samples used were foils of the metal 
1.09, 2.22, and 3.29 g/cm?. Although the resolu- 
tion of the instrument is not sufficient to make 
very precise the determination of the theoretical 
constants, approximate values have been as- 
signed for the two low energy peaks. The two 
sets of constants for the one-level expression 
when added together account very well for 
the observed points. The two sets are: (1) 
Ey=0.635 ev, oo(Eo)'=425, T=0.250 ev; (2) 
Eo=1.35 ev, oo(E»)'=710, T =0.36 ev. 


Gadolinium 


The rare earth, gadolinium, long known to 
have a very large cross section for kT neutrons* 
was measured in heavy water solution. Three 
concentrations of the nitrate were needed to 
cover the various energy ranges measured, 
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Fic. 9. Total cross section of iridium. Cross-section values are multiplied by E} to show 
the 1/v region as a horizontal line. 
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Fic. 10. Total cross section of gadolinium, ¢ in 
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Fic. 11. ¢£*X10 plot of gadolinium data showing resonance 
level symmetrical about Ep. 
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Fic. 12. Total cross section of samarium, ¢ in units of 10-“ cm*/atom. 
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Fic. 13. Samarium cross section for low neutron energies. 

Two false peaks (shown on expanded scale in insert) 

hance correspond to absorption of second and third-order neu- 
trons by the strong resonance level at 0.096 ev. 
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0.0108, 0.00624, and 0.000867-g gadolinium ni- 
trate per square centimeter. The cross section 
(Fig. 10) was found to rise rapidly and con- 
tinuously as the energy of the reacting neutrons 
decreased, and the rate of rise was greater than 
that expected for a 1/v absorber. Plotting the 
data in a different manner (Fig. 11) serves to 
show the resonance level which is superimposed 
on a strong 1/v component of the total cross 
section. Since in this case, with o(EZ)! and E as 
the ordinate and abscissa, respectively, a 1/v 
absorber appears as a horizontal straight line, 
the resonance level is shown as a symmetrical 
peak superimposed upon it. Thus the level in 
this element is shown to be at very nearly 
thermal energy which accounts for its very large 
kT cross section. The Breit-Wigner constants 
for these data are: oo = 44,000 10-* cm*/atom, 
E,=0.031 ev, lf =0.095 ev. 


Samarium 


Measurement of the total cross section of 
samarium made in the energy range between 
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Fic. 14. Total cross section of europium. Slope of curve between 0.035 and 0.10 ev steeper than 1/». 
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0.04 ev and 65 ev, showed the presence of strong 
absorption levels at 0.096 ev, 10 ev, and 33 ev 
and a possible resonance at higher energy (Fig. 
12). As with gadolinium, several heavy water 
solutions of the nitrate and blank cells of dilute 
DNO; were used for the transmission measure- 
ments in the various energy regions. An attempt 
to fit the data theoretically was made only for 
the case of the resonance at lowest energy, its 
constants being: o=1550010-* cm?/atom, 
Ey =0.096 ev, and !' =0.074 ev. Previous indica- 
tion of structure in this cross section at low 
energies was obtained from absorption data and 
showed the existence of the lowest level, although 
no definite energy value could be assigned.” 











A plot of the uncorrected cross-section data at 
energies less than 0.04 ev shows strikingly the 
effect which the high energy, high order com- 
ponents in the diffracted beam can have on the 
observed cross section. The true resonance at 
0.096 ev reappears as two false subsidiary max- 
ima at 0.096/4 ev, when the second-order com- 
ponent of the beam is strongly absorbed, and 
again at 0.096/9 ev when the third-order com- 
ponent is strongly absorbed. The data are shown 
on an expanded energy scale in the insert of 
Fig. 13, indicating that these energies are just 
the expected values for the absorption of second 
and third-order components of the beam by the 
primary resonance at 0.096 ev. 
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Europium 


Measurements of the transmission of a sam- 
ple of very pure europium were made with the 
crystal spectrometer using the solvent method. 
Nitrates containing deuteryl water of crystalliza- 
tion were prepared and dissolved in heavy water, 
preparing thus samples containing 7.42 x 10” and 
2.1610" atoms Eu per square centimeter. A 
survey of this element for neutron energies (Fig. 
14) to 65 ev, showed the presence of resonance 
levels at 0.465 ev, 3.3 ev, 9.2 ev, and 22 ev, and 
possible levels at higher energy. Since in this case 
the maximum cross section observed in the peak 
of the resonance at 0.465 ev represents only a 
lower limit, assignment of a value to ap is neces- 
sarily low. The constants best fitting the ob- 
served points are: o9=5670X10-* cm*/atom, 
E,=0.465 ev, [=0.20 ev. The 3.3-ev level is 
probably due to capture by Eu’ which Gold- 
haber” has shown from boron absorption meas- 
urements to have a level at about 4 ev. From 


self-absorption measurements this level shows 
evidence of being somewhat wider than the usual 
levels measured, and appears here to have char- 
acteristics not inconsistent with this point of 
view. It is, however, not possible to distinguish, 
by present spectrometer measurements at this 
energy between a wide resonance and a group of 
narrow ones with close spacing. 

Crystal spectrometer measurements shows a 
rise in the cross section at energies less than 
0.15 ev at a rate greater than that expected for 
a 1/v absorber. It was considered that this region, 
although outside the range of the crystal spec- 
trometer, merited special emphasis. As a con- 
sequence, measurements of cross section in the 
energy range between 0.007 and 0.2 ev were 
made with a mechanical velocity selector which 
had this range.*** This instrument, the rotary 
shutter of which periodically releases bursts of 
polychromatic neutrons and by suitable delay 
mechanism counts them electronically at meas- 
ured intervals after the release, was able to meas- 
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Fic. 16. Total cross section of dysprosium, cross-section values multiplied by Z!. Negative slope 
indicates resonance level at energy less than zero. 


*° M. Goldhaber, personal communication. 


*** The assistance of G. Arnold, H. Kanner, and R. Nobles in making these measurements with the mechanical 


velocity selector is gratefully acknowledged. 

























776 WILLIAM 
ure transmission as a function of neutron energy. 
These data at low energy, combined with those 
from the crystal are shown plotted in Fig. 15. 
The constants of this low energy level were found 
to be: o9=5570X10-** cm?/atom; E)= —0.011 
ev; '=0.081 ev. 


Dysprosium 


Although the limited supply of dysprosium 
available at the time of the measurements 
shortened the energy range covered by the data 
for this element, sufficient evidence was obtained 
to indicate the existence of a neutron resonance 
in the cross-section curve. As with the other 
rare earths, quartz cells of the nitrate dissolved 
in D,O were used for the transmission measure- 
ments, the surface density of sample used here 
being 5.87 X10” atoms of dysprosium/cm’. 

The o(£)* curve of these data (Fig. 16) shows 
a continuously rising function with decreasing 
energy, which, by its slope, indicates a resonance 
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level at energy less than zero. The cross section 
at kT neutron energy was found to be 1180 « 19% 
cm? atom, a value considerably higher than an 
earlier measurement” of 870 X 10-* cm?/atom at 
this energy. The constants assigned to the the. 
oretical curve are of necessity highly uncertain 
because the only experimental data it is possible 
to obtain in this case is far from the peak of the 
level and represents but one wing of the reso- 
nance. The constants assigned to fit the data, 
assuming a half-width of 0.10 ev, are: —k, 
=1.01 ev, o. =79,000 X 10-* cm?/atom, I =0.10 
ev. 

The writer wishes to acknowledge his indebted. 
ness to Dr. Walter H. Zinn for his valuable 
suggestions and advice throughout the course of 
this investigation and to Mr. Solomon H. Turkel 
for his interest and assistance in many aspects 
of the work. 


21. Meitner, Arkiv. Mat. Astr. Fys. 27A, 3 (1941). 
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Diffraction Pattern of a Circular Aperture at Short Distances 


C. L. ANDREWS 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received January 6, 1947) 


The diffraction pattern of electromagnetic waves was studied at distances from zero to five 
wave-lengths from a circular aperture. Microwaves of 12.8-cm wave-length were employed. 
The aperture was an iris diaphragm that could be varied between one and six wave-lengths in 
diameter. The beam was incident normally upon the plane of the aperture from a 4 ft parabolic 
reflector 24 feet away. The sharpest diffraction patterns were in the plane of the aperture. 
Measurements were made of the intensity of radiation in the electric and magnetic planes 
through the axis. Individual plots were made of intensities along the diameter and along lines 
parallel to the diameter at steps of quarter wave-lengths from the diameter. Measurements 
were also made of the intensities at fixed points along the axis as the diameter of the iris was 
varied from one to six wave-lengths. It was observed that Fresnel zone theory could be em- 
ployed as a rule of thumb in predicting the intensity along the axis even in the illogical case 
of predicting the intensity at the center of the aperture. Checks were made against Kirchhoff's 


theory. 





INTRODUCTION 


HE wartime development of tools and 

techniques for the production of continu- 
ous microwaves has provided a supplement to 
the study of diffraction of x-rays and light. 

The study of the diffraction of electromagnetic 
waves by apertures has been of four kinds: 

1, Fraunhofer diffraction for which the source 
and the field to be studied are both effectively 
an infinite distance from the aperture. 

2. Fresnel diffraction for which the source and 
field to be studied are at large distances from 
the aperture compared to a wave-length but at 
small enough distances to require consideration 
of the effect of phase difference between second- 
ary wavelets from the different points in the 
aperture, even when the source and field to be 
considered are on a normal to the screen.! 

3. Diffraction by apertures of small dimen- 
sions compared to a wave-length.” 

4. Diffraction by apertures of the order of a 
wave-length in diameter observed in and near 
the aperture. 

When the second and third types of diffraction 
are treated theoretically, approximations can be 
made in the solution of Maxwell’s equations or 
in Kirchhoff’s semi-empirical equation. The 
theory of the fourth type of diffraction must be 





'S. A. Shelkunoff, Electromagnetic Waves (D. Van 
Nostrand Company, Inc., New York, 1943), p. 365. 
*H. A. Bethe, Phys. Rev. 66, 163 (1944). 
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general and may, when developed, add to the 
understanding of the other two. 

Most theoretical determinations of diffraction 
patterns since the time of Fresnel have been 
based on the simple but necessarily blind as- 
sumption that the field over the surface of the 
aperture is identical with that of the unper- 
turbed incident wave. 

When a sufficiency of precise experimental 
data has been taken of fields in and near the 
aperture, a sound theory of diffraction may be 
built on the basis of experimentally known 
boundary conditions and Maxwell’s equations. 
Bethe has suggested? that it may be possible to 
determine the boundary conditions by rigorous 
theory such as he applied to small holes. How- 
ever, the generalization of his method of deter- 
mination of the boundary conditions appears so 
difficult that it would not be too unsportsman- 
like to make experimental observations first. 

At present the only self-consistent theory of 
diffraction of electromagnetic waves’ is that for 
the case in which the boundary surface is a 
perfect reflector. The condition of perfect re- 
flectivity is closely approximated for microwaves 
incident upon a metal surface. 


APPARATUS 


The diffraction pattern of electromagnetic 
waves was studied at distances from zero to five 


*B. Baker and E. Copson, Huygens Principle (Oxford 
University Press, 1939), p. 152. 
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wave-lengths behind a circular aperture. Figure 
1 is a sketch of the horizontal cross section of the 
arrangement. Microwaves of 12.8 cm were em- 
ployed. The oscillator was two coaxial, coupled, 
resonant cavities.with a General Electric disk- 
seal 2C39 triode as an integral part of the 
cavities. The parabolic reflector, four feet in 
diameter, was fed by a small paraboloid, with a 
focal length of a quarter wave and diameter of 
one and one-half wave-lengths, having a dipole 
at its principal focus. The three-foot iris dia- 
phragm was mounted in a 6 ft. by 6 ft. sheet iron 
“barn door’ which could be rolled to one side 
when the undisturbed intensity was to be meas- 
ured. To supplement the sheet iron screen, a 
wire-mesh screen was mounted around it making 
the total screen 12 ft. by 12 ft. The top of the 
optical bench was 18 in. below the bottom of 
the iris. A vertical wooden dowel 3 ft. long 
supported the receiving dipole antenna so that 
it could be moved in a horizontal plane through 
the axis of the diaphragm. 

As shown in Fig. 2, the dipole was at the end 
of a 2-ft. coaxial line of silver tubing in the other 
end of which was mounted a silicon crystal 
detector. A quarter-wave choke, coaxial with the 
silver tubing, was mounted behind the dipole. 
A shielded line extended from the detector to a 
50-microammeter outside the field. The com- 
bined crystal and meter were calibrated with a 
wave-guide attenuator. The rectified current 
through the crystal was proportional to the 
square of the potential across it and, therefore, 
to the received power and to the intensity of the 
radiation in the region of the dipole. 

Both the dipole of the transmitter and of the 








j 


receiver could be rotated to the vertical or 
horizontal plane. When they were vertical, the 
diffraction pattern in the H or “magnetic plane” 
could be measured. When they were horizontal, 
the pattern in the E or “electric plane’’ could be 
measured. 

To prevent reflections, the iris was set before 
the center of an open window 11 ft. wide and 
8 ft. high. The intensity at the edges of the 
window and along the floor, ceiling and walls 


was negligible. Standing waves were observed, 


of course, between the transmitter and screen. 
Care had to be taken that the opening and 
closing of the iris did not react on the oscillator 
to change its output. It was found that when 
the antenna was coupled sufficiently loosely to 
the oscillator, the plate current of the tube 
remained constant as the iris was opened and 


closed. 
RESULTS 


In the following set of graphs, all distances are 
expressed in wave-lengths. The intensities are 
relative to the intensity of the undisturbed 
beam when the iris and screen were removed 
from the region. 

Figures 3a, b, c, d, and e show the intensities 
in the plane of the iris along the diameter in the 
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H plane when the diaphragm was one, two, three, 
four, and five wave-lengths in diameter. Figures 
3a’, b’, c’, d’, and e’ were of intensities along the 
diameter in the Z plane for the same sizes of 
diameter. To one accustomed to employing 
Fresnel’s zone theory or Kirchhoff’s theorem in 
optics, with the assumption of constant intensity 
over the aperture, it was startling to find that 
the sharpest diffraction pattern in the whole field 
was in the plane of the aperture. Measurements 
were made of the intensity beyond the edge of 
the iris. In these cases, the antenna was kept at 
a distance of about two millimeters from the 
surface of the leaves of the iris. In the E plane 
the beam was broader than in the H plane. 
Moreover, it was noted that the intensity in the 
E plane was not symmetrical with respect to the 
center. This was possibly caused by the fact that 
contact resistance between leaves of the iris 
diaphragm was variable so that the surface cur- 
rents were not the same on the two sides. 


_ Fic. 3. Intensities along the diameters of the aperture 
in the H and E planes for diameter of from one to five 
wave-lengths. 
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Figures 4a, b, c, d, e, and 4a’, b’, c’, d’, and e’ 
are for a diaphragm one wave-length in diameter, 
and give the intensities in the H and E planes, 
respectively, along lines parallel to the plane of 
the iris, at various distances f from the iris as 
shown. 

Figure 5 gives a similar set of graphs of the 
diffraction pattern of an aperture two wave- 
lengths in diameter. Curves a through k, and a’ 
through k’, give the intensities in the H and E 
planes, respectively, taken along lines parallel 
to a diameter and at distances f from the iris as 
indicated on each graph. 

Figure 6 similarly shows the diffraction pattern 
of an aperture three wave-lengths in diameter. 
Following down the column of graphs the peaks 
of intensity could be seen receding and rising. 
To make it simpler to visualize the intensities 
in the H plane, the three-dimensional graph, of 
Fig. 7 was plotted. The broken lines are lines of 
constant intensity over the surface. The dis- 





Fic. 4. Intensities in the H and £ planes taken along 
lines’ parallel to"the ~y of the aperture when the aperture 
in diameter. 
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Fic. 5. Intensities in the H 
and £E planes taken along lines 
parallel to the plane of the 
aperture when the aperture was 


two wave-lengths in diameter. 








DIFFRACTION PATTERN OF A CIRCULAR APERTURE 


Fic. 6. Intensities in the H 
and E planes taken along lines 
parallel to the plane of the 
a re when the aperture was 
three wave-lengths in diameter. 
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Fic. 7. Three-dimensional plot of intensities in the H plane 
when the aperture was three wave-lengths in diameter. 
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Fic. 8. Contour of intensities in the H plane when the 
aperture was three wave-lengths in diameter. 
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aperture. This plot was shown for the width of 
the aperture. 

For the case of a diaphragm of four wave. 
lengths diameter, the only data taken were of 
intensities in the plane of the iris and are shown 
in Fig. 3d. It would be seen that as the iris Was 
opened from one up to four wave-lengths in 
diameter, an array of peaks in intensity would 
move out from the plane of an iris being like an 
array of ten pins when the diameter of iris 
reached four wave-lengths. 

Figure 8 is a contour map of the intensity in 
the H plane when the aperture was three wave- 
lengths in diameter. Figures 7 and 8 and the 
first column of Fig. 6 gave a complete picture 
of the diffraction pattern in the H plane when 
the diameter of iris was three wave-lengths, 

A simple case to treat theoretically is that of 
the intensity along the axis of the aperture. In 
this case, the intensity is independent of polar- 
ization. Moreover the experimental results can 
be checked more readily with existing theory. 
Twelve graphs of the intensity at fixed points 
on the axis at distances f from the aperture were 
plotted against the diameter of the aperture of 
the iris. These are shown in Fig. 9, with the 
distances f indicated in the upper right-hand 
corner of each graph. The positions of the 
maxima and minima checked approximately 
with the predictions on the basis of elementary 
Fresnel’s zone theory. For each distance f the 
minimum intensities corresponded to apertures 
of an even number of Fresnel’s zones. 

These results, as have been previously noted, 
are surprising when it is noted that in Fresnel’s 
zone theory it is assumed that the intensity 
over the aperture is constant. From Fig. 3, it is 
seen that the sharpest diffraction pattern of the 
circular aperture is in the plane of the aperture. 
Fresnel’s zone theory thus becomes merely a 
rule of thumb that gives results in close approxi- 
mation to experiment. 

When the zone theory is applied to the ridicu- 
lous case of f=0, the zone boundaries are circles 
of radii equal to an integral number of half- 
wave-lengths. The diameters of the aperture for 
maximum and minimum intensities at the center 
of the aperture were correctly predicted. Fresnel 
did not intend that his theory, which assumed 
constant intensity in the aperture, would be 
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Fic. 9. 


Intensities at fixed 


points on the axis of the aperture 
at distances f from the aperture 
plotted against diameter of the 


aperture. 
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used to calculate the diffraction pattern in that 
same aperture. 

The irregularities in the graphs of Fig. 9 were 
studied to see if they could be results of experi- 
mental error. The same irregularities were re- 
peated with each experimental measurement. It 
was noted that some of the irregularities ap- 
peared systematically along the curves. For 
instance, in Figs. 9h to 9m, it was observed that 
a small dip appeared in each curve corresponding 
to diameters of four wave-lengths. Likewise, 
although less pronounced, there appeared small 
humps corresponding to a diameter of three 
wave-lengths and dips corresponding to diam- 
eters of two wave-lengths. That is, dips occurred 
when there was a minimum at the center of the 
aperture. Possibly the currents over the surface 
of the iris were responsible for some of the 
irregularities. All that can be said definitely is 
that the irregularities could not all be accounted 
for as experimental error and that there was 
order in the irregularities. 


THEORY 


It is not necessary to mention here all the 
theoretical papers that have called attention to 
Kirchhoff’s false assumption that the intensity 
over the aperture is the intensity of the undis- 
turbed beam. Figure 3 is experimental confirma- 
tion that the intensity over the plane of the 
aperture is not constant when a plane wave is 
incident normally upon the aperture. 

The sole virtue of Kirchhoff’s theory of diffrac- 
tion lies in its correct predictions and not in its 
false assumptions. It has been used to predict 
correctly the intensity of light at points at dis- 
tances behind the aperture large compared to a 
wave-length. Now that measurements have been 
made near the aperture, the question arises, 
how badly is Kirchhoff’s theory in error in 
predicting the intensity near the aperture? 

Kirchhoff’s theory is applied here to the special 
but important case of calculating the intensity 
along the axis, a case for which the intensity will 
be independent of the polarization in the plane 
of the aperture. It seems desirable to present 
the results for two reasons: 

(1) The results are sufficiently close to experi- 
mental observation to serve as a stop gap until a 
rigid theory is evolved. 
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(2) The solution of Kirchhoff’s theory for the 
general case is found to be simpler than for the 
special case in which approximation is made 
when the distance from the aperture js large 
compared to the wave-length. 

The derivation of Kirchhoff’s theorem is found 
in books on optics‘ and electromagnetic theory,’ 
For the special case of a plane wave incident 
normally upon an aperture, the problem is that 
of integrating the effects of Huygens’ wavelets, 
from the unperturbed wave in the aperture, at 
the point at which the amplitude is to be found, 

uo 2x e~** e-* 
ue=—f {| —a +cos8) —i—cos6 }de, (1) 
4a dN p p? 


where up is the amplitude at the point P resulting 
from the added effects from all points in the 
aperture, % is the amplitude of the undisturbed 
incident wave at the aperture, de is an element 
of area in the aperture, p is the distance from P 
to the element do, ¢ is the phase lag of the 
component from de behind the component from 
an element at the base of a perpendicular from 
the point P to the plane of the aperture, and # 
is the angle between the incident beam and a 
line from do to the point P. 

In books on physical optics, the second term 
is always omitted since p>d in all of the cases 
treated in light. However, in the study of ultra- 
high frequency radio waves, useful applications 
occur when p is of the same order as the wave- 
length. Note the direction factor (1+cos9), 
which was discovered by Kirchhoff. A Huygens’ 
secondary wavelet is not uniform over the whole 
sphere but has a maximum amplitude in the 
direction of propagation of the wave that gen- 
erated it, half maximum at right angles and 
zero in the direction opposite to the primary 
wave. A Huygens’ wavelet for the second term 
of the integral has zero amplitude at right angles 
to the direction of the primary wave and equal 
amplitude and opposite phase in the directions 
with and opposite to the direction of the primary 
wave. 

Let the study be restricted to the determina- 


4M. Born, Optik (Verlagsbuchhandlung Julius Springer, 
Berlin, 1933), pp. 147-151. 

SJ. A. Stratton, Electromagnetic Theory (McGraw-Hill 
Book Company, Inc., New York, 1941), pp. 424-428, 
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Fic. 10. (a, b, c, d) Spirals, 
calculated from Kirchhoff's the- 
ory, for the determination of 
amplitudes at fixed points on 
the axis as the diameter is 
increased. (a’, b’, c’, d’) Intensi- 
ties, calculated from Kirchhoff's 
theory, at fixed points on the 
axis plotted against diameter 
= compared with experimental 

ta. 
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tion of the disturbances along the axis of the 
aperture, and let f be the distance of P from the 
aperture. Let the increments of area be con- 
centric rings. For any ring, p, ¢, 0, and the radius 
r of the ring will be constants for the given ring. 
R, is the radius of the aperture. It will be most 
convenient to express all the variables p, @, r, 
and do in terms of ¢. Applying the transforma- 
tion 


orA+2xf 2xf 
p=————_,_ cos0=————-, da = Apdo 
2r ov\+2nrf 
to Eq. (1) after combining the cosine terms, 


uo 8 1 4 2nf 
2 fellate) 
2 +5 A godA+2rf/ orX+20f 


1 
+ |e, (2) 
r 





6 is the phase of the component from the 
outer ring. Express f in wave-lengths. Let f=md 
where is any positive real number. Thus 


B 1 2nn 
e{(1- ) +1 a9. (3) 
0 o+2an/ o+20n 


Note that this solution for the general equation 1s 
simpler than the solution when the second term is 
omitted. 


B+41n 
up= ud (sind-+4 cos) (——~) -i|. (4) 
28+42n 


Uo 
up=— 





We wish to express the intensity along the axis 
in terms of the diameter of the aperture. Let d 
be the diameter. Then 


d/=[(8/r)?+4nB/x }. (5) 


Figures 10a, b, c, and d are plots of real versus 
imaginary components of the amplitude for 
values of 8 from 0 to 7x for distances f along the 
axis from the center of the aperture of 0.5, 1.0, 
2.0, and 5.0 wave-lengths. The ratio of the ampli- 
tude of the diffracted wave to the amplitude of 
the undisturbed wave is measured by the distance 
from the origin to the point on the spiral corre- 
sponding to a particular 8. The broken circle is 
the asymptote approached by the spiral. 
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The square of the ratio of the amplitudes is 
equal to the ratio of the intensity of the dif. 
fracted wave to the intensity of the undisturbed 
wave I/Ip. From Eq. (5), d/X the diameter 
expressed in wave-lengths was computed for 
values of 8 from 0.0 to 9.07. Figures 10a’, b’, ¢’ 
and d’ are plots of J/Io against diameter of 
aperture for four distances f along the axis from 
the center of the aperture of 0.5, 1.0, 2.0, and 
5.0 wave-lengths. The envelope of the curye 
approaches asymptotes indicated by broken lines 
at I/Io equal to 0.25 and 2.25. 

The points indicated by circles are of experi. 
mental data. As a _ theoretical development, 
Kirchhoff’s theory does not take into account 
currents over the surface surrounding the aper- 
ture. It assumes constant intensity, constant 
phase, and constant direction of polarization over 
the aperture, none of which is experimentally 
true. Since Kirchhoff’s theory contains so many 
false assumptions, the deviations in Fig. 10 of 
experimental data from the results of Kirchhoft’s 
theory are to be expected. These results of 
Kirchhoff’s simple theory are still sufficiently 
precise to be useful. 

Any correct theoretical prediction of the dif- 
fraction pattern of the aperture would be made 
on the basis of the known distribution of fields 
in the plane of the aperture and currents over 
the surface surrounding the aperture. The distri- 
butions along the line of intersection of the H 
plane and the plane of the aperture shown in 
Figs. 3a, b, c, d, e appeared simple to analyze. 
However, because the resistance between the 
leaves of the iris was variable, the distribution 
in the E plane was not symmetrical and, there- 
fore, was difficult to analyze. Before the distribu- 
tion in the E plane can be analyzed, data will 
have to be taken for a set of apertures in solid 
metal sheets. 

It is a pleasure to acknowledge the support 
throughout this study of W. C. White and E. D. 
McArthur of the General Electric Research 
Laboratory, to thank Professor J. R. Collins of 
Cornell University for reading the paper and 
Dr. Emory Cook of Johns Hopkins University 
for a written status report on the theory of 
diffraction. 
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It is shown generally that the Riesz method of analytic continuation and the d-limiting 
process give, in general, the same results for both the classical electromagnetic potential of a 
point source and all its derivatives with respect to the coordinates of the field point. 





I. INTRODUCTION 


[ is well-known that the d-limiting process 
developed by Wentzel, Dirac, and Pauli! 
leads to divergence-free results for classical wave 
fields of a point source. In the case of the 
electromagnetic field, the \-limiting process gives, 
for a field point not on the world line of the 
particle, the Lienard-Wiechert retarded field just 
as in the ordinary theory. On the other hand, 
for a point on the world line, where the ordinary 
mathematical treatment gives divergent results, 
the \-limiting process gives one-half of the radia- 
tion field, i.e., half the difference of the retarded 
and advanced fields, which is finite on the world 
line. Dirac? studied in detail the field strengths 
of the radiation field of a point charge. General- 
ized investigations for the fields of point sources 
possessing higher multipole moments were later 
made by Bhabha and Harish-Chandra for both 
the electromagnetic and meson fields.’ 

An alternative mathematical treatment of the 
problem was introduced by M. Riesz.‘ The 
application of Riesz’ method has recently been 
developed by Gustafson® for the quantum theory 
and by Fremberg® for the classical theory. 
Fremberg has shown that, for the wave fields of 
a point charge, the Riesz method gives the same 
results as the A-limiting process for the following 
field quantities: (i) the retarded potentials and 
field strengths for a point not on the world line; 


1G. Wentzel, Zeits. f. Physik 86, 479 and 635 (1933); 
87, 726 (1934); P. A. M. Dirac, Ann. de I’Inst. Poincaré 
9, 13 (1939); W. Pauli, Phys. Rev. 64, 332 (1943). 

?P. A. M. Dirac, Proc. Roy. Soc. 167, 148 (1938). 

*H. J. Bhabha and Harish-Chandra, Proc. Roy. Soc. 
A185, 250 (1946); Harish-Chandra, Proc. Roy. Soc. A185, 
269 (1946). 

*M. Riesz, Conference de la Réunion internat. des 
math. tenue a Paris en Juillet 1937 (Paris, 1939). 

*T. Gustafson, Kgl. Fys. Sallsk. i Lund Forhandl. 15, 
No. 28 (1945); 16, No. 2 (1946). 

*N. E. Fremberg, Meddelanden fran Lunds Univ. Mat. 
Seminarium 7 (1946); Proc. Roy. Soc. 188A, 18 (1946). 


(ii) half the radiation field strengths for a point 
on the world line. 

Unlike the retarded field quantities given by 
the formula of Lienard and Wiechert, the values 
on the world line of the derivatives of the radia- 
tion field potentials, with respect to the coordi- 
nates of the field point, cannot be derived from 
the corresponding values of the potentials by 
simple differentiation. The agreement of the two 
methods found by Fremberg therefore raises the 
question of the equivalence of the two methods 
in general. It is the main purpose of the present 
investigation to give a general proof of the 
equivalence of the two methods for the electro- 
magnetic field generated by a point source. The 
case of the meson field can be studied in a 
similar manner. 

It follows from the general equivalence of the 
two methods that the Riesz method provides us 
with an additional method of calculating the 
radiation fields. This is of practical importance 
for the sake of computation, as the Riesz method 
is usually simpler and more straightforward for 
evaluating the higher derivatives of the radiation 
fields, which appear in the equation of motion of 
a particle having a multipole moment. 

Before giving the proof, we shall first derive 
in Section II a new expression for the Riesz field, 
which will help bring out the connection between 
the two methods. Explicit results for the po- 
tential and its first two derivatives on the world 
line will be given in Section III, which will be 
seen to agree with those of the A-process. A 
general proof of the equivalence of the two 
methods will be given in Section IV. 


Il. THE RIESZ FIELD 


We begin by considering the electromagnetic 
field generated by a source distribution J which 
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may be a scalar, vector, or tensor of higher 
rank, and the field potential A determined by 


the wave equation 
OA=4rJ. (1) 


The Riesz potential is defined as 


A@(P) = A(Q)R*-dQ, = (2) 





4\a@ Ds? 


where H,(a) denotes the factor 


Hi(a) -r2-r(=)r(“-1), (3) 
2 2 


Ds? is the domain of integration bounded by the 
retrograde cone with its vertex at the field point 
P and a space-like surface S, and R is the 
hyperbolic distance between the points P and Q. 
In the case of a point source moving along a 
world line Z(r), the source density is of the form 


J(X)= f S(1)8(xo—20)8(x1—21) 
X 5(x2—22)5(x3—23)dr, (4) 


where 7 is the proper time of the particle. In 
order to simplify our considerations we assume 
that there is a single source which begins to 
radiate at an instant r4 and we take S to be a 
surface which cuts the world line at the point 
Z(ra). It will also be assumed that there is no 
incoming field present. Under such conditions, 
the Riesz potential becomes 


A@(x)=——™ f " §(2)Re-"dr, (5) 
a a ee ee 7) R*-*dr, 
H,(a+2) TA 


where 1, is the retarded proper time and R is the 
length of the vector X—Z. We denote the time 
and space components of R by 7o and r and 
denote the length of r by 7, so that 


R=(re—?*)!. (6) 


Equation (5) holds for the region x9—2o0(ra) >0, 
(K—Z(ra))?>0, i.e., the region where the 
Lienard-Wiechert potential does not vanish 
identically. The derivatives of A®(X) with re- 
spect to the coordinates of the field point X are 
of the form 
004 
F(odur-++(X) =— —.- - A(X). (7) 


OX, OX, 
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The field quantities are to be calculated firgt for 
a value of a large enough for the integrals to be 
convergent, and the final results are obtained by 
analytic continuation to a=0. 

It will be convenient to base our later investi- 
gation on a somewhat different expression for the 
Riesz potential, namely, 


Ra 
A(X) =a f A™(X,R)R-UR. 


Here R plays the role of a variable of integra. 
tion, Ra is the length of the vector X—2Z(z,), 
and A™*(X, R) denotes the value of the potential 


A™(X, R)=S/(V, R), (9) 


at the retarded time when the vector X—Z(r) 
has the magnitude R, and V=(v, v) denotes 
dZ/dr. The retarded potential given by Eq. (9) 
differs from the ordinary Lienard-Wiechert po- 
tential in that the variable R is in general not 
zero. 

To derive Eq. (8) from Eq. (6), we note that 
(6) may be written in the form 


A(X e ” Sind 
)"FatD f er 


xf" d[ro— (r?+p?)*] 


re 





where pa is the length of X—Z(r,). This may be 
verified with the help of the general formula 
which states that 


f : f(»)5Lg(p) jdp 


’a sf (p) 
of ela 
i} PIR Y TE akties 


ion f (p1) 
(dg/dp) po: 


according as g(p) has a zero p; between 0 and pa 
or not, f and g being two arbitrary functions 
with g subject to the condition dg/dp>0. We 
carry out the integration with respect to r im 
(10) by applying (11) with p replaced by 7 and 


or 0 (11) 
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using the relation 


vo(r?+ p*)#—(v, r) 
(r?+ p?)* 

_W,R)e-, 

(P+ pt 





d 
—[(r?+p?)*-ro]= 
dr 


(12) 


The result is 


Pa Ss 
A(X) = af 5 = p*—'dp, 
0 ; —_ 


which is just Eq. (8) with A(X, R) given by 
Eq. (9). The factor 


4nr 4ra 


Hi(a+2) 4 a P 
vedn(>+1) 


has been replaced for brevity by a, as the ratio 
of these two factors becomes unity when a=0 
and so its omission does not make any difference 
in the final result. 

It follows from Eq. (8) that the derivatives of 
A(X) in Eq. (7) are given by 





(13) 


F(X) =af : Fret(X, R)R*"dR, = (14) 


where Fret(X, R) denotes the field quantity that 
can be derived from the formula of Lienard and 
Wiechert by differentiation with respect to the 
x, with R kept fixed. We have neglected in Eq. 
(15) terms arising from the variation of R4 with 
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the x,. These terms are independent of the 
variable R and therefore vanish when a is put 
equal to zero. 

-The analytic continuation to a=0 is governed 
by the following general theorem: For any func- 
tion f(x) which can be expanded into a Laurent 
series of the form 


fx)=> a,x" (m=0,1,2,---), (15) 


we have 


x 
af f(x)x*—'dx—va (16) 


when analytic continuation is performed to a=0, 
irrespective of the magnitude of the upper limit 
X. The coefficient a may be regarded as the 
finite part (“partie finie”) of the function f(x). 
It follows from this general theorem, and Eqs. 
(8) and (14), that the Riesz method consists in 
replacing the Lienard-Wiechert potential and its 
derivatives by their finite parts when R is taken as 
the independent variable. 

It was pointed out by M. Riesz* that the use of 
the square of R as the independent variable 
brings the retarded potential to a very simple 
form. Writing 


R=R (17) 
and 
S(r) = —4$dT/dr, (18) 
we have from Eq. (9) 
A™*=dT/d@. (19) 


Ill. EXPLICIT EVALUATION OF THE FIELD QUANTITIES ON THE WORLD LINE 


Equation (9) may be written for brevity. 


Att=S/x, (20) 


where 


xk=(V,R). (21) 


The derivatives of At™*t can best be calculated by a method used by Dirac in his classical theory of 
radiating electrons. The first two derivatives of A™* are 


Sr+ 
cae =-— —), (22) 


OX, xk dr 





a 


*I am indebted to Professor Riesz for this remark. 





3? 1 d Sr*r’ 
et) SEE m 
Ox,0X, drixdr k 
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For the explicit evaluation of the radiation field quantities at a point X on the world line, it is 
. convenient to take as the independent variable ¢=1)—1, ro being the time when Z(r) =X. Since 


RdR= kde, 
we have - (24) 


F(X) =af KF R*—*do, (25) 
0 


with F given “by Eqs. (20), (22), and] (23) for. the potential and its first two derivatives, o, being 
the value of ¢ corresponding- to ra. Using the expansions 


-0-(D) AZ): “eer a 
R=X-Z= Voe——(— “de +-(S) #3 ~(=) ae tm 


mii {1-2 (= ) += (Ss ate], (28) 
dr be dr? , 
pneu {You see ...| - 
24 dr dr’ 


where the suffix 0 denotes the value at time 79, we obtain from Eq. (25) and the general relation 
of Eq. (16) 
A(X) = —(dS/dr)o, (30) 


0A(X) > qd as 
2D (2) oof) e no 
ry ite dr? 


and a rather lengthy expression for 0°A/dx,0x,, which will not be written down here. These results 
are just half the corresponding quantities of the radiation field given by the A-limiting process, 


namely, 
A(X) = 3 (S/x) ret +(S/k)aav}, (32) 


Fe--(X) 10 20 ~) 4 Ss | - 
2 ax, ax, (: we (~)..} ae 


The calculation here is simpler than the previous calculation of the radiation field because the 
Riesz field has finite values on the world line from the beginning, so that we do not have to consider 
its values in the vicinity of the world line. 


IV. GENERAL PROOF OF THE EQUIVALENCE OF THE TWO METHODS 


We give now a general proof that the two methods give the same results for the classical electro- 
magnetic field. It is obvious that the two methods are equivalent for a point not on the world line, 
since at such a point the retarded field quantities given by the formula of Lienard and Wiechert 
are finite, and therefore the finite parts given by the Riesz method are just the retarded field quantt- 
ties, in agreement with the \-limiting process. 

This is no longer the case for a point on the world line where the Lienard-Wiechert potential 
becomes infinite. To show the equivalence of the two methods in such a case we use the expansion 
method introduced by Dirac and developed in detail by Bhabha and Harish-Chandra. 
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Let us consider a field point X near the world line at a distance ¢ from its “contemporary point” 
Zc, i.e., the point at which «=0. Using the suffix c to denote the value of a quantity at the contempo- 
rary point and writing r—7r-=Ar, we have, as shown by the above authors, 


s- s.+(= =) art3(=) (ants-(—=), (Ar)? 


1//dV\? 
k=—(1- c)ar+heo"(ar)*+-( (—) sels cee, 
6 dr/. 


R?= —&+(1—«')(Ar)?— iany'——(( -) +40" Janes si (35) 


where k.’, Ke’, etc. denote (R, dV/dr)., (R, d°V/dr?);, etc. From Eq. (35) with R set equal to zero, 
Ar=te(1+ 3K’ +3x-/2) + 3x e®+(1/24) (dV/dr).268+---, (36) 


where the upper and lower signs are to be taken for the advanced and retarded times, respectively. 
From Eqs. (34) and (36), 


S S. dS_ Sx’ Sa’? cf” /4V\" @ dS e«d’S 
(2) ELEY D-Ses]e m 
alm «¢ dr P 8e 3 dr/ 24 dr 2dr*l, 
The value of the same quantity at the advanced time is obtained by writing —e for ¢ in the right- 
hand side. It follows therefore, 


irsS Ss dS 1 dS 
TOO Nn scan 
2L\ x J ret K J sav dr/. 3 


With the help of the formulas 














OT. v 

7 ( (39) 
OX, 1 — x’ e 
or? 

r-(— (40) 

ax, 1—x«’ 
de? 
—=-2r/, (41) 
Ox, 


Bhabha and Harish-Chandra reached the conclusion that the expression in Eq. (38) and all its 
derivatives are of the form 


co 
ie a,e", 
0 


where the a, are functions of r, and R, which are finite on the world line. It follows therefore that 
the radiation potential and all its derivatives are finite when e=0, i.e., when the point X lies on the 
world line. 

For the Riesz field R is not necessarily zero. For the purpose of evaluating the finite parts, however, 
it is sufficient to consider values of R in the neighborhood of zero. We can then take R to be a small 
number and solve again Eq. (35) for Ar, the result being 


Ar=+(1+gn' + $x) (PO +R*) + ne” (+R?) + (1/24) dV /dr) (+R +: -- (42) 
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for the advanced and retarded times respectively. Equation (42) differs from Eq. (36) only in having 
— instead of ¢. Similarly, 


S. ds Sx’ 
KJ ree (€8+R?*)! dr 2(e+R?)!J. 


X” ss dS one : 
(= ) = =| +E an(P+R)" +z ba(2+ Rem, 


+s (saa 3 


dr? 
(43) 
where the b, denote functions similar to the a,. Using the formulas (39), (40), (41) and 
He+R) 2 a, al 
OX, OX, 
we see that both the expression in Eq. (43) and all its derivatives are of the form 
¥ an(P+ RITE belt R)+4, (45) 


when m is either zero or a positive integer. 


The expansion in Eq. (45) contains the results of both methods for the field quantities on the world line. 
The A-limiting process consists in putting R=0, taking the mean value 


and finally putting e=0. On the other hand, the Riesz method consists in putting e=0 and replacing 
the expansion 


¥ a,R*+>D b,R**! (47) 
v0 -—m 
by its finite part. We see therefore quite generally that the two methods are equivalent. 


I am indebted to Professor W. Pauli for a discussion which suggested the present investigation, 
and for his advice. 
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The problem of the relativistic invariance of a first-order wave equation with matrix coeffi- 
cients 8 is examined. It is found that it is intimately connected with the structure of the en- 
veloping algebra of the 8-matrices. In particular the center of this algebra can contain only 
elements of a very restricted type. This fact provides a powerful means for the investigation of 
the abstract §-algebra, as is illustrated in the general s-dimensional case of the Dirac and the. 
Duffin-Kemmer matrices. Relativistic invariance imposes severe restrictions on the spurs of the 
B-matrices and their multiple products. Conversely, these restrictions insure relativistic 
invariance. The theory is applied to the case of particles of higher spin. 





1, INTRODUCTION 


HE wave equation for a particle of arbitrary 
spin can be written in the form 


iB'aW~+xv=0, (0,=0/dx*) (1) 


where x is a constant determining the mass of the 
particle, and 6*(k=0, 1, 2, 3) are square matrices 
operating on y which itself is a one-column 
matrix. The object of the present paper is to in- 
vestigate the conditions which the #’s must 
satisfy in order that Eq. (1) may be relativistically 
invariant. The conditions so obtained are also 
found to be largely sufficient. They insure the 
invariance of Eq. (1) under proper Lorentz 
transformations. 2 

The method of this paper seems to provide a 
powerful means for the investigation of the 
algebra generated by the §’s, and therefore for 
the determination of their irreducible repre- 
sentations. In order to find all the irreducible 
representations of a semi-simple algebra one has 
to obtain a linearly independent basis for its 
center. It will be shown, however, that only 
quantities of a very restricted type can belong to 
the center, and therefore the task of obtaining a 
basis becomes much easier. Moreover, once such 
a basis has been obtained, one can immediately 
decide whether or not a particular representation 
is equivalent to its complex conjugate, transposed 
or Hermitian conjugate. The question of the 
invariance of the representation under reflection 
can also be answered at once. It will be shown 
also that considerable information concerning the 


' The center of an algebra consists of the set of elements 
which commute with all elements of the algebra. 
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spurs of the §’s and their multiple products can 
be deduced from the relativistic invariance of 
Eq. (1). Conversely, if the spurs of the 6’s and 
their multiple products satisfy certain conditions, 
the invariance of Eq. (1) for the full Lorentz 
group is ensured. 

As an illustration of these general results, the 
algebras of the Dirac and the Duffin-Kemmer 
matrices in the general case of s dimensions are 
discussed briefly. Although most of our results 
have already been obtained by Schrédinger* and 
Kemmer,’ they serve to illustrate the power of 
the new method. 

In the scheme of equations proposed by Dirac‘ 
one must impose certain other subsidiary condi- 
tions on y in addition to Eq. (1). However, as is 
well known, these subsidiary conditions in 
general become inconsistent with Eq. (1) when 
the electromagnetic interaction is introduced in 
the usual manner. A way out of this difficulty has 
been indicated by Fierz and Pauli,’ but their 
method, though logically consistent, is rather 
cumbersome and consequently lacking in ele- 
gance. It therefore seems worth while to investi- 
gate alternative methods of approach to the 
problem. Bhabha® has recently proposed a form- 
alism in which the subsidiary conditions are 
completely eliminated, so that the interaction 
can be introduced in the usual manner without 
difficulty. But the particle in general no longer 


( ons) Schrédinger, Proc. Roy. Irish Acad. [A] 48, 135 
1943). 
’N. Kemmer, Proc. Camb. Phil. Soc. 39, 189 (1943). 
‘P. A. M. Dirac, Proc. Roy. Soc. [A] 155, 447 (1936). 
( +) Fierz and W. Pauli, Proc. Roy. Soc. [A] 173, 211 
1939). 
*H. J. Bhabha, Rev. Mod. Phys. 17, 200 (1945). 
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obeys the second-order wave equation 
n.0*Y + xy =0 (2) 


in the force-free case. The physical result of this 
circumstance is that the particle described by 
Bhabha’s equation would possess several different 
mass values. However, on grounds of simplicity, 
it seems desirable to try to discover a scheme in 
which the difficulty concerning the subsidiary 
conditions may be avoided without abandoning 
the second-order equation. Such an attempt is 
made in this paper. 

As shown by Pauli,’ the connection between 
spin and statistics is really based on the rela- 
tivistic invariance of the wave Eq. (1). It there- 
fore remains valid in the present case, as indeed 
in every other where the wave equation is 
relativistically invariant (e.g., for the equations 


proposed by Bhabha).*® 
2. ALGEBRAIC THEORY 


For the relativistic invariance of Eq. (1) it is 
necessary that for every Lorentz transformation 
(t') there should exist a non-singular matrix T 
satisfying the condition 


By’ =t'B' = TB, T-. 


However, the problem considered in this paper 
may be stated somewhat more generally. Jf 
By (R=1, 2, ---, s) be an irreducible® set of s mat- 
rices, what are the conditions which the B's must 
satisfy in order that, for every transformation of the 
real orthogonal” group, there should exist a non- 
singular matrix T such that 


By’ =t,'B:= TB,T—, (3a) 

yeti" mn = Leis (3b) 
0 k+l, 

BN) k=l; @) 


and are these conditions sufficient to insure the 


7W. Pauli, Phys. Rev. 58, 716 (1940). 

*H. J. Bhabha, Proc. Ind. Acad. Sci. 21, 291 (1945). 
* ®In order to insure that Eq. (1) describes an “elemen- 
tary” poreee. which should be an entity of the simplest 
possib: kind, we assume that the set of matrices f is 
irreducible. 


10For the moment only the real orthogonal group is 
one but 4 = be shown Rncaner 3 that Ger 
corresponding problem concerning other real grou or 
example, the Lgcunte group) associated with indefinite 
ground forms, is intimately connected with the present one, 
and has an almost identical solution. 
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invariance" of the B-representation under the orthog. 
onal group? 

Let (8) denote the enveloping algebra of the 
8-matrices,” i.e., the algebra consisting of all 
quantities of the type 


C=CBET CB BitC™BBiBmt+--+, (5) 


where c**** are numerical coefficients belonging to 
the field § of complex numbers."* Since the 6; are 
irreducible, (8) is simple. Let u, (p=1, --., h) 
be a linearly independent basis for Y. As elements 
of (8), u», must be expressible in the form (5) 
and therefore one may write u,=1u,(8). On re- 
placing each & by #;’ in u,(8) one obtains 
u»(B’) =u,(t8). Now u,(t8) is also an element of 
% and therefore is expressible linearly in terms of 
the basis u,(8) ; that is, 


Up(tB) = dig Tp%u,(8), (6) 


where"! 7,%e(§, and may be regarded as coeffi- 
cients of a transformation T in the u-space. Thus 
every transformation ¢ in the B-space induces a 
transformation T in the u-space. It is easy to 
verify that the correspondence t->T is a homo- 
morphism with respect to addition and multi- 
plication and therefore the set of transformations 
T is a representation of the orthogonal group 
O(s). Now, it is well known that every repre- 
sentation of O(s) is fully reducible and therefore 
the u-space can be decomposed into a direct sum 
of irreducible invariant sub-spaces. Let U denote 
the entire u-space and let the irreducible in- 
variant sub-spaces be denoted by V(1), V(2), 

-+, V(m) such that!® (V(m):§) =m. Now all 
irreducible representations of O(s) are known." 
Any such single valued representation of the 


1 The property of the existence of a non-singular matrix 
satisfying Eq. (3) corresponding to every transformation 
t.' of a particular group G, can be briefly denoted as the 
property of the “invariance of the 8-representation under 
G.” The full (real) orthogonal group in s dimensions will be 
denoted 4 O(s) while O*(s) will denote the corresponding 
proper orthogonal group. 

2 The whole set of matrices 6:(k=1, 2, ---, s) is con- 
veniently denoted by 8. In the same way it is also convenient 
to abbreviate Eq. (3) to 8’ =#8, where ¢ is now to be looked 
upon as a transformation operator acting on the set 8 and 
transforming it to 6’. 

1% The summation convention for contracted indices is 
used throughout in this paper. 

4 This bolism indicates that the coefficients T;,* 
belong to the field §. 

16 If V is a vector space over §, then (V: §) denotes the 
dimension of V over #. i 

16H. Weyl, The Classical Groups (Princeton University 
Press, Princeton, New Jersey, 1939). 
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entire group O(s) is equivalent to the set of 
transformations of a tensor of suitable rank and 
symmetry properties. It is therefore possible to 
find a linearly independent basis “v7, (p=1, ---, 
hm) for V(m) such that 7 transforms “», as 
components of a tensor of rank f,, with suitable 
symmetry properties. Thus, one can define a 
tensor “vk;---kym (Of suitable symmetry) such 
that its 4, independent components are equal to 
the “™v,’s and 


Vie) (48) = tye! - 9, (8). (7) 


Here the superscript (m) has been omitted for 
simplicity, and the symbol [k] is used to denote 
the entire set of tensor indices ki, ---, Ry. tu!" is 
merely a convenient abbreviation for tiy'tthy!*: - « tk’ 
At this point we need the following lemma. 
Lemma 1. If Dix(8) ts a set of quantities be- 
longing to U(8) which transform as a tensor of rank 
f for all orthogonal transformations of the B-space so 
that 
Dy (t8) = tp!" - Di (8), (8a) 
then 


Diy(8) = day (g) +4 win(g) -B" 
+d witite(g)-B"8%+---, (8b) 


where dix)...(g) are tensors constructed out of gx 
alone i.e. 


dx..-(g) -  ® a(a, B, 7, 6, +> * )Saphys° — (9) 


where a(a, B, y, 5---)ei} and the sum is to be taken 
over all permutations a, B, y, 5-+- of the indices 
ki++ky, «++ occurring on the left hand side of 
Eq. (9). The sum in Eq. (8b) is, of course, finite. 
When the total number of indices ki, ---k,;, «++ is 
odd, a relation such as Eq. (9) is impossible and so 
in this case d,x)...(g) =0. 
Proof: Put 
D=§&").Di(8), 
where 
gle = ECA... EDA, 


and ¢‘”* are indeterminates transforming contra- 
variantly”’ to 6,. Clearly D is invariant under the 
transformations of O(s). Since O(s) is a topo- 
logically compact group, it is possible to perform 
the process of averaging over the entire group 
manifold.'* Let (D) be the average of D over O(s). 

Although the difference between co-variant and 


contra-variant vectors disappears for orthogonal trans- 
formations, it is convenient to recognize it formally. 


RELATIVISTIC WAVE EQUATIONS 


From Eq. (5), D must be of the form 


D= El"). {de 1nB%+d en wB"B4+---}, (10a) 


where dix)... are suitable numerical coefficients. 
Therefore 


(D) = &#). {(d) wing 
+(d) hineB"B4+---}, (10b) 


where (d);x)... are the averaged coefficients. Con- 
sider the forms 


* D(n) = §)- fd wing 4 


+d [k [tylen? in a4 eee , (11a) 
and 


(D)() = &)- {(d) ei? 4 
+(d) Ue ty ton) Wg aew eed, (11b) 


where »‘”* are also indeterminates transforming 
contravariantly to 6,. Notice that Eqs. (11a) and 
(11b) are derived from Eqs. (10a) and (10b) by 
replacing the first 8 in each multiple product of 
the 6’s by n*, the second 8 by »*, and so on. 
Also (D)(n) is the average of D(n), and, therefore, 
clearly (D)(n) is an absolute invariant of the 
orthogonal group O(s). Since x*y,=x*y'g,; is the 
only typical basic invariant of O(s),'*(D)(m) must 
be expressible as a rational integral function of 
the scalar products of the é’s and y’s. This, in 
turn, means that the coefficients (d),4)... are of the 
form of Eq. (9). Now, because of the invariance 
of D, D=(D) and therefore Eq. (8b) holds. 

Lemma 2. If Eq. (8a) is true at least for all 
transformations teO*(s), then 


Dui (8) = day (g, €) +d win(g, €)-B"+-++, (12) 


where d;x)...(g, €) are tensors constructed out of gx: 
and ¢,---k,'* only. 

The proof is similar to that of Lemma 1 with 
the difference that Ot(s) now has two typical 
basic invariants x*y, and €&y---k&yx(DP ig Aa. « « (oho, 
Since the product of two e’s can be expressed in 
terms of the g’s, we need consider only terms 
linear in «. 

From Eq. (7) and Lemma 1 it follows directly 
that 

v4) (8) =r 4) (8, g) (13) 


i.e., as in Eq. (8b) v44;(8, g) is a tensor composed 
of the #’s and g’s only. Now “»,4)(8, g) taken to- 


18 ¢k)...k, is anti-symmetrical in every pair of indices, and 


€12.-..= i 
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gether for all values of m form a complete basis 
for %(8). Hence the product of any two v’s is 
again expressible as a linear sum of the v’s so that 


™ 974) (8) - ™v¢ (8) 
+Yo me (m, m mT yy ° (mo 59 (8) 


where &™ ™? (+) T(x}, !! are certain numerical coef- 
ficients. (Here the summation convention is used 
for the j’s.) On imposing suitable symmetry 
properties on the I’s corresponding to those of 


(14) 


the v’s in respect to permutations of the j’s, they. 


are uniquely determined by Eq. (14). Now from 
the similarity of the matrices 8 and #8 it follows 
that the correspondence $8 is an auto- 
morphism of %(8). Therefore 


(™) 914) (4B) - 0, (¢8) 
De OO ome T pay) = 0, (4B). (15) 


So that from Eqs. (7), (14), and (15) one finds on 
omitting the m’s that 


bay!) -tey Pay’) =e Puan. 


Put 

Prat tay = Sta) “Tas, (16) 
where 

Bla) (o] = LarbiZarbe: + * Lardy. 
Now put 


P= EQ) Pegi enter, 


where {‘”* is a third set of indeterminates trans- 
forming contravariantly to B;. Then I is an 
absolute invariant of O(s) and, therefore, 


LP peicaitay = T peje tai (g) (17) 


where the right side denotes a tensor constructed 
out of the g’s only. Equation (14) can therefore 
be written as 


™y4)(8, g) v4 (B, g) 


= Dome OO eT paercatsa(g) * ™ vt(B, g). (18) 
Equation (18) has the form 
Ckike-+-(B, g) =0. (19) 


Since “;)(8, g) form a complete basis for 2(8), 
Eq. (18) describes the structure of %(8) com- 
pletely. Therefore the relations determining com- 
pletely the structure of an algebra can always be put 
in the tensor form of Eq. (19), provided the corre- 
spondence B—>t8 is an automorphism of U(8) for 
every orthogonal transformation t. 

Conversely, let one or more equations of the 
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type (19) be given to define B,. From Eqs. (19) 
and (3b) it is clear that 


Chika ---(1B, g) = thy "Mtko!*~+ -cirta---(B, g)=0, (29) 


so that /8 satisfy the same relations as Eq. (19). 
Hence 6-8 is an endomorphism of (8). Simi. 
larly, B—>t—' is also an endomorphism. The prod- 
uct of these two endomorphisms is the identical 
automorphism $->8. Hence each of the two 
endomorphisms must be an automorphism. Next, 
assume that the defining relations of Eq. (19) are 
such that the abstract algebra 2(@) is finite and 
simple. Now every automorphism of a simple 
algebra is an inner automorphism. Hence there 
exists a unit (i.e., a regular element) 7(8)e%(g) 
such that 

8 = T(8) -B-T—*(8). (21) 


This shows that in any representation of (8) the 
matrices 8 and é8 are similar. Thus in order that an 
trreducible B-representation be invariant under the 
orthogonal group O, it is necessary and sufficient 
that the algebra A(B) be finite and simple and its 
structure be completely described by relations of the 
type (19). 

Notice that if e(8) is an element belonging to 
the center of 2(8), then from Eq. (21) 


e(t8) = T(8) + e(B) -T—*(B) =e(8). 


Hence, e(8) is invariant under transformations of 
the orthogonal group, i.e., it transforms as a 
scalar. Hence, from Lemma 1, 


e(8) =e(g, B), 


that is, e(8) is a scalar constructed by contracting 
all the indices of 6* in the multiple products 
among themselves. As is well known, the center of 
a simple algebra over the field § of complex 
numbers consists of §}. 1 where 1 is the unity 
quantity of %(8). Hence 1=e(8) where e(8) is a 
‘scalar.’ One can therefore allow a multiple of 1 to 
appear in Eq. (19) without spoiling the tensor 
character of the equations. 

Since it is not usually easy to see from a set of 
equations such as (19) whether the algebra (8) 
is simple or not, we drop the restriction about the 
simplicity of 2. Assume, therefore, that a set of 
equations such as (19) is given, defining an 
algebra %(8) of a finite order. We wish to in- 
vestigate the invariance of an irreducible repre- 
sentation of %(8) under O(s). Let B,—'‘B, be an 
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irreducible representation of %{(8) by the matrices 
‘B,. Since p—t"'8 is an automorphism of Y%, 
tp-'B is also an irreducible representation of 2. 
Hence B—'8 and $—+'8 are both irreducible 
representations of %. Choose 


te'=Sn'+ex', (€x1= — xx), (22a) 


where ¢,: is an infinitesimal quantity (¢,' corre- 
sponds to an infinitesimal rotation). Then 


BiB, (22b) 
Bi ‘Bet ex! ‘Bi, (22c) 


are both irreducible representations of YU. Now 
every finite algebra has only a finite number of 
irreducible representations apart from equiva- 
lence. Therefore from continuity, Eqs. (22b) and 
(22c) must be equivalent representations. Since 
all proper orthogonal transformations can be 
generated from the infinitesimal transformations, 
the representations 6—'‘8 and B-—?t's where 
tteOt(s) are equivalent. Hence every irreducible 
representation of (8) is invariant under the proper 
orthogonal group O*. 

Let 9t(8) be the radical of Y. Clearly the radical 
is invariant under any automorphism of Y. 
Hence exactly as before, one can find a tensor 
basis “™nx)(8, g) for Nt. In any irreducible repre- 
sentation of % all elements of Jt are represented 
by zero. Hence all the irreducible representations 
of 9 are the same as those of 6(8), where S(8) is 
defined by Eq. (19) and the further equations 
(™)n1x)(B, g) =0. Clearly the radical of S(8) is zero 
and therefore © is semi-simple. Also @(8) is 
defined by proper tensor equations, so that 
p—#tB is an automorphism of ©. Now every semi- 
simple algebra is expressible as a direct sum of 
simple algebras.'* Therefore, 


© (6) =i(8) @A2(8)B---@A-(8) (23) 


where %..(8) (w=1, ---, r) are all simple. This 
decomposition of © is unique apart from the 
order of %,, in the decomposition. Hence under 
any automorphism @ of ©, the simple algebra 
goes over into some Y,,’ or in symbols 


A.0=W,’. (24) 
In particular let @ be the automorphism p—+#8, 
where ¢ is given by Eq. (21a). Then as before it 


*A. A, Albert, Structure of Algebras (Am. Math. Soc. 
Coll. Pub., New York, 1939). 
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follows from continuity considerations that u’ = 
in Eq. (24), and so the transformation B—* for 


t+eO*(s) isan automorphism of %,,. Corresponding 
to Eq. (23) one has 


B=BP+BH+---+8, BweM,. (25) 


Hence %, is the algebra generated by 6™. The 
automorphism B—t+8 of S corresponds to the 
automorphism p™—+p™ of Y,. Since the latter 
must be an inner automorphism of %, there 
exists a unit 7+, such that 


itp™ as TH) BCH) (T+0s))—1, (26) 
Put 
Tt =THO47+@4...4 7+), 
then 


(T+)-! = (TH) “14 (T+). 4 (THO)—1, 
since if a,e%, and b,e%,, then a,b, =0 if uv, and 
1=e,+---+¢,, (27) 


where 1 is the unity quantity of S and e,eW%, is 
the unity quantity of U,. Thus 


Hp=T+B(T*)—. (28) 


Now if c(8) is an element belonging to the center 
of S, then c(t+s) = T+c(8)(T+)-'=c(B). Hence 
c(8) is invariant under proper orthogonal trans- 
formations. Therefore, from Lemma 2, 


c(B) =c(B, g, €). (29a) 


By the improper orthogonal transformation J, 
which reverses the direction of the x-axis, c(8) 
goes over into 


c(JB) =c(B, g, —e). (29b) 


The ceater of S is obviously invariant under all 
automorphisms of ©. Hence c(J8) also belongs to 
the center of S. Therefore both 


Cec(B) = $c(B, g, €-) +4$c(B, g, —e), 


and 


Cye(B) = $c(B, g, €) —3c(B, g, —€), 


belong to the center of S. c,. is the “‘scalar’’ part 
of c, while c,, is the ‘ pseudo-scalar’’ part. It is 
therefore sufficient to consider only those quanti- 
ties of the center which are either pure scalars or 
pure pseudo-scalars. 

The number of inequivalent irreducible repre- 
sentations of © is r, which is the number of simple 
sub-algebras in the decomposition Eq. (23). Also, 
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(é:, -+ +, é) form a linearly independent basis for 
the center of ©. Hence the number of irreducible 
representations is the same as the dimension of 
the center over §§. One has, therefore, merely to 
obtain all scalar or pseudo-scalar quantities in © 
which commute with 6; and choose from them a 
linearly independent set. Let (c,(8), -- +, ¢-(B)) be 
such a set. In an irreducible representation of © 
these quantities are represented by multiples of 
the unit matrix, i.e., by numbers in §. In general, 
two irreducible representations of © are equiva- 
lent only if c,(8)(u=1, ---, 7) are represented by 
the same numbers in the two representations. In 
particular, suppose the c’s are so numbered 
that ¢:(8), c2(B), ---, ¢,(8) are scalars, while 
Cu+1(B), -**, €r(8) are pseudo-scalars. Then an 
irreducible B-representation is reflection-invari- 
ant only if ¢44:= +--+: =c,=0 for it. 

From the point of view of the physical applica- 
tion it is desirable that the matrices 8; should be 
equivalent to their transposed 6,* and Hermitian 
conjugate B;'. Hence B;? also satisfy Eq. (19) ie., 


Cei...(B*) =0. (30a) 
Taking the transposed of (30a) one obtains 
Cxi...7 (B) = 0, (30b) 


where ¢x:...7(8) is derived from c;,:...(8) by re- 
versing the order of the factors in the multiple 
products. Hence from Eqs. (19) and (30b), 


Crr---(B) +Cx2...7(8) =0, (31a) 
Cx1...(B) —Cxr...7 (8) =0. (31b) 


Since Eq. (19) follows from Eqs. (31a and b), the 
latter define @(8) completely. Hence the defining 
equations of S can always be put in the form 


Cx1...(B, g) =0 (32) 


where Cx:...7=-+c,:.... From the equivalence of 
Bx, Bx and @;,' it follows that 6; is equivalent to 
its complex conjugate §,*. Hence from Eq. (32) 
Cxu...(B*, g)=0. Taking the complex conjugate 


one finds 
Cei...*(B, g) =0 (33) 


where ¢x:...* is derived from c;,:... by applying the 
automorphism i——i of § to the numerical 
coefficients in c;,;.... Hence from Eqs. (32) and (33) 


Cer-.-(B, g)+cxr...*(B, g) =0, (34a) 
iCx1...(B, g) —4Cx1...*(B, g) =0. (34b) 
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Equation (32) follows from Eq. (34), and so (34) 
defines S(8) completely. Therefore the defining 
equations of S can be written in the form 


Cx1...(B, g) =0 (35) 


where Cxu...7 = Cx... ANd Cy...* =Cy.... Thus (35) 
is an equation with real coefficients. If the algebra 
©(8) so defined by Eq. (35) be simple, then apart 
from equivalence it is capable of only one 
irreducible representation, say 6.6, where 6, 
are irreducible matrices. Clearly §,' and 6,? are 
also irreducible and satisfy Eq. (35), so that they 
also constitute a representation of S. All these 
representations must therefore be equivalent. 
Hence if (35) defines a finite simple algebra and 
Bx be an irreducible set of matrices satisfying (35), 
then B, 6°, B*, B', and tB (for any teO(s)) are all 
equivalent. 

In the general case when ©(8) is semi-simple, 
Eq. (35) is not sufficient to ensure the equivalence 
of B, 6", and g'. Of course, as shown above, it is 
still necessary ; i.e., if for every irreducible repre- 
sentation of S, 8, 8' and g are equivalent, then 
the defining equation’s of S can always be put in 
the form Eq. (35). To find under what conditions 
the irreducible sets 8, 8, and 6" are equivalent we 
go back to the general criterion for deciding 
the equivalence of two representations. Take 
(c:(8), --*, ¢-(8)) to be a linearly independent 
basis for the center of ©. Clearly c,7(8) and 
c.*(8) also belong to the center. Therefore, for a 
particular irreducible representation, 8 and 6 are 
equivalent only if c,(8)=c.7(8) for this repre- 
sentation. Similarly, the equivalence of 8 and 8 
requires that c,(8) =c.7*(8). 

From Eq. (26) it follows in the usual way that 
there exist elements J;;“™ in %, such that 


i™= —Iy,™ and 


[BeS™, Dim ] = ge iBm'” —gimBi™, (36a) 
CLei™, Ima’ ] = —geml in +2 iml en 
+ginl im —ginlim™, (36b) 


where [A, B]=AB-—BA . It is easy to prove that 
the choice of J;,;“ is unique since if J: were 
another set satisfying the same equations then 


Cri™ =I —Te™ 


would commute with 6,“ and would therefore 
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belong to the center of W.. Hence 


o=(Jei™, Tana ]— [Jei™, Jan ] 
=gimC in” —gimCkn™ —gienC im +2 inCem™ 


from which it follows that Cy:“ =0 unless s=2. 
So assuming that s#2, Jy: is uniquely de- 
termined. Put 


Teir=Tpei ++ ++ +Lei™. (37) 
Then Ji:= —Ju and 
(Be, Lim] =£%iBm —ZimB1, (38a) 


[Tes ra) = —Ziml int Ziml kn 
+gknlim—Zinlim. (38b) 


Ix: is also uniquely defined by Eq. (37), since 
from Eq. (23) every Je: is decomposable uniquely 
in the form (37), and Jx:“, in consequence of 
Eqs. (37) and (38), must satisfy Eq. (36) and 
therefore are uniquely determined. J,; is an 
element of @(8) and hence is expressible in the 
form of Eq. (5). Write it as I;:(8). Consider the 
automorphism B—t-'6 of ©. It is then easy to see 
from Eqs. (38) that 


D1 =te™ti"Imn(t-'B) 
satisfy the same equations as (38). Hence 
Txi(tB) =t,"ti"Imn(B) 


so that I,:(8) transforms as a tensor. Applying 
Lemma 1, we immediately get , 


Tei(8) = T(Bx8: —BiBx) 
+(D P**(8)BxQ°-:(8)B8.R...(8)]-, (39) 


where J is a numerical constant, P**’, Q:::, R... are 
tensors constructed from 8 and g only, and the 
dots denote contracted indices. }> denotes a sum 
over similar terms and the minus sign at the end 
means that terms with / and & interchanged are 
to be subtracted. The first term J(8,8:— 8.8) is 
of the simplest possible kind and has been 
purposely separated out. Bhabha’s theory*® is 
based on this simplest substitution 


Iki= I(Bx8:—BiBx). 


By an entirely similar consideration of the 
other two automorphisms of © one concludes 
that J, :(8) = — I:7(8) and Ix:*(8) = Ix1(8). Hence 


Tyi(8) = 1(8181—8 Bx) 
+$LX P***(8)BxQ:::(8)BiR...(B) 
— P™*"(B)BxQ7:::(8)B:R™...(8)]_ (40) 
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where the numerical coefficients in Eq. (40) are 
all real. 


3. EXTENSION TO INDEFINITE GROUND FORMS 


Let us now consider instead of O, other groups 
G of linear transformations with real coefficients 
which leave the indefinite ground form 


guix*x'=0, (41a) 
gei=0, R¥l, ger=+1 0r —1 (41b) 


unchanged. As usual, we use gy: to raise or lower 
the tensor indices. For every co-variant tensor 


Ak,---k, define 
Ak) -+-ky = (4)"Aky-+*ky 


where r is the number of indices k, among 
(ki-++, Rx) such that gkmktm = —1. In this way for 
every co-variant tensor in the x-space one con- 
structs a co-variant tensor in the x-space. In 
particular, gx: is now the same as defined by 
Eq. (4). We regard gi: as the metric tensor of the 
x-space and use it (instead of g,:) to raise tensor 
indices in the x-space. It is then easy to see that 
the tensor corresponding to Abk2---B"" is 
Akyke---B’**"""2 where the dots denote contracted 
indices. Thus every tensor equation can be taken 
over from the x-space to the x-space by just using 
boldface type. However, a real transformation ¢ 
of the x-space corresponds in general to a com- 
plex transformation of the x-space. But it is well 
known that the necessary and sufficient condition 
for the invariance of a set of matrices 8, under G 
is that there should exist matrices [4;= —I 
satisfying Eqs. (38) (now with the altered mean- 
ing of ge:), and further that there be a non- 
singular matrix R such that 


RBR™ = JB, (42) 


where J is the transformation in G corresponding 
to the reversal of the direction of the x-axis. But 
these conditions are equivalent to 


[ Bx, Lim] = Se18m—amG2, (43a) 


CIer, Inn] = —Gimlint+imlin 
+inlim—Zinlim (43b) 
and 
RGR = JB, (43c) 


since J = J. Equations (43), however, are just the 
conditions that the §-matrices be invariant under 
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O. Thus it is sufficient to examine the invariance 
of § under O. This we have already done. Hence 
all our results in section 2 are still valid, provided 
that each tensor Ak--- is replaced by Az,---. On 
going back from the x-space to the x-space we 
have to revert to the normal type for the symbols. 
Therefore, all our conclusions which are expressed 
as tensor equations remain unaltered in form also 
when the new meaning is attributed to gy. 


4. THE SPURS 


The invariance of the B-representation under G 
leads to important restrictions on the spurs of the 
B-matrices and their products. Put 


Sp(Bki- + + Bkn) = Ski---kn, 


where clearly sk;---k, are numerical quantities. 
Under the transformation B-—’ =i8, Eq. (44) 
goes over to 


Sp( Ber’: + * Bken’) = thy"!+ + +thy™S11---tn. 
But from Eq. (3a) 


Sp(Bir’ + - - Bien’) 
= Sp(TBk1: + -BknT—") = Sp(Bkr° - - Ben), 


(44) 


so that 
tky")- . + ten'™S1y+++ly = Ski-* kn. 
From this we deduce as before that 
Sky-+-ky = S(g)k1---kn. (45) 


Since sk---k, cannot be constructed from gx: 
when 2 is odd, it follows that 


Sp(Bkr- - * Bkony1) =0, 


and in particular 


(46) 


Sp(Bx) =0. 
Now if 6, and ;* are equivalent, 
Sp(Bkr° . * Bk n) = Sp(Bkn*- ' - Bk*) = Sp(Bra° + + By) 


i.e., 
S(g)k1-+-kn = S(g)kn---k1- (47) 


Also further, if 8, and @,' and, therefore, 6,* are 
equivalent, the numerical coefficients in Eq. (45) 
are all real. Conversely if 8, are an irreducible set 
of matrices satisfying Eqs. (45) and (47) with 
5(g)k1-:-kn =S*(g)k1---kn, then the 6-representation 
is invariant under G and By, 6:', Bx’, and B;,* are 
all equivalent. To prove this it is sufficient to 
show that if there exists a correspondence 8,—8;’ 
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between two sets of irreducible matrices such that 
Sp(Bir ++ + Bk yn) = Sp(Bir’ ++ +Bkn’), 


then 6, and £;’ are equivalent. Consider the 
algebra 2(B) generated by the matrices B, which 
are the direct sum of ; and 8,’ i.e., 


Bi | 0 





By= 





(48) 











0 | Be’ 





Clearly 2(B) is semi-simple, and therefore either 
it is simple or a direct sum of two or more simple 
algebras. In the first case it is capable of only one 
irreducible representation apart from equiva- 
lence. Since B,—x, and B,—f;’ are both repre- 
sentations of 2{(B) they must be equivalent. On 
the other hand, in the second case 


W(B) =%(B) OA2(B) @---, 


where %;, Ws, --- are all simple. In case 6; and ,' 
are not equivalent, they are equivalent to the 
regular representations of two different simple 
sub-algebras, say %1 and %2, respectively. Con- 
sider the unity quantities e,(B) and e.(B) of 4, 
and Y2, respectively. We have, e:=0. In the 
8-representation @, is represented by the unit 
matrix and é: by the zero matrix, while in the 
8’-representation e,=0 and e:=1. Hence 


Sp(er(8)) # Sp(er(6’)), (49) 


where ¢:(8) and e;(8’) are both obtained from 
é:(B) by replacing B; by 6; and 8;’, respectively. 
Equation (49) therefore contradicts the hy- 
pothesis. Hence 8; and f;’ are equivalent. 


5. THE DIRAC AND THE DUFFIN- 
KEMMER ALGEBRAS 


As an illustration of the above general theory, 
let us consider the algebra generated by the 
Dirac and the Duffin-Kemmer matrices in the 
general case of s-dimensions. It is well known 
that in both of these cases 


Tyr = I(8x81—B Bx) 


where J is a constant. Therefore, as pointed out 
by Bhabha,* any representation of these algebras 
is a representation of the orthogonal group 
O(s+1). Since every representation of the orthog- 
onal group is fully reducible, it follows that‘every 
representation of the Dirac (or the Duffin- 
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Kemmer) algebra is also fully reducible. In par- 
ticular, the regular representations of these alge- 
bras are fully reducible, and so the two algebras 
are semi-simple. 
The Dirac matrices a, are defined by the 
commutation rules*® 
Opt Oe = 2QK1 (50) 
so that 
a,a* =§. (51) 
In any scalar quantity all the indices are con- 
tracted, and on using Eq. (50) the contracted 
indices can always be brought adjacent to each 
other. Therefore, from Eq. (51) every scalar 
quantity is a multiple of J. In the same way it is 
easy to see that every pseudo-scalar quantity 
must be a multiple of 


w =—e*l* Fequey > + + Oks. (52) 

s! 
Now, w commutes or anti-commutes with a, ac- 
cording as s is odd or even. Clearly J and w are 
linearly independent, since a pseudo-scalar can- 
not be linearly dependent on a scalar. Thus a 
basis for the center of the algebra consists of 
(1, w) or J according as s is odd or even. Thus the 
algebra has one or two inequivalent irreducible 
representations according as s is even or odd. 
Hence, for s even, a, a’, a', a* are all equivalent. 
For s odd, a and Ja are inequivalent since 
w*=(—1)#¢-) 0. Also w7 = (—1)*#*“-Yw. Hence a 
and a are equivalent only if }s(s—1) is even, i.e., 
if s is of the form 4n+1. In a particular 
representation w is represented by -two where 
wo=t*-). Therefore, a and a* are equivalent 
only if wo=wo*, i.e., if s=4n+1. a and a’ are 

always equivalent. 
The Duffin-Kemmer matrices are defined by 


BB Bm + BmB iB = 18m + Sm iBrx- (53) 
Put 8.6*=B, so that 


8.B+ BB, =(s+1)£:, (S4a) 

that is, 
8,.B=(s+1—B)B,, (54b) 
B*B 1B. =B1. (55) 


Two contracted indices in a scalar expression can, 
from Eq. (53), be brought either adjacent to each 
other, or immediately to the left and right of a 


* For convenience we take gx: as defined in Eq. (4). 
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8x, so that from Eqs. (54) and (55) every scalar is 
expressible as a polynomial f(B) in B. Further, 
if f(B) belongs to the center, 


f(B) Bi =Bif(s+1—B) =B,f(B). (56) 


Since the algebra, being semi-simple, has a 1 
which must be expressible as a scalar e(8), it 
follows from Eq. (56) on multiplying on the left 
by e*(8), where e*(8)8,=e(8), that 


f(s+1—B) =f(B) =}{f(s+1—B)+f(B)}. 
Put 


s+1 
}(s+1)—B=A. and (—-a) = (A). 


Then every scalar quantity of the center is 
of the form }¢(A)+4¢(—A) =y(A*) where ¥(A*) 
is a polynomial in A*. Now A*=}(s+1)? 
— B(s+1—B). Put 


6=B(s+1—B). (57) 


Then every scalar quantity in the center is ex- 
pressible as a polynomial in @. 

In a similar way, it follows that every pseudo- 
scalar quantity is of the form wf(B), where f is a 
polynomial in B and 


w = eft FeBe,+ + + Bhs. (58) 
Now 


Br yw = eke | BE ?Bi-g- - + Bka — Bk BeBe Bes: * - Bk, 
+ Bk BkoBesBk iBks- - - Bk 
+++++(—1)*BkiBko- + - BegBer} 


where there is no summation over the index &). 
Using Eq. (53) and remembering that ki, ke, ---, Re 
are all different, we get 








s+1}]. 
arna—| ; Jerse ° *Bks, (59a) 
where [m ]denotes the integral part of n. Similarly, 
s+1 
wBky -| 2 Jerse + + Bk Bk. (S9b) 


Remembering that 8,78;=8;(1 —8,?) for k#1, one 
finds from Eq. (59) that 





Bio — wBx 
0 if s is odd, (60a) 
= 4fst+l 
; Jette -- Bk, if s iseven. (60b) 
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Also, from Eq. (59b), 


s+l1 . 
wk? = |—|or-sn + + Bk Bk. 


Hence 


1 
wB=> wse:t=|—. (61) 


ki 


Therefore any pseudo-scalar quantity wf(B) is 
really a numerical multiple of w. Now it is clear 
irom Eqs. (60) that w belongs to the center only if 
s is odd. Hence the center contains pseudo-scalar 
elements only if s is odd and then they are nu- 
merical multiples of w. 

Consider the expression 


B, => (k1---kr) Bky?- + + Bey? 


where the summation is to be taken over all 
possible values of ki, ---, &, such that no two of 
the k’s are the same. Notice that 


B,B= LD hi -+ kr) Bk? . + Bee? Dokryt Bkr+1”. 


Now 8,? and 8,? always commute, and 6,‘ =8;,’. 
Hence the right side is 


Dd (hr + -kr41) Ber? + + Bee?Bk py? 
+7 (k1-+-kr) Bkx®> + + Bk,?, 


(62) 


so that 
Brat =B,(B-r). 


But B,=B and, therefore, 
B,=B(B-1)---(B—r+1). 


+, key cannot all be 


(63) 


Clearly, B,4:=0 since ki, -- 
different. Hence 


B(B—1)---(B—s) =0. 


Notice that B=B,, Bo, ---, B, are all line- 
arly independent and so B, B’*, ---, B* are 
also linearly independent over §. Put f(B) 
= B(B—1)---(B—s). Then f(B)=0 and so be- 
longs to the center. Therefore, f(B) = f(s+1—B) 
=0; i.e., 


(s+1—B)(s—B)---(1-—B) 
=(—1)**"(B—1)---(B-s) 


(64a) 


X(B—(s+1))=0. (64b) 
Hence from Eqs. (64a) and (64b), 
(B—1)(B—2)---(B—s) =0. (65) 


Equation (65) is the minimum equation for B, 
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since otherwise B, B*, ---, B* would not be 
linearly independent. From Eqs. (65) and (57) it 
follows that the minimum equation for @ is 


(@—1-s)(@—2(s—1))--- 
x (6@—[4(s +1) 4(s+2)]) =0. (66) 


Hence 1, 0, 6’, ---, 0[3(s+1)]—1 forma linearly 
independent complete scalar set in the center. 
When s is even, this set is a basis for the center. 
On the other hand, when s is odd, w must be added 
to it to obtain a complete basis. Therefore the 
dimension of the center over § is [$(s+1)] or 
[$(s+1)]+1, according as s is even or odd. 
Hence the algebra has $s or }(s+3) irreducible 
inequivalent representation, according as s jis 
even or odd. 

From Eq. (61) it follows that w is zero in all 
representations except those in which B can have 
the eigenvalue [4(s+1) ], i-e., in which @ has the 
eigenvalue [3(s+1) ][$(s+2) ]. But @ is in the 


center and so 
s+17[s+2 
aes ea 
2 2 


for this representation. 
Proceeding as in Eq. (59a), one easily finds 
that 
s+1 1s 
BkoBk iw = |= Jen seks 
2 2 
X (1 —Bk1?)Bk2*Bks---Bks, (Rixk), 


and, similarly, 


s+l1 TI's 
Bk Bks—1** -bw-|——]|-]. oe 
2 2 


2 
[Jerre — Bks—1”) Bks—2"(1 —Bks_s*)-**, 


where hi, ---, &, are all different. Therefore 


w? = (— 1tmee-of os -] 
2 2 


X P(k1...k)Bks2Bky—2?* - - (1 —Bks—12) (1 —Bks-s*)***5 
(67) 


where P;---%,) denotes a sum over all permu- 
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tations of (1, 2, ---s). It is easy to see from deduce that 


Eqs. (67) and (63) that 


socom fT -L 
Hed 


5 al 
X | Bacto! |= Busan Brye+s)) 


GET) 


3! 





+++++(—1)49B,}. (68) 


Using Eqs. (63) and (61), it follows that 


encom EL 


Now w’? is a scalar and hence a polynomial in B, 
as is also clear from Eq. (68). Hence w? can be 
written in the form 


w* = P(6)+Q(0)-B 


where P and Q are polynomials in @. Since w? 
vanishes for all representations except those for 
which 6=6o, it follows that 


w=9(6){p+¢B}, (70) 


where 


tomer ZPD) 


In view of Eq. (66), » and q can be taken to be 
constants. From Eqs. (61), (70), and (66) one 


while on multiplying Eq. (70) by (s+1—B) one 
gets 
\— 
w? 
2 


From Eqs. (69), (70), (71), and (72) it is easy to 





|-e@te6+1)+06-pa}. (72) 
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x|—+1~81, (seven), (73a) 


w’?=-¢ 
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{. yim) 
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X (3(s—1)!)4 (s odd). (73b) 
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Now notice that it follows from Eq. (59a) that 








rs+17 
BkywBk, = a etka: -keBe Bee » «Be Big 
Ps 15 
=| —— Jett Beeo- + Bk (1 — Bki?) Bey 
C2 J 
=0, (s even). (74) 


It will presently be shown in Section 6 that 
Eq. (74) necessarily implies 


8x08: +B8w8.=0, (s even). (75) 


Now consider w8,w. It is a vector and hence must 
be expressible in the form 6,f(B) where f(B) is a 
polynomial in B. Then 


8.208 .w = 8 f(B) =6if(B) = wBw. 


But the left side is zero from Eq. (74), when s is 
even. Hence 
wB.w=0, (s even). (76) 


Further it is easy to see from Eq. (61) that 
Bk;?w+whky2=w, (Ss even). 

As shown in Section 6, this implies that 

(8:8: +B Bx)w+w(BcBit+BiBx)=2wgu, (s even). 


But w is a pseudo-scalar and therefore commutes 
with J, i.e., 


(8x81 —B Bx)w —w(8x8:—B Bx) =0. 
Adding the two equations, one gets 
BiBiwtw8 i 8.=wgn, (Ss even). (77) 
Moreover, Eq. (73a) gives 


¢(8) 
Baw? + w'Be= (—1)UO-M1(§ s1)¢ By. (78) 
$(9o) 
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Then wo’=wo. The only irreducible representa- 


tions in which w»+0 are those in which @= 6». In 
these representations 


, 4s+1-—B (seven), (79a) 

@ = 
’ 1 (s odd). (79b) 
In view of Eqs. (75), (76), (77), and (78), one can 
put 8,41 =wo, (s even) for the representation @= 0, 
and allow the indices to run from J to s+1 in 


Eq. (53). Also notice that for this representation 
the matrices y, defined by 


Yr =Biwot wok (80) 


satisfy the same commutation rules as Eq. (53). 
In fact 





ao= 


R= (1+ wo — wo) Be(1 + wo — wo”) ; (81) 


and since (1+wo—w»?)?=1, y and 8 are equiva- 
lent. Since 1-+-wo— wo? is a pseudo-scalar quantity 
it commutes with J,;. Also clearly it commutes 
with B and wo; therefore 


VEYI-V1Vk =BiBi—BiBe, yv*ve=B*Bx, 
firs Kany, + Vky = efi kaBy, » + Bhs. 


Since @= 6* = 67 it follows that for even s, any 
irreducible f-representation is equivalent to 
Jp, 8*, 6, B'. For s odd, 6, JS, 6*, 6, B' are 
equivalent for the representations for which 
6605, w=0. In the other two 6= 4, 


w= bille—D14 (54-1) (4(s—1)!)?. 


Clearly here 6 and J are inequivalent. 8 and 6* 
are equivalent only when s is of the form 4n+1. 
The same holds for 8 and §*, while 8 and 6 are 
always equivalent. 


6. PHYSICAL APPLICATION 


We now come to the physical application of the 
above abstract theory. The tensor indices now 
run from 0 to 3 only and goo= —gi1= —ge2 
= —g33;=1. Consider the wave equation (1). Let 
the minimum equation for 8» be 


(Bo)*+a1(Bo) -_ +a2(Bo) o-3 4 sat +d, =0, (82) 
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with a’s in §. Write it in the form 


(Bo)" + go0d2(Bo)"~*+ (go0)*a4(Bo) "4+ - - - 
+ {a1(Bo)"-"+-gooas(Bo)"-* +... } =0, 


Since 6 and #8 are equivalent, ¢ now being a 
Lorentz transformation, it follows that 8,’ =t'B, 
has the same minimum equation. Put ty =y, 
Then 


y* yx =to*to'gnt=Z00 = 1 
and Bo’ = y*8,. Therefore 


{ (y 1g 1) ” +azy"y"*gk ike (y 1g 1) n—2 
+-a4(y**y**geik2)*(y'B1)"-4+ ++ +} 
+ {ai(y'B1) "> +asyiy"egtita(y'Bi)"*+- +} <0, 


Now the y’s are arbitrary except for Eq. (83), 
Hence put 


(83) 


z' g! 


EEE 


m (gmaz"2")? Z ; 


where the 2’s are independent indeterminates, 
Then on multiplying by Z* one finds 


{ (24B1)"+ae2"!z**gk ke(2'B ) m2 4 eee } 
+Z{ai(z'B,)"""+---}=0. (84) 


This must be an identity in z. Therefore from the 
irrationality of Z in z it follows that 


a;(2'B,)""!+---=0. 
Put 2°=1, z'=0, 1~0. Then 
a1(Bo)*""*+--- =0. 


If a1, a3, *** are not all zero, this is an equation, 
of a degree lower than the minimum equation. 
Since this contradicts the definition of the mini- 
mum equation, it follows that a; =a;=a5=---=0. 
Hence only the terms in the first bracket in 
Eq. (84) remain. Equating coefficients one gets 


Pky -++kn)(Bkr° ° * Ben +GogkikoBks- + + Bin 


+ dagkkogksksBkg* > -Bkent+---)=0, (85) 


where P¢;---k,) denotes a sum over all permuta- 
tions of ki---k,». From Eq. (85) it follows that 


{ (0'B1)"+a20*0,(0'B,)"-? 
+a4(9*dx)?(9'B1)"-*+ +> -jy=O. 
On making use of Eq. (1), this can be written as 


{ (tx)"+a2(ix)"-*0*d; 


+a4(tx)"-*(d*0x)*+---}y=0. (86) 
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Equation (86) is the differential equation (with 
numerical coefficients) of the lowest order which 
the most general solution ¥ of Eq. (1) would 
satisfy if no subsidiary conditions were imposed. 
Now, if we demand that (86) be the usual second- 
order equation, we must have"! a4=as= --- =0 
and a2= —1. Equation (85) then becomes 


Pier: n)(Bk1* + * Ben — gkikeBks * *Bkn) =O. (87) 
Equation (87) is completely equivalent to 
(2'B:)"=2*2,(2'B.)"-’, (88) 


where 2's are indeterminates. Equation (88) is the 
natural generalization of the well-known cases of 
the Dirac and the Duffin-Kemmer matrices, for 
which 
('B,)? =s*2, 
and 
(2'B1)* =2*z,(2'B), 


respectively. 

Since Eq. (87) is already in a proper tensor 
form, one might at first sight think that the 
invariance of an irreducible representation under 
the proper Lorentz group is insured. This, how- 
ever, is in general not the case, since for n23 
Eq. (87) does not generate a finite algebra. Some 
other stronger tensor condition compatible with 
Eq. (87) is needed to make the algebra finite. An 
example of such a condition for n =3 is provided 
by the Duffin commutation rules 


BxB Bm + BmB Bx — Ze ibm — ZmiBx =O. 


In the Dirac and the Duffin-Kemmer case it is 
possible to choose a representation in which #’s are 
Hermitian (i.e., 8) =8o' and 6, = —8,', k=1, 2, 3). 
But this is clearly impossible when n> 3, since if 
8 is Hermitian it can be brought to the diagonal 
from. Now, its minimum equation is 


Bo”~*(8o?—1) =0, (89) 


and hence its eigenvalues are +1, 0. But then, 


"It is assumed here that x0, so that one can divide 
out (ix)*-, 
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since Bo is diagonal, it clearly satisfies the 
equation 


Bo(Bo? — 1) =0. 
Hence n& 3. 

I have not been able to establish the existence 
of an irreducible relativistically invariant set 6; 
satisfying Eq. (87) for n>3, but not for n <3. 
All the obvious commutation rules which are 
compatible with Eq. (87) and with the general 
restrictions laid down in Section 3, and which 
make the resulting algebra finite, either turn out 
to be inconsistent or are found, on closer in- 
vestigation, to degenerate into the already known 
cases of the Dirac or the Duffin-Kemmer rules. 
On the other hand, there seems to be no obvious 
way of proving that such a representation is 
impossible for n>3. If it be true that no such 
8-representation could exist for »>3, our results 
should be interpreted to mean that no higher 
spins are possible. 

It is comparatively easy to find an explicit 
expression for 8; with »=4; e.g., 


1 
hen iale)Xots- 1Xax(})} 


where a,(#) and a,(4) are matrices constructed in 
the way proposed by Bhabha,’ which satisfy the 


minimum equations 
(ao? —4)(ao?—9/4)=0, and (ao?—})=0, 
respectively. w is uniquely** defined by w* = 1 and 
an(})w+wax($) =0. 


a,(4) are really } times the usual Dirac matrices. 
X denotes Kronecker product. It is easy to see 
that this representation is relativistically in- 
variant (including reflection) and is equivalent to 
its Hermitian conjugate. But since it is not 
certain that it is irreducible, one cannot rule out 
the possibility that the irreducible part may 
satisfy an equation such as (89) with m3. 


* Since the manuscript was sent to the press, one such 
representation has been found corresponding to n=4, 
*% Apart from ambiguity in sign, which is immaterial. 
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Absorption Line Width in the Rotation Spectrum of Atmospheric Water Vapor* 


ARTHUR ADEL 
Mc Math-Hulbert Observatory of the University of Michigan Ann Arbor, Michigan 


(Received February 28, 1947) 


Grating measurements at oblique incidence are employed as a means of obtaining the narrow 
slits essential to the direct observation of line width. Correction for the finite width of slit is 
small. It is accomplished by empirical extrapolation to zero slit in one instance, and by nu- 
merical estimation of slit influence in another. The absorption coefficient half-width at half- 
height for the atmospheric water vapor lines at 18.64 mu and 15.99 mu is found to be 0.12 cm=! 


and 0.11 cm™, respectively. 





R a spectrometer of unlimited resolving 
power the fractional transmission of radia- 
tion through an absorption line is given by the 
expression, 
Teo) +t, 


and the fractional absorption by A=1—T. By 
the width of the line is meant the width of the 
absorption coefficient k(v). In the infra-red, 
broadening by thermal collisions is the principal 
agent in shaping an absorption line, and the 
absorption coefficient has the Lorentz dispersion 
form 





J Do 
k(v) =— ’ 
T (v— vo)? +De? 


where vo marks the center of the line, and 


J= k(v)dv 


is the total intensity of the line. The distance 
from vo to the point at which k(v) falls to half 
its maximum value is Do, the half-width of the 
line, that is, half the width at half the height of 


Sut —_ —_— WIDTH 


uv 


Fic. 1. The intensity distribution within the spectral 
interval transmitted by the exit slit of a spectrometer 
with equal entrance and exit slits. 


*This paper is based on work done for the Office of 
Scientific Research and Development under Contract 
OEMsr-1361 with the University of Michigan. 
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the curve of absorption coefficient. J and Dy are 
important in the problem of radiative transfer. 
In an absorption line of moderate depth the rate 
at which energy is absorbed is proportional to 


(JDo)', 


where ¢ is the length of path. 

For the purpose of measuring widths, the 
curve of absorption coefficient can always be 
obtained from the absorption line following cor- 
rection of the latter for the finite slits of the 
spectrometer. This correction can be made 
directly for a line of any strength by using a 
graduated series of slits and extrapolating the 
curve of apparent line width versus slit width 
to zero slit width. The correction can also be 
made indirectly by application of the simple 
theory of slit influence to a line of moderate 
depth.? Both methods will be applied to atmos- 
pheric water vapor—the empirical to the line 
at 18.64 mu, the numerical to the line at 
15.99 mu. 

Infra-red lines are generally exceedingly nar- 
row. It has been shown, for example, that the 
directly observed absorption line half-width in 
N,O, at a pressure of 30 cm Hg, even prior to 
correction for slit width is but 0.09 cm~'.? The 
need for narrow slits is evident, since one may 
hope to explore a line significantly only if the 
slit is narrower than the line. 

The spectral width of a slit depends more im- 


‘portantly upon the available dispersion than 


upon actual slit opening. The intensity distribu- 
tion within the spectral interval transmitted by 


1D. M. Dennison, Phys. Rev. 31, 503 (1928). 
2 Arthur Adel and E. F. Barker, Rev. Mod. Phys. 16, 
236 (1944). 











or* 


Dy are 
ansfer, 
he rate 
al to 


is, the 
ays be 
ng cor- 
of the 

made 
ising a 
ng the 
width 
ilso be 
simple 
»derate 
atmos- 
he line 
line at 


ly nar- 
iat the 
idth in 
rior to 
12 The 
1e may 


if the 


ore im- 
1 than 
stribu- 
‘ted by 


*hys. 16, 











SLIT WwiOTH cw 








36 PA 10.4 . 16o 19.5 204 MY x 
! ! | ! 
10 30 90° _ 9) 








Fic. 2. Effective slit width as a function of wave-length 
or grating angle for the spectrometer \, = 20.764 sin@, with 
2400 line grating and slits of angular width one minute 
of arc. 


the exit slit of a spectrometer with equal entrance 
and exit slits is shown in Fig. 1. This is the fre- 
quency pattern received by the thermopile at 
the frequency setting v;. The width of the slit 
is the width of the slit distribution function at 
half-height. The functional dependence of spec- 
tral slit width on slit opening and dispersion can 
be obtained from the spectrometer equation 


A=2a cos}¢ sind, 


where a is the grating constant, ¢ is the angle 
between the incident and diffracted beams, and 
8 is the angle between the grating normal and 
the bisector of ¢. 

For entrance and exit slits of equal angular 
width dg the spectral slit width is given by 


10*-do 
———(cotd —tan}$¢) cm, 


di,= 
2K, sind 


where 
K,=2a cos} gu. 


dé, must not be confused with dic, the spectral 
interval which sweeps past a point in the exit 
slit when the grating turns through a small angle 
di=dy; for the grating sweeps the spectrum 
across the exit slit faster than would a mirror 


ABSORPTION LINE WIDTHS 


a ray of light, that is, 
10‘dd 
K, sind 





cot’ > 2dF_, 


ic= 
when dd=dy, the slit width. di, may also be 
written 


10‘-dy cos(d+4¢) 


2K, sind sind ‘cost 





dv, 


from which it is clear that di,=0 at grazing 
incidence, that is, when 3+}$g¢=$2, and re- 
mains very small at oblique incidence. 

In the experiment described in this paper a 
grating with 2400 lines per inch was used, thus 
making 2a=127/6 mu. From the spectrometer 
constant, K,=20.764 mu, it is found that 
$¢=11° 20’. Figure 2 is a plot of di, versus 8 
(or A) for a grating with 2400 lines per inch and 
an angular slit width of one minute of arc, 
according to the equation 


cotd — 0.2004 
cm, 





di,=0.07 
sind 


The technique of oblique incidence has been 
applied to the atmospheric water vapor lines at 
15.991 mu, 16.667 mu, 18.018 mu (74—6o), and 
18.64 mu (6;—5;). All are extremely narrow and 
essentially of the same width. Only the two at 
15.991 mu and 18.64 mu will be treated in detail. 
The line at 18.64 mu was examined with a 
variety of slit widths ranging from 0.62 cm™ to 
0.061 cm. The reduction appears in Fig. 3, 
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Fic. 3. Observed dependence of absorption line width 
on spectrometer slit width. 
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5 Line half-width is here plotted against slit half-width. 
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Fic. 4. The atmospheric absorption line at 18.644 and 
its absorption coefficient. 


which has the form indicated by theoretical 
considerations.‘ The extrapolation to zero slit 
is short and can certainly lead to a value no 
greater than 0.15 cm™ for the half-width at 
half-height, at the observed fractional absorption 
of 56 percent. The line, and the curve of k(v):t 
deduced from it, appear in Fig. 4. k(v)-¢ and k(y) 
possess the same half-width, since ¢ is a constant 
factor. The half-width at half-height of R(v) is 
0.12 cm, 

It also seems worth while to correct for finite 
slits by computing the effect of slit-width on 
apparent line width, in contrast to the empirical 
extrapolation just concluded. The elementary 
theory of slit influence has been applied to the 
problem of line width*® in N,O. Figure 5 is 
taken from this work. Here the apparent half- 
width of the line, in units of true half-width, is 
plotted against the apparent half-width, in units 
of slit half-width. The abscissa is an observable. 
Thus D/Dy may be read from the curve, D 
being known, Dy follows. It is assumed in this 
analysis that the line is weak so that line half- 
width and absorption coefficient half-width are 
sensibly equal. The atmospheric water vapor 


* See Fig. 6, reference 2. 
5 The author takes this opportunity to correct a small 
error in the work on N,O. dic rather than dé, was used to 
compute slit width, and, as a result, the line half-width was 
overcorrected to 0.07 cm™. By using the true slit half- 
width of 0.02 cm, the apparent line half-width of 0.09 
cm for N;O is reduced to 0.08 cm, true line half-width. 
* For an observed central absorption of 39 percent the 
line is of moderate strength rather than weak, and the 
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Fic. 5. Infra-red ors = a line width in units of 
absorption coefficient width as a function of absorption 
line width in units of slit width. 
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line at 15.991 mu was observed to have a half- 
width at half-height of 0.145 cm™, at a time 
when the maximum absorption was 39 percent, 
and the slit half-width was 0.075 cm. Since 
D/S=2, D/Do=1.35, and Do=0.11 cm, in 
substantial agreement with the value obtained 
for the line at 18.64 mu. 

Since this work was done, in the summer of 
1944,’ several investigators have determined the 
half-width of the water vapor line (6.5;—5_,) at 
41.35 cm in the microwave region. See, for | 
example, G. E. Becker and S. H. Autler,® who 
found values of 0.087 cm=! and 0.107 cm for | 
low and high density of water vapor, respec- | 
tively, in good agreement with the results given | 
in this paper. | 


computed half-width will be conservative, that is, slightly | 
too large. 

7 Div. 14, NDRC. Report No. 320, October 10, 1944; 
and an unpublished letter to Professor G. E. Uhlenbeck, 
dated February 1, 1945. 
( 045) E. Becker and S. H. Autler, Phys. Rev. 70, 300 
1946). 
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Finite Elastic Strain of Cubic Crystals* 


FRANCIS BIRCH 
Harvard University, Cambridge, Massachusetts 


(Received February 24, 1947) 


Murnaghan’s theory of finite strain is developed for a medium of cubic symmetry subjected 
to finite hydrostatic compression, plus an arbitrary homogeneous infinitesimal strain. The free 
energy is developed for cubic symmetry to include terms of the third order in the strain compo- 
nents. The effect of pressure upon the second-order elastic constants is found and compared 
with experiment, with particular reference to the compressibility ; the pressure-volume relation 
in several approximations is compared with the measurements to 100,000 kg/cm*. The simplest 
approximation is shown to give a satisfactory account of most of the experimental data. The 
results are also compared with some of the calculations based on Born's lattice theory. 





1. INTRODUCTION 


URNAGHAN! has published an exact 
treatment of finite deformation of elastic 
solids, with application to several questions 
involving isotropic media. The theory is not 
restricted to the case of isotropy, and it will be 
developed below for the case of cubic symmetry. 
Finite elastic strain is in practice limited to 
finite hydrostatic strain plus superimposed in- 
finitesimal strain of an arbitrary type. Excluding 
such special materials as rubber, non-hydrostatic 
strains large enough to require treatment as 
“finite” are generally beyond the elastic limit. 
On the other hand, hydrostatic strain of what- 
ever magnitude is “‘elastic’’ and may be accom- 
panied by elastic infinitesimal strain of any kind. 
The hydrostatic strain affects the response of the 
medium to the infinitesimal additional strain, 
giving rise to an effect of pressure upon the 
elastic constants. Such a strain is of fundamental 
interest in many problems concerning the interior 
of the earth, to a few of which the theory has 
been applied.? A large body of experimental data, 
almost entirely the work of Bridgman, may be 
used to find the change of compressibility with 
pressure ; there are also a few determinations of 
the effect of pressure upon the rigidity of quasi- 
isotropic materials. 
Born’ has treated this kind of deformation 


*Paper No. 100 published under the auspices of the 
Committee on reuhaantal Geology and Geophysics, and 
the Division of Geological Sciences at Harvard University. 

1F, D. Murnaghan, Am. J. Math. 49, 235 (1937). 

*F, Birch, Bull. Seis. Soc. Am. 29, 463 (1939); K. E. 
Bullen, ibid. 30, 235 (1940). 

* Max Born, J. Chem. Phys. 7, 591 (1939); also a series 
of papers in Proc. Camb. Phil. Soc. 


from the point of view of lattice theory, and has 
carried out with several collaborators detailed 
computations for various monatomic cubic lat- 
tices. The results are of great interest, showing 
plausible pressure and temperature coefficients 
for the elastic constants. Born appears to have 
accepted as general the theorem of the classical 
theory of elasticity, which holds that the stress 
components are the partial derivatives of the 
free energy density with respect to the corre- 
sponding strain components. Murnaghan has 
shown that this is valid only for infinitesimal 
strain, and in consequence, some of the calcula- 
tions according to Born’s method are subject to 
a correction which increases with the amount of 
compression. 

Another paper on this general topic has been 
published by Fiirth,t who combines certain re- 
sults of Born and Misra® concerning the third- 
order terms of the elastic potential with an 
incorrect application of Murnaghan’s stress- 
strain relation. Except in the special case of 
hydrostatic pressure, Fiirth’s calculations are 
consequently wrong, and with the same excep- 
tion, they deal with unrealizable examples of 
finite elastic strain. 

In the following sections, the theory is devel- 
oped for a medium of cubic symmetry subjected 
to a strain consisting of (1) a hydrostatic com- 
pression of any amount, plus (2) a homogeneous 
infinitesimal strain of general type. The exact 
stress-strain relations for this case are found and 


*R. Firth, Proc. Roy. Soc. A180, 285 (1942). 
§ Max Born and R. D. Misra, Proc. Camb. Phil. Soc. 
36, 466 (1940). 


809 











810 


decomposed into a hydrostatic pressure plus a 
general stress related to the infinitesimal strain. 
A general form for the free energy is given, to 
the third order in the strain components, and 
used to find the variation of the elastic constants 
with pressure. Particular attention is given to 
the effect of pressure on compressibility and on 
volume, and to the relevant experimental data. 
Since there apparently exists some misunder- 
standing concerning the application of Murna- 
ghan’s theory, most of the developments are 
carried along according to both of the alternative 
methods which Murnaghan has formulated ; this 
naturally lengthens the treatment, but it is 
perhaps worth while to show clearly that no 
ambiguity arises. The particular case of hydro- 
static pressure could have been dealt with more 
directly, and indeed, since cubic crystals are 
isotropic with respect to such a stress, this case 
has already been covered in the studies of 
isotropic media."* The general theory is given 
for the sake of comparison with Born’s lattice 
theory, and with experimental data which should 
be obtainable. 


2. THE STRAIN 


Throughout the following, rectangular Car- 
tesian coordinates are employed, the axes taken 
to coincide with the axes of fourfold symmetry ; 
the same coordinate system is used for the 
description of the initial and final configurations. 
The usual notation of Cartesian tensors is used, 
repetition of an index signifying summation over 
the values 1, 2, 3 for this index. The coordinates 
of a point before strain are denoted by (a1, a2, a3), 
after strain by (x1, x2, x3). These points determine 
a displacement vector with components (x1—a1), 
(x2—@2), (xs—@s). 

We now consider a special strain composed of 
two parts: (1) a hydrostatic compression by 
which every line in the crystal is shortened by 
the factor (1—a), where a, intrinsically positive, 
may have any value from 0 to 1; (2) a subsequent 
homogeneous strain expressed in terms of six 
coefficients 8,,, which will be treated as infini- 
tesimals, whereas all powers of a are retained. 
By the hydrostatic strain, a point initially at 
(a1, @2, @3) is brought to the point (x1°, x2°, xs°), 


*F. Birch, J. App. Phys. 9, 279 (1938). 
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given by x,°=(1—a)ay, with p=1, 2, 3, The 
additional displacements resulting from the 
strain B,, will be 


uy’ =Bi1X1°+Bi2r2"+-Bisx;”, 
Us’ = B21X1°+B22%2" + Bosx;, (1) 
U3’ = B3iX1°+Bs2X2° +Bs3x;3°, 


with 8,,=8,,. The final coordinates are x,=x,0 
+u,’, or 


x1 = (1—a@)[(1+811)a1 +B 1202+ i905], 
X%2= (1—a)[B2101+(1+B22)a2+Beeds], (2) 
x3 = (1 —a) [3101 +Bs2d2+(1+ss)a]. 
From these we find the partial derivatives: 
Ax,/Aa, = (1—a) (5re+Bre); 
04,/ OX. = (Sre— Bre) /(1— a); (3) 
r,s=1, 2, 3. 


Here 6,, is equal to 1 for r=s, to 0 for r#s; 
squares of 8,, have been neglected. 

Putting po for the initial unstrained density, 
p for the density in the strained state, with V, 
and V for the corresponding specific volumes, 
we have the mathematical identities, 








Oa; Oa; Oa, 
Ox, OXe OX3 
0a2 Ode 0a2 
p/ po= Vo/ V=|— eager ines 
Ox, OXe OX3 
0d3 0a3 0a3 
Ox) OXe OX3 

Ox, OX, AX, 

0a, 0a2 0a3 

OX2 OX OX2 

ait — =e 

0a, 0a2 0a3 

OXs OX3 OX3 

0a, 0a2 Oa3 








Using the relations (3) and neglecting squares 
and higher powers of the 8,,, we obtain 


Vo/ V=(1—B811—B22—Bss)/(1—a@)*. ~— (8) 








(1) 


x, =x,! 


J, 
J, (2) 


wr rs; 
lensity, 


vith Vo 
olumes, 


-|. (4) 





squares 


(5) 








For all #’s equal to zero, we have the exact 
formula, V = Vo(1—a)*. 

Murnaghan has shown how finite strain may 
be treated rigorously from either of two view- 
points; by analogy with the similar situation in 
hydrodynamics these are termed the “La- 
grangian” and the “Eulerian.” In the ‘La- 
grangian” scheme, the strain is described with 
the initial or unstrained coordinates (ay) as the 
independent variables ; in the “Eulerian” scheme, 
the coordinates (xp) of the strained state are 
taken as independent. The two alternative sets 
of strain components are then, using Murna- 
ghan’s notation, 


the Lagrangian strain components, 
OXp OX» 


tre = 4 a, a in) ’ (6) 


or the Eulerian strain components, 


0a, 0a 
oe oe a4 Y (7) 
OX, OX, 





Except for a factor of $ for rs, the 9,, are 
identical with the strain components given by 
Love for finite strain. They are also equal to 
4t,,, where £,, are the strain components used 
by Born and Misra and by Fiirth. For infini- 
tesimal strains, the 7,, and €,, become identical 
and (except for the factor } for rs) equal to 
the strain components of the classical theory. 

For the particular strain defined by (2), these 
become 


tre = ben +Bre(1 +2n), 
and (8) 
€ra = 5.€+Brs(1 = 2e), 


where 9 and ¢ are defined by 
(1—a)*=(1+2n)!=(1—2¢)-!. (9) 


The relations between the quantities a, 7 and « 
and the density for hydrostatic compression, 
are tabulated for convenience in Table 1. 


3. STRESS-STRAIN RELATIONS 


Murnaghan’s fundamental contribution is the 
discovery that for finite, elastic, isothermal 
strain, the exact relation between the stress 
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TABLE I. Com tive values of ¢, a, and » for 
hydrostatic compression. 

p/ pe V/Ve —e a -1 
1.0 i. . = 0 0 
1.1 .909 0.0328 0.0313 0.0308 
1.2 833 .0646 .0590 .0572 
1.3 .769 .0955 .0837 .0802 
1.4 .714 -1257 .1061 -1005 
1.5 .667 1551 .1264 .1185 
2.0 .500 .2937 .2063 1850 
3.0 .333 -5400 .3067 .2596 
4.0 .250 .7599 .3700 .3016 
5.0 .200 .9620 4152 .3290 
x 0 2 1.0000 .5000 








components 7,, and the strain is of the form 





Od OX, OX, 0 a, Aa, 
T= p— ——=—2p———; (10) 
One Ip OAg Ojpe OX OX, 
r,s, Pp, g=1, 2, 3. 


Here jrs=(0a,/0x,)(0a,/8x,), and tension is 
taken positive. In developing these expressions, 
all nine partial derivatives are to be formed; in 
this procedure we ignore the symmetry relations, 
Np¢=Nep- Furthermore, it is supposed that the 
rotation has been removed from the strain, so 
that the matrix 0(x)/d(a) is symmetric. 

The function ¢ is the Helmholtz free energy 
(=E-—TS) per unit mass; for a given material, 
it is assumed to depend only on the strain 
components 7,, and the temperature. The second 
equality in (10) is the Eulerian form of the first; 
the two are identically equal. ¢ is equivalent to 
the function A of Born, who, while using finite 
strain components equivalent to 7,,, has used 
the “‘infinitesimal’’ form of (10) for finding the 
stress. 

In an isotropic medium, ¢ depends upon the 
strain components by way of the strain invariants 
(for arbitrary rotation of coordinate axes) of the 
first, second and third orders and their various 
products. In this simplest case, the Eulerian 
form becomes 


do dd 
T+. = p| ——2e-5— }. (11) 


Gre O€ns 


Even in this case the stress matrix does not in 
general possess the same principal axes as the 
strain matrix 7,.,, whereas the principal axes of 
the matrices T,, and ¢,, coincide. Fiirth’s pro- 
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cedure amounts to taking the expression (11) 
for the isotropic medium, substituting 7,, for ¢,., 
and finally applying this result, incorrect even 
for the isotropic case, to the problem of cubic 
symmetry. Only for the degenerate case of 
hydrostatic compression does this give a correct 
result. 


4. THE FREE ENERGY FOR CUBIC SYMMETRY 


In a crystalline medium, the form of the free 
energy ¢ must conform with the symmetry of 
the crystal; in other words, ¢ must be invariant 
for rotations of the coordinate axes corresponding 
to the covering operations. The forms of ¢ 
appropriate for the various classes of crystals 
have been developed for terms of the second 
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order in the strain components, and Constitute 
the basis of the classical theory of crystal elas. 
ticity. For the present work, it is desirable to 
extend the development to include terms of the 
third order in the strain components. The details 
of this process are given in the Appendix, where 
it is found that the third-order terms involye 6 
new elastic constants for the most symmetrical 
classes of cubic crystals, 8 for the cubic clasges 
of lower symmetry, and 56 for triclinic symmetry. 
Most actual cubic crystals belong to the class of 
highest symmetry and the further developments 
pertain to this case; the required modification 
for the other case is slight and easily made. 
The free energy, including third-order terms 
but omitting terms independent of strain, is thus: 


pod = (11/2) (911? + 22? + 9337) +C12( 911922 + 922933+ 139711) +C44( m2” + 9212+ N28? + 1327+ ms? + N31”) 
+ Cy1a(911°+ 922° + 33°) + Cire { 911°( 22+ 933) + 9227(911 + 933) + 938°( 911+ N22) } 
+ Ci2s9119227383 + Caso( 012928981 + 021713782) + (Ciaa/2) { 911923" + 032") + 922(m137+ 031°) 
+ ss(m12? + 921) } + (Cie6/2) { (11 +022) (127+ 921) + (922+ 033) (nes? + 0327) 


For ease in forming the relations (11) the 
equalities 77.=.r have been ignored in (12). 
Rather than ¢, the product pod has been ex- 
panded so that the second-order coefficients will 
equal those of Voigt, which are derived from 
the energy per unit volume. The C’s with triple 
indices are the new ‘“‘third-order’’ coefficients. 
For a given material, the c’s and C’s are functions 
only of the temperature. This differs from 
Born’s development of the free energy in which 
the c’s are functions of both pressure and 
temperature. 

This expansion is based on the assumption 
that there is always a state of zero strain under 
no external forces for any temperature at which 
the crystal can exist. Energy dependent upon 
the temperature alone may be included in an 
additive term which, being independent of the 
strain, plays no part in the following treatment 
of isothermal processes. Thermal expansion does 
not constitute a strain in the sense of the present 
theory, and no special significance is attached to 
absolute zero. 

The third-order terms of the free energy have 
been given by Born and Misra® for a cubic 


+ (ns3+ 11) (137+ n31?)}. (12) 





Bravais lattice with central forces between 
particles at absolute zero. In this special case, 
Cauchy’s relation holds: ¢::=¢4, and by com- 
parison with (12) the following relations among 
third-order coefficients are found : 


Cyse=8Ci23;  Crsa=2Ci2z; Cise=4Cire. (13) 


Cauchy’s relation is rarely satisfied by real 
crystals, and it seems probable that the relations 
(13) are also rarely valid. 

When, however, (13) is assumed to hold, the 
third-order coefficients may be expressed in 
terms of the 5 lattice sums computed by Born 
and Misra, as follows: 


Cin = 8C(3), 
Cizs=48C(1,1,1), 
Cire = 24C (2,1). 


C11 =8€(2), 


C12 = 8C (1,1), (14) 


The expansion of ¢ in terms of the matrix jn 
might prove awkward in the general case, but 
for the particular strain (2) it is possible to 
express the j,, in terms of the Eulerian strains, 
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€rs. DiNce 
0a» Od, Ody 
; =——, and —=(1—a)-(5-—Bp,), 
Joa OX, OX, OX, + 
we find 


joa= (1 —a)~*(5pqg— 2B pq) =Spq—2epg. (15) 


Assuming that in this case, ¢ is expressible in 
terms Of. €pg, we may replace 0¢/Ajp_ by 
—}(8¢/d€p~). Also ¢ must have the same sym- 
metry in terms of the e,, as in terms of the 7,., 
and the second-order coefficients must be the 
same, since the expressions must coincide for 
infinitesimal strain, where ¢,, and 7,, become 
indistinguishable. Thus we obtain an expansion 
for pod in terms for ¢,, exactly like (12), except 
that we write ¢,, instead of 7,,, and introduce 
new third-order coefficients, Di,,; and so on in 
place of the corresponding C’s. | 

These expressions for ¢ must be identical, and 
hence there must be relations between the C’s 
and the D’s. One such relation is readily found 
by considering a hydrostatic strain, letting 
8,,=0, and putting both expressions in terms of 
a. We have first, in terms of n, 


pop =3n*CotnCetn'Cat:--; 
in terms of e, (16) 


pp =FC2+EDst+ Dit ---; 
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where 


C2=Cut+2¢12; C3=3Cu1+Cies+6Ciia; 


Ds=3Diu+Di23+ 6Dirz; 


and C,, D,--- are possible coefficients of higher 
order terms. On substituting a for « and 9 
according to (9), and identifying the coefficients 
of like powers of @ in the two expressions, we 
obtain : 


D;=C3+6C2; Dyg=Cy+6C34+18Co, etc. (17) 


Evidently the identity of the two formulations 
can be assured to any power of a by the addition 
of terms of higher order. Where C; has been 
calculated from lattice theory or otherwise, D; 
may be found by (17). Other relations may be 
found by taking other types of strain. Thus if 
we take 8120, we find 


Dyat 2D 166 = Crss t+ 2Cre6 + 16644+4011+-8cr2. (18) 


In deriving this identity, care must be taken to 
include second-power terms of 812; the approxi- 
mations (8) should not be used. 

These formulations of the free energy are not 
dependent upon the form of the potential 
between particles, except insofar as the numerical 
values of the coefficients, c’s, C’s, and D's, are 
concerned. 


5. EFFECT OF PRESSURE ON ELASTIC CONSTANTS 


The symmetry of a cubic crystal is preserved under hydrostatic compression, so that the relations 


between the additional stresses and the infinitesimal additional strains must have the same form as 
for a cubic crystal under infinitesimal strain at zero pressure, but the effective elastic ‘‘constants”’ 
will depend upon the pressure. These relations now follow ; the algebra, which is tedious, is omitted. 
By substituting in (10) and noting that the density is given by (5) and (9), we find, for example, 


Ti: = (14+2n)-4(14+-811 —B22 — Bs) (06/911), (19) 


neglecting higher powers of the A's. 
We now set 73,:= —P+T7))’ and write 7,;’ in terms of the 8's and a new set of elastic constants, 
C11’, C12’, and C44’, which depend upon P. Expanding (19) and identifying the coefficients, we find 


Ty’ =C11'B11 +612’ (B22 +83) ; (20) 
—P=9(1+2n)“*(Cut+2cizt+aCst ---); (21) 
Cx’ = (1420) { C11 +-9(3e11 + 2¢12 + 6Cir1 +4Crr2) +29°(15 Cr+ Cres $+ 14Cir2) + + + +}; 


(22) 
rz’ = (1+2m)-8 {C12 +-9(Ci2s+-4Cir2 — C11) +.9?( C123 + 2Cir2— 3Cinr) + «+ «} 
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Similarly, expanding T12= 712" = 2c44’B1o, 
Qua’ = (14+-20)—4{ 2e44t+ (4C44+ 2611 +412 + Cras +2C ies) 
+ 9?(6Citr+2Ci23+12Cii2+2Ciss+4C 66) + - + +}, (23) 
To the first power in 7 these become 
C1’ =Crrt(2e11+ 2612 +-6Ci11 +4C iz) 


C12’ = Cia t+ (Cizs+4Cirze— C11 — C12) ; (24) 
' Can’ =Caat (Cast O11 +212 + Crsa/2 + Cite) ; 
also 
C11’ +2€12' = C11 +2¢12+29C3 = C2+2nCs. (25) 


The three parameters ¢11’, Ci2’, and ca’ determine completely the response of the compressed 
crystal to the additional small strains. The constants for the bending, extension, or twisting of 
variously oriented bars or plates under pressure are found by substituting these parameters jn 
Voigt’s formulas in place of the ¢11, C12, and ¢44. The equations of small motion in the compressed 
medium are found from T,, in the usual way, the displacements being the u,’ defined in (1). For 
uniform pressure, the stress equations of equilibrium involve only the derivatives of T,,’. 

Parallel development in the Eulerian form leads to: 


T= —P+T i’ =(1—2¢)*!2(1 — 3811 —B22— B33) (0¢/ex1) ; (26) 
—P=e(1—26)5*[¢i1+2¢12.+eDst+ :- +]; (27) 
1! = (1—26)¥!?{e41 — (Serr +612 — 6D 111 —4D 12) — €*(21 Darr + 3D 123+ 26D 112) + + ++}; 

C12’ = (1 —2e)9/? {42 — €(€11 +4012 — Dias —4D 112) — (3D + 3D 123 + 14D 12) + «+ +}; 

2e44= (1 —2€)8/?{ 2e44 — €(4044 +2611 +4012 — Daa — 2D 1066) 


— €(6Di11+2D 123 +12D 112 —2D 144—4D 166) +--+}. (28) 


. To the first power in e, 
11 =C11 — €(10€11 +6612 —6D 111 —4D 112), 


C12’ = C12 — €(€11 + 9€12 — D123 —4D 112), | (29) 
Cael = Caq— €(7 C44 +011 +2612 —Di4a/2 —D 06) ; 
also 
C11 + 2€42" = C11 +2012 — €(12C2—2D3) = C,—€(12C,—2D3). (30) 


In this form it is clear that the c’’s will all increase with the pressure if the D’s are sufficiently small. 
Since for small compressions, a, 7 and ¢ are all equal to the first order, we obtain further identified 


among the C’s and D’s: 
Dy23t+4D 112 = C123 + 4Ci12 + 8612; 


3Diurt2Dire=3 Cin +2Ci12 +6611 +4612; (31) 


as well as the relations (17) and (18) already given. Since it appears that the C’s are usually negative, 
while the second-order coefficients must be positive, the D’s will be numerically smaller than the C's. 

The relations between pressure and volume, and between pressure and compressibility, are of 
special interest as they offer the possibility of rationalizing a large mass of experimental data. For 
hydrostatic strain, we may set all 8,,=0, or we may put B=811=822=8s3s, with the other §’s zero. 
The former method is more direct; the latter may serve as a check. With the 6’s all zero, we have 
from (21), with Vo/V=(1+2n)-, 


P=}3[(V/Vo)-*—(V/Vo)*][C2—Cs/2+(Cs/2)(V/Vo)*+--- J. (32) 
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Denoting the ordinary compressibility, —(1/Vo)(9V/dP)r by xo and its reciprocal, the “bulk 
modulus” or ‘‘incompressibility’’ by Ko, we find for the limiting case of zero pressure, 


(0K o/8P) poo = 1 — (2C3/3C,). (33) 

In terms of ¢, with Vo/V=(1—2e)!, 
P=$[(Vo/V)"8—(Vo/V)**][C2+D3/2 —(D3/2)(Vo/ V)**+ - - -]. (34) 
(0K o/OP) poo =5 —(2D3/3C:). (35) 


But since Ds=C3+6C2, (33) and (35) are identical. 

The same result is of course obtained if we put 8,,=85,,, where 8 is a small constant. In this case, 
we have —AP=Ti' =B(C11'+2¢12’), and AV/Vo=38(1—a)*=38(1—2¢)-9; Ko= —(VoAP)/AV, so 
that Ko=}(C2—€(15C2—2Ds3)+ oes ] = C2/3+P(5—2D3/3C2)+ -++, Thus in the limit, for P=0, 
we have again (35). The combination 3/(¢:’+2c.2’), gives the “instantaneous” compressibility, that 
is, the compressibility referred to the volume at the pressure P, rather than to the volume at zero 
pressure. 

We may remark that all of these expressions become identical with those of the isotropic medium® 
on setting ¢11.=A+2y, Cie=A, and Cyuy=p. 

It is also noteworthy that (33) and (35) are unchanged by the introduction of higher order terms 
in the free energy; terms contributed by powers of the strain higher than the third, disappear in 
passing to the limit of zero pressure. Equation (33) has also been found by Fiirth.‘ 





6. RELATIONSHIP TO OTHER THEORIES tion. Thus we have 


Since @ is the Helmholtz free energy, it Ky=— Vo(aP/aV) r= Vo(d%/dV?), 
must satisfy the thermodynamic relation, —P 
=(d¢/9V)r= —p*(d¢/dp)r, when the stress sys- and 
tem degenerates to a hydrostatic pressure P. In (8K o/8P) pao= — Vo(d*d/8 V*)o/(8%/AV)o; 
this case, the stress-strain equations (10) take 


the form or if we take a linear dimension such that 
V=r', we find 
0g or) 
tata Oh ataliaeeet  Rat (36) (0K o/dP)o = 2 —19(0*¢/dr*)o/3(0°/dr*) o. (37) 
711 €11 
But The partial derivatives with respect to the vol- 


=f(1+2 1+2 1 ial ume may be obtained from (12) or its equivalent 

o/po=C (1+ 2m1)(1 +222) (1 +230) J in terms of ¢ or 7 by using (5) and (9); we 
=[(1—2e11)(1—2¢22)(1—2es3) ]#. return to the relations already found. 

Hence If, however, the form of ¢ is known or as- 

8p/den= —o(1+291)-", sumed in terms of the interparticle distances and 

o/Onu ott +2) force constants, additional relations may be de- 


and 
Op/den1 = — p(1—2e1:)—". rived. With the power-law favored by Born, 
a where the potential between two isolated parti- 
° cles is of the form —ar-"+5r~-*, (37) leads at 
86/An11 = (86/8p)(8p/dn) 11 once to the well-known expression (dK»/dP)o 
and =(9+m-+n)/3; this is strictly valid only for 
O6/€11 = (06/8p) (Op/de11), T=0. Thus we find that D;/C;.=3—(m+n)/2 


and C;/C.=—3—(m-+n)/2. Conversely, if 


We might also have proceeded to find the (8Ko/8P)o is known, the sum (m+n) is deter- 
compressibility and its change with pressure by mined. This method has been applied by Fiirth’ 


successive differentiation of this general equa- 7 R. Firth, Proc. Roy. Soc. A183, 87 (1944-45). 


we verify at once that —P = —p*(0¢/dp)r. 











TABLE II. Relative values of third-order coefficients for 


power-law potential. 











(m+n) Ci/C2 D3s/C2 (8Ko/dP)o 
3 —4.5 +1.5 4.0 
4 —5.0 +1.0 4.33 
5 —5.5 +0.5 4.67 
6 — 6.0 0 5.00 
7 —6.5 —0.5 5.33 
8 —7.0 —1.0 5.67 
9 —7.5 —1.5 6.00 

10 — 8.0 —2.0 6.33 
11 —8.5 —2.5 6.67 
12 —9.0 —3.0 7.00 
14 — 10.0 —4.0 7.67 
16 —11.0 —5.0 8.33 
18 — 12.0 —6.0 9.00 








in an effort to show that consistent values of m 
and can be found from several independent ex- 
perimental quantities. 

As m and n are intrinsically positive, Ds will 
always be smaller than C3; values of the ratios 
of these coefficients to C, for various values of 
(m+n) are shown in Table II. Since, as it will 
appear later, the most common experimental 
values for (0Ko/0P)o are in the neighborhood of 
5, it may be justifiable for many purposes to 
neglect D;/C: where the neglect of C3;/C: would 
introduce serious error. It is thus advantageous 
to treat most of the questions involving large 
compressions in terms of the Eulerian strain; 
the same results may be obtained with the 
Lagrangian form, but at the price of retaining 
terms of higher order. 

Pressure-volume relations are furnished by 
Born and co-workers for a number of special 
cases, but the results are given in numerical 
form and comparison is not immediate, except 
for zero temperature. At T=0, however, it is 
readily found that the power-law potential leads 
to the general form, 


P=constant 
x c( Vo/V)8+™/3— (Vo/V) (3+n)/3)], (38) 


This equals (34) for m=2, n=4, in which case 
D;=0. For other values of m and n, these 
expressions will be equal only with the inclusion 
of higher order terms in (34). 

The higher order coefficients may be found 
for the special case of hydrostatic pressure by 
developing the power-law potential in terms of 
e and of ». Putting Vo/V=(r/r)*, we have 
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(ro/r) =(1—2e)*=(1+2n)-!. Then if 


non dl (>) (2) 
m\r n\r 
it follows that 
pod =constant+A {(m—m)e?/2! 
—[("—2)(n—4) —(m—2)(m—4)]e/3! 
+[(n—2)(n—4)(n—6) 
—(m—2)(m—4)(m—6)Je/4!— ++}; (39) 
and also that, alternatively, 
pod =constant+A {(m—m)n?/2! 
—[(m+2)(m+4) — (m+2)(m +4) ]n*/3! 
+[(n+2)(n+4)(n+6) 
—(m+2)(m+4)(m+6) |n*/4!—---}. (39a) 


Whenever m and n are both even integers, the 
development in ¢ becomes simply a polynomial; 
if either m or m is odd, (39) gives an infinite 
series which ultimately diverges for very high 
compressions (see Table I). It is hardly to be 
expected that the assumed power-law potential 
remains valid for such compressions, regardless 
of its suitability for small volume changes. The 
expression in 9 is always convergent, since |y| 
never exceeds 3, for compression; but also it is 
never a polynomial. The terminating property 
of (39) possibly accounts in part for its usefulness 
in dealing with reasonably large compressions. 
Except for ionic crystals, good approximations 
to the potential can probably always be obtained 
with even values of m and n (within the scope of 
usefulness of the power-law form of potential) 
and it seems unlikely that the larger of these 


TABLE III. c44’/c44 as function of volume. 








cad’ [cus 





3 —2e Vo/V (Birch) (Gow) 
m=6, n=12 
0 0 1.000 1.000 1.000 
0.05 0.0164 1.0247 1.125 1.193 
10 .0323 1.0488 1.247 1.405 
5 .0477 1.0724 1.363 1.635 
m=4,n=8 
0.05 0.0247 1.0375 1.189 1.207 
10 .0495 1.075 1.378 1.436 
5 .0818 1.125 1.625 1.711 
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need ever exceed 12. For n=12, the expansion 
terminates with the sixth power of ¢. Smaller 
values seem more suitable for most of the metals. 
The expansion in ¢ is given below for several 
pairs of values of m and n. 


For m=2, n=4; pod =constant+Ae. 

For m=2, n=6; pod=constant+A’e(1— je). 

For m=4, n=8; pod =constant+A”e(1 —2e+e). 

For m=1, n=8; a good approximation for most 
of the alkali halides,* po?=constant+A’”’é 
X(1-e+de+he+- ++). 


These expressions are found from the power- 
law potential and hence possess the same restric- 
tion to zero temperature, but it appears probable 
that the differential effect of temperature upon 
the coefficients of the various powers of « will not 
be great, at least for moderate temperatures, and 
hence that the values for such ratios as D;/C2, 
D,/C; and so on, which are immediately obtain- 
able from the above expansions for zero, will be 
approximately valid at higher temperatures. We 
notice that for m=2, n=4, we have D,; and all 
higher order D’s equal to zero, and the pressure 
is given by the relatively simple expression, 
P=constant X[(Vo/V)7*—(Vo/V)**]. Ds de- 
pends only upon (m+m) but the other terms 
depend upon the individual values of m and m. 
Thus for m=2; n=6, D, and higher coefficients 
are zero, while for m=1, n=7, D,/C.=9/8, 
D;/C:=3/8, and so on. 

In treatments which deal only with the case of 
hydrostatic pressure, development of the_po- 
tential in terms of strain components is of course 
unnecessary; any function of the volume may 
serve, subject to certain restrictions. Thus Slater 
expands the free energy in terms of (1—r/ro). 
This is equal to 1—(1+2)* which becomes 
approximately —y for small compressions. If, 
however, the individual elastic constants are of 
interest and not merely the compression and its 
various derivatives, it is necessary to expand the 
free energy in terms of the strain components, 

as has been done by Born and his associates, 
using the equivalent of 7,.. It is possible that in 
the general case, as well as in the special case of 
hydrostatic compression, the use of the é,. 


* J.C. Slater, Introduction to Chemical Physics (McGraw- 
Hill Book Company, Inc., New York, 1939), p. 393. 
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TABLE IV. Calculated effect of pressure on ¢4’ for 
several substances. 











1 dead 
cu OP 
Substance 7 +(¢11 +2c12) /cas (P in bars) 
Na 9.85 70.2 X 10-* 
Cu 12.60 2.97 
Ag 13.86 4.63 
Au 20.2 3.83 
Al 15.28 6.59 
WwW 13.02 1.41 
Fe 11.46 2.22 
MgO 10.0 2.15 
NaCl 12.74 17.2 
KBr 14.22 31.8 
CaF», 14.52 5.7 








instead of the 7,, might lead to more rapidly 
convergent results. 

A paper by Miss Gow’* gives the variation with 
pressure (and temperature) of all of the elastic 
constants for face-centered monatomic cubic 
lattices, using a potential of the form just con- 
sidered ; the lowest values of m and n for which 
calculations are presented are m=4, n=8. The 
results are given in numerical tables, in terms of 
a parameter £ related to the volume by (1+) 
=(V./V)‘*-™*, Comparison with the results of 
this paper can be carried out by computing the 
changes of the corresponding elastic constants 
for the same changes of volume. This has been 
done for c4, in Table III, where we have taken 
the approximate form of (29), neglecting the 
third-order coefficients Diy, and Dyge: 


Cag’ /Cos=1—€L7 + (C11 +2€12)/Caa]. (40) 


As shown in Table IV, the coefficient of € in (40) 
varies between 10 and 20 for various materials, 
most of the values falling near 14. In Table III, 
the figure 15.28, the value for aluminum, has 
been used. 

Miss Gow gives 


C4a/Po=(Vo/ V)[As/ut+B.T/6], 


where A;/u and B, are tabulated as functions of 
£; Po, @ are constants, T the absolute tempera- 
ture. This cy is equivalent to our cq’. For 
comparison with (29), we form the ratio for 
T=0, Cas’(€)/Caa(0) =(Vo/ V)[As(€)/As(0)]; this 
ratio is not greatly affected by including the 


®* Margaret M. Gow, Proc. Camb. Phil. Soc. 40, 151 
(1944). 
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TABLE V. Experimental values of change of compressi- 
bility with pressure: elements (pressure in bars).* 











Element °C 10% 10%  —-2b/a*=(8Ko/AP)o 
W 30 0.32 . 01 2 
75.32 a. --4 
Pt 30 ~—-.363 
Rh 30 ~—-.364 1.5 23 
Cr 30.525 9 7 
Ni 30.535 9 6 
75 S41 8 5 
Pd 30.534 9 6 
Co 30.546 8&8 5 
B 30 ~—-.558 - Ff 
Au 30.584 2.0 12 
75 577 8 5 
Cb 30.577 9 § 
Fe 24 «5941 83 4.7 (to 30,000) 
(Ebert)** 20 .603 1.37 7.5 (to 5,000) 
V 30.617 1.35 7.1 
Cu 30. ~—-.729 16 6.0 
.7182 9 4 (to 30,000) 
(Ebert)** 20 731 2.6 9.7 (to 5,000) 
n 30 ~—-.803 4.2 13 
Ag 30 1.002 3.7 74 
(Ebert)** 20 1.025 1.9 3.6 (to 5,000) 
Al 30 1.365 49 5.3 
1.328 3.7 4.2 (to 30,000) 
(Ebert)** 20 1.370 3.1 3.3 (to 5,000) 
(Birch) 23 1.378 7.2 7.6 (to 10,000) 
301.435 48 4.6 
Th 30 1.850 12 7 
Pb 30 2415 196 68 
2.340 10.2 3.7 (to 30,000) 
(Birch) 20 2.354 143 5.1 (to 10,000) 
(Ebert)** 20 2.448 145 4.9(to 5,000) 
Pr 30 3.451 12.2 2.0 
La 30 3.578 13.9 2.2 
Ca 30 5.805 65 3.9 see Table VIII 
Sr 30 8.278 103 3.0 
Ba 30 10.39 = 133 2.5 








* Unless otherwise noted, the measurements are by P. W. Bridgman, 


and the range is 12,000 kg/cm*. Several corrections have been applied 


to the values published in the original papers and in The Physics of 
High Pressure (Bell and Sons, London, Pamy me 1931), notably for 
6 for iron given in Proc. Am. Acad. Arts Sci. 74, 11-20 

(1940). See also Proc. Am. Acad. Arts Sci. 70, 285-317 (1935), and 
Rev. Mod. Phys. 18, 1 (1946). Corrected values, expressed as above 
in bars rather than in the ice are tabulated in Handbook 
ts) ent & Constants, ae 36, published by the Geological 
America, Section Values of a Po 6 for the new measure- 


mar 4 
1935). Those deny “pie h” are from F. Birch and R. 
ull. Geol. Soc. Am. 46, 1219 (1935). 
** Single crystal. 


temperature terms. Values of the ratio, as a 
function of volume, are given in Table III, for 
m=6, n=12, and for m=4, n=8. Again it is 
clear that better agreement would exist for 
smaller values of m and n. 

The magnitude of the neglected terms in (29), 
(Diss t+2D166)/2, relative to the quantity re- 
tained, (7¢44+¢11+2c¢i2), may be estimated by 
using Fiirth’s evaluation‘ of the third-order 
coefficients, and the identity (18). The C’s and 
the c’s are given in terms of the quantity €(2) 
for certain combinations of m and n. In this way, 
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it may be estimated that for the face- -centered 
cubic lattice, with m=4, n=8, the neglected 
terms would increase the pressure coefficient by 
about 14 percent; for m=4, n=6, the neheeail 
terms would reduce the coefficient by about 4 
percent. The few experimental determinationsw 
of the effect of pressure on rigidity support the 
calculated values of Table IV: the measured 
values of the relative change of rigidity with 
pressure are 2.7-10-® for copper, 7.6-10-® for 
aluminum, and 2.36-10-* for steel (P in bars), 
However, these are for polycrystalline materials, 
and are not strictly comparable with the calcy- 
lated values for a particular direction in the 
crystal, although it seems probable that the 
effect will not differ greatly for different orienta- 
tions. 

The correction to Born’s theory arising from 
the stress-strain relation (10) for finite strain 
will apply only to the expressions for the elastic 
constants; the pressure-volume relation escapes 
this correction, since it is based on the exact 
thermodynamic relation for the pressure as 
volume derivative of the free energy. The relation 
(10) should be regarded as the exact extension 
of this thermodynamic relation to the general 
finite strain. 


7. COMPARISON WITH MEASUREMENTS 
OF COMPRESSION 


The compression of a large number of ma- 
terials crystallizing in the cubic system has 
been measured over pressure ranges sufficiently 
great to establish a definite departure from 
linearity. Most of the measurements have 
been made by Bridgman, who has extended the 
working range from about 12,000 kg/cm? in his 
early work to 100,000 kg/cm? in recent publica- 
tions. The results for the lower ranges are usually 
given in the form, —AV/Vo=aP—bP"’, where 
AV is the change of volume from the initial 
volume Vo at zero pressure, to the pressure P; 
a and b depend upon the temperature. Thus, at 
constant temperature, (0Ko/dP))=2b/a?. Table 
V includes the measurements for all of the 
elements of cubic symmetry except for those of 
high compressibility (Table VII and VIII), and 
several (Mo, Ta, Ce) which show anomalous 


” P. W. Bridgman, Proc. Am. Acad. 63, 401 (1929); 
F. Birch, J. App. Phys. 8, 129 (1937). 
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negative values of b. Table VI includes measure- 
ments for a number of cubic compounds. Meas- 
urements made before 1940 have been corrected 
to the new value determined for iron at that 
time; this correction amounts to a reduction of 
the former values for b by about 1.5-10-" and is 
thus of importance for the materials of low 
compressibility. 

Most of the values of 2b/a? in both Table V 
and Table VI lie between 4 and 7; if the distri- 
bution is considered to represent a normal 
probability curve, the most probable value for 
the elements is about 4.5 with a half-width of 
about 0.5, an amount not inconsistent with the 
experimental uncertainty when reasonable allow- 
ance is made for such defects of material as 
initial strains and porosity. Occasional deviations 
seem to be of this character, as for gold, where 
an erratic value occurs for 30°C, a normal one 
for 75°C. The values for rhodium and for manga- 
nese are almost certainly erratic. The most 
probable value for the compounds is between 5 
and 6. The only significant deviation. is for the 
sulfides and perhaps certain fluorides, where 
2b/a? is about 12. It is noteworthy that the latest 
and most precise measurements, to 30,000 kg/ 
cm’, show even lower values for 2b/a? than the 
earlier ones for such relatively reliable materials 
as Al, Pb, and Cu. Figures lower than 5 mean 
positive values of D;, and correspond to very 
small exponents in the power-law potential 
(Table 11). Low values occur both at the be- 
ginning and end of Table V, which is arranged 
in the order of increasing compressibility. Not 
too much weight should be attached to the low 
value for tungsten, where the experimental error 
must be relatively large ; but the same reasoning 
would suggest that the values for Ca and Sr 
should be very reliable. Other data for Ca in Table 
VIII demonstrate a more “normal” variation. 

Slater" has concluded that 5 is probably the 
best general value of 2b/a? for the metals; this 
leads to a figure of 1.833 for Gruneisen’s con- 
stant, as compared with about 2.4 preferred by 
Fiirth. An outstanding feature of Tables V and 
VI is the evidence that (0Ko/dP)») has approxi- 
mately the same value for such diverse materials 
as metals and various types of ionic compounds. 


uJ. C. Slater, Phys. Rev. 57, 744 (1940). 
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_TAaBLE VI. Experimental values of change of compressi- 
bility with pressure: compounds. Measurements at 30°C 
(P in bars).* 














Compound 10*a 10!2% 2b/a*® 
TiN 0.334 0.87 16 
TiC 477 8 7 
Pyrope (garnet) .545 9 6 
Fe;O, .547 .82 5.5 
MgO .598 1.0 5.5 
Andradite (garnet) .673 .86 3.8 
Fe.S .68 .87 3.8 
CoAsS 767 1.88 6.4 
CdF. 1.120 7.5 12 
ZnS 1.303 1.28 1.5 
1.258 5.4 6.8 (to 50,000) 
CaF», 1.226 6.49 8.7 
LiF (S) 1.53 7.6 6.5 
1.52 5.5 4.8 
SrF;, 1.607 9.4 7.3 
PbS 1.869 7.43 4.3 
Cu,0 1.943 19.3 10.1 
BaF. 1.964 14 7.3 
NaF 2.11 17 7.7 
CaS 2.321 39 14.5 
SrS 2.426 38.8 13.2 
CuCl 2.508 13.3 4.2 
Cul 2.802 24.8 6.3 
BaS 2.946 52.8 12.1 
Ag-S 3.265 49.5 9.4 
KF (S) 3.30 32 5.9 
LiCl (S) 3.40 32.3 5.6 
NaCl 4.260 51 5.6 see Table IX 
LiBr (S) 4.30 51 5.6 
NaBrO; 4.401 76 7.8 
NaClo; 5.033 95 7.5 
NaBr (S) 5.07 63 4.9 see Table IX 
KCI (S) 5.62 75 4.8 
CsCl 5.94 99.5 5.7 see Table IX 
Lil 6.00 110 6.1 
KBr (S) 6.70 105 4.7 
CsBr 7.05 149 6.0 see Table IX 
RbBr (S) 7.93 137 4.4 
KI (S) 8.53 155 4.3 
CsI 8.565 209 5.7 see Table IX 
RbI (S) 9.56 204 4.5 








* All measurements by P. W. Bridgman (see references for Table V) 
cnet gg marked “S,"’ which are by J. C. Slater, Phys. Rev. 23, 
488 (1924). 


This is understandable if the ratio (—D;/C2), 
which is the part of (@K»/dP)) which depends 
upon the individual peculiarities of the inter- 
particle forces, is generally no greater than 
unity, or in other terms, (m+) no greater than 
8 or 9. 

The same conclusion may be reached by 
examining the striking body of data for com- 
pressions up to 100,000 kg/cm* given in two 
papers by Bridgman, as well as some earlier 
results to 45,000. The data are given as —AV/ Vo 
for various pressures. In a number of instances 
the density was more than doubled, and in many 
cases, raised by 25 percent or more. Taking the 
pressure-volume relation in the form (34), we 
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ression of the alkali metals to 45,000 


TaBLe VII. Com 
oe and values of F/P. 
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TaBLeE VIII. Compressions to 100, 000 k 2 
of F/P.* g/cm? and Values 








P 
hy Li Na K Rb Cs 
em? V/Vo F/P V/Vo F/P V/Vo F/P’V/Vo F/P V/Vo F/P 


ee 


P 
in 108 Li Na K Rb 
kg/ Ca 
em? V/Vo F/P V/Ve F/P V/Vo F/P V/Ve F/P V/V, Rip 
Pe 





5 0.957 641 0.929 113 210 0.826 374 0. 
10 .926 599 .883 107 .817 202 .767 301 .729 397 
15 .899 584 .852 98 .770 196 .721 280 .674 387 
= as 581 .818 100 .732 194 .684 271 .628 395 
30 =.835 S77 = .767 100) =.671 197 629 261 .568 392 
35 .816 581 .746 100 .647 199 .607 259 .542 400 
40 .798 591 .727 101 627 199 .587 259 .519 409 
45 .782 S97 .710 101 .604 206 .569 260 .499 417 








Bridgman, Phys. Rev. 57, 237 (1940); Proc. Am. Acad. Arts 
Sci. 72. woo7 (1938). 


see that if it is permissible to neglect D3/C2, 
we obtain the simpler form, 


P= FC,/2, 
=[(Vo/V)"*—(Vo/V)**). 


The ratio F/P is in this case a constant, 2/C2 
=2x0/3. The following test has been applied to 
the measurements of Tables VII, VIII, and IX: 
F has been formed from the measured values of 
Vo/ V and divided by the corresponding measured 
pressures. The test is extremely sensitive since 
F varies rapidly with Vo/V. All, or nearly all, 
of the measurements of inorganic cubic materials 
have been included except those which show 
polymorphic changes in the pressure range, a 
fairly large fraction of the whole. 

Let us first examine the alkali metals in the 
45,000 kg/cm? range, shown in Table VII. F/P 
is constant within the precision of the measure- 
ments for all of these metals except Rb, where 
there is an abnormally high initial compressi- 
bility ; the first value for Li also deviates widely. 
For the others, the initial compressibility can be 
obtained with an error of less than about 10 
percent from any single volume determination. 
There is no particular trend which would suggest 
that improved agreement would result from the 
next approximation, with a non-zero D3. 

The measurements on the alkali metals to 
100,000 kg/cm? (Table VIII), do show a definite 
trend, with F/P increasing beyond 50,000, by 
roughly 50 percent for Li and K, nearly 100 
percent for Rb, and about 20 percent for Na. 
The values above 50,000 can be fitted by in- 
cluding a term in D3, as in (34), but D; turns 
out to be large and positive and far beyond any 


with 
(41) 


10 0.926 598 0.767 301 
20. 875 S81 .818 100 .732 194 .684 271 (997 451 











30 835 S75 .767 100 .671 197 629 261 861 449 
40 .798 S91 .727 101 627 199 .587 259 339 443 
50.763 621 .690 104 581 216 .538 288 ‘805 44 
60 .728 668 .657 108 .543 232 499 313 789 454 
70 «=€=.695 720 .628 112 .508 252 461 351 (Transi- 
80 .664 776 .603 116 478 272 425 401 tion) 
90 .634 842 581 120 .451 293 .392 459 
100 4.606 913 .560 124 .427 316 .362 527 
| 
*P. W. Bridgman, Proc. Am. Acad. Arts Sci. 74, 425 (1942); or 


Phys. Rev. 60, 351 (1941). 


value compatible with the measurements below 
50,000. It seems fairly clear that the determina- 
tions for these two ranges are not mutually 
consistent. Bridgman has discussed the possible 
sources of fairly large errors in the higher range, 
In view of later determinations to 100,000, we 
may conclude that the abnormal behavior of the 
alkali metals above 50,000 is somehow connected 
with the method of correcting the data, which 
was improved in the later work. For example, 
as shown in Table IX, if the data for indium 
(a tetragonal crystal) given in the two papers 
(1942 and 1945) are reduced in the same way, 
we find that F/P increases in the former, 
though much less than for the alkali metals, 
whereas it decreases in the latter. 

A decrease or approximate constancy of F/P 
is characteristic of the 1945 determinations, 
shown in Table IX. The decrease is readily 
accounted for by restoring D3, which may be 
evaluated by plotting P/F versus (V)/V)}; 
according to (34) this should be a straight line. 
In a fairly typical case, we find D;/C,= —0.8 for 
CsI. This corresponds to a value of 5.53 for 
2b/a*, whereas from Table VI, we find 5.7; and 
from Table II, (m+n) between 7 and 8. 

With the retention of the constant Ds, the 
measurements of Table IX can all be reproduced 
with an error not exceeding about 3 percent on 
the pressure. This is probably less than the 
experimental error of pressure determination for 
this range, and there is some indication from 
the data that the experimental pressures may 
be perhaps 10 percent lower than supposed. We 
reach this conclusion by observing that the ratio 
(F/P) 10,000 to (F/P) 100,00 is close to 1.1 for 9 of 
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the 14 materials of Table IX; it is nearly con- 
stant for AgCl and CsBr and equal to about 1.2 
for CsCl and TINOs. There appears to be no 
correlation with the amount of compression. 

The approximate form (41) thus gives the 
pressure corresponding to the extreme compres- 
sion in terms of the measured compression at 
10,000 with an error not exceeding 20 percent, 
and generally less than 10 percent for the entries 
of Table 1X. This is probably not far from the 
limits of experimental precision for this range of 
pressure. 

Bardeen” has calculated the zero-point energy 
as function of density for Li and Na by approxi- 
mate quantum-mechanical methods, obtaining 
the energy in the form, ¢=Ay*+By’— Cy, with 
y=(Vo/V)*. A, B, and C are constants which he 
calculates and compares with the experimental 
values. This leads to a formula for the pressure 
which is approximately, 


P=constant[(Vo/V)*/*—(Vo/V)**], (42) 


equivalent to the form for a power-law potential 
with m=1, n=2. Comparison with experiment 
requires a somewhat uncertain reduction to 
absolute zero, and the numerical agreement for 
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the constants is not very close. If, however, the 
constant of (42) is determined from the measured 
initial compressibility, the agreement with the 
measured compression of the alkali metals to 
45,000 is very nearly as close as for (41). Fiirth’ 
has applied the power-law potential to some of 
the measurements to 100,000, with values of m 
and m larger than are implied by (41); the 
agreement with the experimental values is less 
good. It appears that larger values of the expo- 
nents in the force law are required to satisfy the 
experimental determinations of temperature ef- 
fects than to fit fhe measurements of pressure 
effects. It may be questioned whether this is 
mainly the consequence of an inadequacy of this 
form of potential, of some approximation in 
Born’s treatment of the entropy of the lattice, 
or of more fundamental causes, as suggested by 
Bridgman in a discussion of the implications of 
Schottky’s theorem." 


8. CONCLUSION 


Besides the application of the theory of finite 
strain to find the effect of pressure upon the 
elasticity of cubic crystals, we have attempted a 
clarification of certain questions regarding the 


TABLE IX. Compressions of cubic compounds and values of F/P. Compression of indium.* 








P 
in 10° F/P V/Vo Dal 


V/Vo F/P V/Ve F/P V/Vo F/P V/Vo F/P V/ 
TIC! AgCl 


Ve F/P 
NaBr CsBr NH.«Br 





10 0.962 279 0.952 362 0.962 279 0.979 145 0.957 321 0.947 405 0.951 370 


20 .932 270 914 359 929 285 
907 264 .882 359 -901 285 
885 260 856 353 877 284 
865 258 834 348 -856 282 
848 254 816 338 838 279 
832 253 801 329 823 274 
817 253 -788 319 -809 270 
803 252 777 = 309 798 262 
790 252 767 = 301 -787 258 


SSSsssss 


— 


960 149 922 319 905 406 912 369 
942 150 893 315 870 409 880 367 





TINO; Nal CsI 


926 149 868 313 840 411 .853 364 
910 152 847 309 814 416 .830 361 
896 152 .829 303 792 413 811 354 
883 153 814 296 773 «411 .795 346 
871 152 801 288 -757 405 .781 338 
-860 151 -789 282 742 402 .768 332 

778 = =6277 .728 400 -756 327 

Indium 
Tu 1940 1942 1945 





10 0.963 271 0.944 429 0.935 513 0.950 379 0.978 155 0.977 166 0.9760 170 


20 932 270 902 424 -887 508 
30 905 271 868 417 849 505 
40 882 266 840 41i 818 501 
50 863 264 816 406 792 496 
60 848 255 795 404 770 490 
70 835 246 777 = 398 -751 484 
80 824 238 761 393 734 479 
90 815 228 -747 = 388 719 473 
100 .806 222 734 383 706 467 


910 379 958 156 955 168 9558 164 
877 378 940 155 936 167 -9381 162 
848 381 925 152 919 167 
823 384 913 146 903 166 9068 157 


9219 159 


802 382 888 167 8928 157 
7184 378 874 168 8798 158 
768 373 .860 170 .8676 156 
755 366 847 172 -8561 158 
744 356 835 173 8451 157 








* P. W. Bridgman, Proc. Am. Acad. Arts Sci. 76, 1 (1945), (Halogen compounds and indium); 74, 425 (1942), (Indium); 74, 21 (1940), (Indium). 


® J. Bardeen, J. Chem. Phys. 6, 372 (1938). 


1% P, W. Bridgman, Rev. Mod. Phys. 7, 1 (1935). 
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scope of this theory. It is shown that there is no 
ambiguity: the alternative methods developed 
by Murnaghan lead to identical results for com- 
parable approximations. The choice of method is 
a matter of convenience; for studies of large 
compressions, the Eulerian form furnishes sim- 
pler, more compact relations. Lattice theory does 
not conflict with this method; in cases where it 
can be sufficiently developed, it furnishes numer- 
ical values for the elastic constants in terms of 
atomic parameters, but many relations of a 
general nature are obtainable without a detailed 
specification of lattice structure and forces. The 
present method is valid for elastic strain at any 
constant temperature, and hence avoids the 
complexities which arise in lattice theory when 
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temperatures other than zero are considered, It 
is shown that a relatively simple approximation 
accounts for practically all of the measurements 
at high pressure. The fact that this first approxi- 
mation is also a good one, is shown to imply 
values for the exponents in a power-law form of 
potential considerably lower than have com. 
monly been adopted. 
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APPENDIX 


Derivation of Third-order Terms of the Free Energy for Cubic Crystals 


The general methods of taking into account the symmetry of a crystal for various properties are described in Voigt," 
Wooster! and other texts on crystal physics, and the application to the second-order elastic coefficients is given in con- 
siderable detail in Love.'* The symmetry of all classes of cubic crystals is such as to reduce the 21 possible second-order 
coefficients to 3. There are 56 possible distinct third-order coefficients, and this is the actual number for a triclinic crystal. 


.Using the customary short notation, where we write 1 for the pair of indices 11, 2 for 22, 3 for 33, 4 for 23 and 32,5 


for 31 and 13, and 6 for 12 and 21, the 56 third-order terms are as follows where evidently Cygr = Cprg= Cgpr and so on: 


Crrrm?® + Crran?net+ Crrami'ns + Cren?ns + Crsntns + Crenrne 
+ Cronin? + Cresninans + Crzeninene + Crssninens + Crises 
+Cissnin?? + Ciseninsnst Cissninans t+ Cisennane 
+Crsannd + Crsninenst Crusoninens 


+Cissmnse + Cisennsns 
+ Cisenné 


+Ca2an2* +Cr2anans +Ca2ena*ng +Caasna’ns + Co2ens’ns 
+ Coassnn?? + Casenenanet Cassnanans t+ Casenenane 
+Crumme? + Casnmenst Cums 


+Carssnms? + Casennsns 
+ Cosenne? 


+Casms® +Casena’ng +Cassns’ns + Casenene 
+Csaansnd + Coasnsnanst Coronsnane 


+Cassnns? + Casensnens 
+ Crsensne? 


+Cume +Cusnens +Cumens 


+Casnans’ + Casonenens 
+ Cusnané 

+Cossn® + Cosens?ne 
+ Cosens’. 


The crystallographers have amused themselves by inventing a variety of ways of designating the 32 crystal classes 
and of describing the elements of symmetry. Table X shows in parallel form a few of the proposals with regard to the 
5 classes of cubic or isometric symmetry. Only two of these classes are characterized by the tetragonal symmetry about 
the three cubic axes which one might expect to be a common property ; the other three classes possess digonal symmetry 


4 W. Voigt, Lehrbuch der Kristallphysik (B. G. Teubner, Leipzig, Germany, 1928). 
W. A. Wooster, A Textbook on Crystal Physics (The Cambridge University Press, New York, New York, 1938). 


New York, 1937). 


16 A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity (The Cambridge University Press, New York, 
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TABLE X. Designations and elements of symmetry of the isometric classes. 














———— 
Name or symbol of class 
' Hermann- Elements of 

Dana (Groth-Rogers) Miers Lewis Mauguin Schoenflies symmetry 
Hexoctahedral (galena) Ditesseral central Cubic II <3 2 ow C, 344, 44s, 6A; 
Gyroidal (cuprite) Tesseral holoaxial Cubic I 4 2 O 3A4, 4As, 6A2 
Hextetrahedral Ditesseral polar Cubic V 4q m T® 4As, 3A, 
Diploidal (pyrite) Tesseral central Cubic IV 2 3 T® C, 4As, 3A2 
Tetartoidal Tesseral polar Cubic III as T 4A3, 3A:2 














Note: In the last column, C denotes a center of symmetry, A« an axis of fourfold symmetry, As an axis of threefold symmetry, A: an axis of 
twofold symmetry. The numbers before the A’s denote the number of such axes. This notation is used in Dana and Hurlbut’s Manual of Miner- 
alogy (John Wiley and Sons, Inc., New York, 1941) and is by far the most intelligible and useful to the reader unversed in cry phy. The 
Hermann-Mauguin symbols are to be found in the first volume of International Tables for the Determination of Crystal Structures ( traeger, 
Berlin, Germany, 1935), The other notation has been extracted from Love. 


about these axes. All classes possess trigonal symmetry with respect to the 4 cube diagonals. This is expressed analytically 
as an invariance of properties under cyclic permutation of indices. Thus if the three original axes of coordinates are 
X,X2X3, coinciding with a set of cubic axes, then rotation of 120° about a trigonal axis defines a new set of 
axes, X;’ X2' X;’, say, such that X3’ coincides with X,, X3’ with X2 and X,’ with Xs, for example. This corresponds to 
the scheme: 
| XxX,’ xX,’ xX,’ 

Xi | 0 1 0 

X:| 0 O 1 

X3 1 0 0. 


The strain components referred to the new axes, 9,.’, will be related to the strain components referred to the original 
axes by the general transformation equation for second-order tensors, rs’ = apratgspq, P, 7, 7, S=1, 2, 3, where ay; is the 
cosine of the angle between X; and Xj’, and so on. For this transformation, this becomes simply 





m'=n3 ne’ =n6 
a’ =m ns =n 
ns’ =: ne =n. 


TaBLE XI. Effect of the elements of symmetry upon the Now the free energy must possess the same form with re- 
third-order elastic constants. spect to the new axes as to the old; thus we must have 

















Trigonal Diagonal Tetragonal Crm? + Crsmnet + ++ = Cum’ + Crm’ +>, 
aT — Tae ay a the C’s being necessarily the same on both sides. Substitu- 
Catamt C C tion for the 9,’ of their values in terms of 9, as given 
aah om _ os vt above, leads to equalities among the C’s. By this process, 
Cus = C332 = Cire Cus =Ci2 the 56 constants are reduced to 20, as shown in Table XI. 
114 = Coon = Core =0 In addition to the trigonal symmetry, all 5 isometric 
Ge Cn Cone = classes possess digonal symmetry about the 3 cubic axes. 
nos ae Care =0 The corresponding transformation scheme for rotation of 
Cus= Coea= Coss =0 180° about the X; axis is: 
Cus= Coes=C =0 ’ ’ ’ 
Cin= a “ =0 __| Xi! Xa! Xs 
Coa Cram Cone = Xi} -1 0 0 
126 = C135 = Cox ™ xX 0-1 O 
Cisa= Cree = Coss Crs Cras X, 0 0 1 
Cosa = Cie = Cass 166 166 ? , 
3144 = Coes = Ciss 156 = Cres so that 
qe= ao se — mn’ =m n= — 
146 = C356 = Caas = ! = ‘=-_— 
A wel =0 m2 = ™ = ms 
Ci2s Ciss Criss "3 =": "6 n6- 
456 Cass Cass Application of this transformation, taking into account 
Non-zero distinct the equalities already found, gives, for example, Ci:e:*n« 
— P ‘ = —Cien*n, and hence Cy14=0. By this process, the 20 


constants are reduced to a maximum of 8. 
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Finally, for two of the classes, the cubic axes are axes of tetragonal symmetry. For a rotation of 90°, we have the 
typical scheme 
| Xi’ Xs’ X,’ 
X,;| 0 -1 O 
X2 1 0 O 
X;| 0 S 8. 





leading to 
m’ =n2 n= —5 
ne’ =m m= 
ns’ =n3 ne =— He 


We find two additional equalities, C112=Ci1s, and Ciss=Ciss; the number of constants is thus reduced to 6. This is the 
minimum obtainable from considerations of symmetry alone, although further reduction may result from special assump. 
tions about lattice structure and forces, as Born and Misra have shown. 
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Radiations from As’** 


Attan C. G. MitcHet, Ep. T. JuRNEY, AND MARGARET RAMSEY 
Indiana University, Bloomington, Indiana 
May 1, 1947 


PRELIMINARY investigation of the radiations 

from As” has been made with the help of the coin- 
cidence counting apparatus available in this laboratory. 
Only meager information exists in the literature about this 
element. It is shown to be a positron emitter of approxi- 
mately 26 hours half-life. It was of interest to this labora- 
tory because it disintegrates to Ge”. The authors have 
recently studied the energy levels of Ge” when it is formed 
by beta-ray emission from Ga”, and it is to be expected 
that a study of the radiations of As™ would shed further 
light on the levels of Ge”. 

As” was prepared by bombarding gallium with 23 Mev 
alpha-particles from the cyclotron. A chemical separation 
was made, using arsenic carrier, and the active material 
showed a 26-hour period from As” and a small amount 
(<2 percent) of an activity with a period of several days, 
probably from the 16-day As” formed from the other iso- 
tope of gallium. 

The positron end point was obtained by measuring the 
range of the positrons in aluminum and applying the 
Feather rule. The value obtained was 2.78+0.10 Mev. 
Absorption curves were taken at various times after the 
end of the irradiation to make sure that this end point was 
that of As” and not caused by the small amount of As”*. 

The energy of the most energetic gamma-ray was ob- 
tained by placing an aluminum radiator in front of the 
source and measuring the range of the Compton electrons 
passing through two counters arranged in a coincidence 
circuit.1 The range of the Compton electrons was 1.10 
grams/cm* which, if one uses the curve of Curran, Dee, 
and Petrzilka,? corresponds to a gamma-ray energy of 
2.4 Mev. 

Positron-gamma and gamma-gamma coincidence ex- 
periments were also made. In the first instance the source 
was placed between two counters. With enough aluminum 
between source and counters to stop all positrons, the 
total non-particle coincidence rate was determined. This 

rate arises from annihilation radiation, true gamma-gamma 
coincidences and gamma-annihilation radiation coinci- 
dences. In the present experiments it served merely as a 
base for measuring particle-gamma coincidences. With the 
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source between the two counters and enough aluminum 
between the source and one counter to stop all positrons, 
positron-gamma coincidences were measured as a function 
of the thickness of aluminum between the source and the 
other counter. The results show that positron emission 
probably does not lead to the ground state of Ge” and that 
there is more than one group of positrons. 

In another experiment, lead-lined counters were used to 
measure gamma-rays and gamma-gamma _ coincidences. 
With the counters placed in line with the source between, 
the number of coincidences observed from the source was 
the sum of annihilation-gamma, gamma-gamma, and anni- 
hilation coincidences. When the line between the source 
and one of the counters was perpendicular to the line be- 
tween the source and the other counter, the coincidences 
due to annihilation radiation were greatly reduced. With 
the counters in line with the source, the ratio of coinci- 
dences per recorded gamma-ray was 3.89 X 10-*; when they 
were out of line, the ratio was 1.06 107°. 

The fact that a 2.4 Mev gamma-ray is present in the 
radiations from Ga™ as well as in those from As™ shows 
that possibly one level of Ge™ is excited by either disin- 
tegration. 

* This research was supported by a grant from the Office of Naval 


Research. 
sae tioge Mitchell, E. T. Jurney, and M. Ramsey, Phys. Rev. 71, 
2S. C. Curran, Dee, and Petrzilka, Proc. Roy. Soc. 169, 269 (1938). 





Erratum: The Analysis of the Vibration- 
Rotation Band w; for C!*O,’* and C"O,'* 


[Phys. Rev. 68, 173 (1945)] 

Atvin H. NIELSEN* AND Y. T. Yao** 
Mendenhall Laboratory of Physics, The Ohio State University, 
Columbus, Ohio 
May 5, 1947 


N re-examination of the data from which the band 
analysis for w; of C¥%O,"* and C'0O,"* was made, it 
was discovered that a mistake had been made in stating 
that the frequencies of the rotation lines had been reduced 
to vacuum. Through an oversight the lines given in Tables 
I and II are uncorrected. The constants have been recal- 
culated in order to determine how they would be affected 
by this error. The only constants given in Table III 
which are affected are wo for both molecules. For C"O,'* 
the corrected value of wo is 2349.33 cm which is now in 
excellent agreement with the value quoted by Herzberg! 
in his compilation of the CO; data. The corrected value of 
wo for C40,'* is 2283.56 cm™. Booo, Boor, and a3 are not 
affected because the slopes of the straight lines in Figs. 3-6 
are not altered measurably. 
To obtain the corrected frequencies the following table 
of Av may be applied: 


»(cm~) Av(cm™) 
2050 —0.532 
2100 —0.545 
2200 —0.571 
2250 —0 
2300 —0.597 
2350 —0.610 


* The University of Tennessee, Knoxville, Tennessee. 

** The National University of P Kunming, China. | 

1G. Herzberg, Infrared and Raman a of Polyatomic Molecules 
(D. Van Nostrand Company, Inc., New York, 1945) p. 274. 





Preliminary Determination of the Neutron 
Absorption Cross Section of 
Long-Lived I'** 


Seymour KATCOFF 
University of California, Los Alamos Scientific Laboratory, 
Santa Fe, New Mexico 
April 30, 1947 


HE following chain of mass number 129 is known to 
occur in fission: 


’ 4.2-hr. Sb->72-min. Te—long-lived I—stable Xe 


t 
32-day Te. 


An attempt was made by J. Seiler' to detect by radio- 
chemical means the long-lived I'*. He isolated the iodine 
from a uranium slug which had undergone irradiation in 
the Clinton pile for 123 days. The sample decayed essen- 
tially to background, with a half-life of 8.0 days which is 
characteristic of I'**. This indicated that the half-life of 
I was extremely long, or that its radiations were ex- 
tremely soft, or both. 

In the present experiment this same iodine sample was 
irradiated simultaneously with two normal iodine samples 
for 12.8 hours in the Los Alamos homogeneous pile, just 
outside the BeO reflector. Normal iodine, which is known 
to consist only of stable I’, should give on slow neutron 
irradiation only 25-min. I’. However, the sample which 
also contained some I'** yielded several thousand counts 
per minute of 12.6-hr. I'® after irradiation. All three sam- 
ples were purified radio-chemically, mounted as Agl, 
and counted with a mica-window Geiger counter. Absorp- 
tion of the radiations of the 12.5-hr. component in alumi- 
num indicated two beta-rays with maximum energies of 
0.54 and 1.02 Mev. This compares favorably with the 
values for I'® reported in the literature of 0.61 and 1.03 
Mev. 

The neutron absorption cross section of I* could be 
calculated from the data of this experiment if the number 
of I®* atoms in the samples were known. Unfortunately, 
this can be calculated only roughly because the isolation 
of the iodine from the uranium by Seiler was not done 
quantitatively. Two separate methods of making this 
estimate were used. The first is based on the fission cross 
section of uranium, the fission yield of the chain of mass 
129 (0.7 percent), the weight of uranium, the time of 
irradiation, the approximate neutron flux in which the 
uranium slug was irradiated, and the assumptions that 
about 30 percent of the iodine in the slug was extracted and 
that the chemical yield was also about 30 percent. This 
leads to a value of 7.610" for the number of I atoms 
and to a value of 2.5 barns (10-4 cm?) for the cross section. 
The other method of estimation is based on an approxi- 
mate measurement of the 8.0-day I'" activity which was 
present at the end of the 123-day irradiation in the iodine 
sample. This assumes that the bombardment of the 
uranium was steady, especially near the end. The fission 
yield of I'** was taken as 2.5 percent. This leads to a value 
of 1.910" I° atoms and to a value of 10 barns for the 
cross section. This is thought to be more reliable than the 
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first value, so that a weighted average of the two gives a 
final value of about 8 barns, which is very probably correct 
to within a factor of ten. 

The neutron flux passing through the iodine samples 
during their irradiation at the Los Alamos pile was calcy. 
lated from the activity of the 25-min. I and the cross 
section of stable I” (6.25 barns). In the calculations, gor. 
rections were made to 100 percent counter geometry, to 
100 percent chemical yield, to zero absorber, and for decay. 

A half-life of 10° years, or longer, for I is consistent 
with the results of Seiler’s experiment. Such a long period 
would make it possible for I'** to exist in nature. In fact, 
the two normal iodine samples used in this experiment did 
show a small amount of residual activity after the decay 
of the 25-min. I (see Fig. 1). Some of this may be due to 
13-day I"* formed by an (nm, 2m) reaction, and to impuri- 
ties. However, there is some indication of a 12.9-hr. 
period which may be due to 12.6-hr. I’. This sets an upper 
limit of 3 parts per million for the amount of I"** that might 
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Fic. 1. Decay of 12.6-hr. I! produced by neutron 
bombardment of long-lived I%*. 


be present in normal iodine if a cross section of 0.8 barn is 
taken. The limit is only 0.3 part per million if the most 
probable cross section of 8 barns is used. Nier reported? 
an upper limit from mass spectrographic data of 25 parts 
per million. This problem can be investigated further by 
bombarding larger samples of normal iodine in a higher 
thermal neutron flux and then subjecting them to more 
extensive purification. 


* This document is based on work performed at Los Alamos Scien- 
tific Laboratory of the University of California under Contract No. 
ete litt Mamie ees A aS a 

ap n Report , May 29, . 
It will also appear in Division V of the Manhattan Project Technical 
Series as part of the contribution of the Los Alamos Laboratory. 


1 Manhattan aig Progress Ri CN-1998, p. 2. 
2A. O. C. Nier, Phys. Rev. 52, 937 (1937). 
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Solar Absorption Lines between 2950 
and 2200 Angstroms 


E. Duranp, J. J. OBERLY, AND R. TouseEy 
Naval Research Laboratory, Washington, D. C. 
April 23, 1947 


OLAR spectrograms were taken by the Naval Research 
S Laboratory on the V-2 rocket flights of October 10, 
1946! and March 7, 1947. The spectra obtained at 55 km 
and 75 km on the respective dates have been studied to 
jdentify the Fraunhofer absorption lines in the region be- 
tween 3000 and 2200A. Some 300 observable absorption 
features were compared to a master finding list of 1100 of 
the principal classified lines of the arc and spark spectra of 
elements 1 to 30. The heavier elements will be added as 
time permits. The finding list was prepared in collaboration 
with Dr. Charlotte Moore-Sitterly, who generously made 
available the unpublished multiplet lists which are being 
compiled as an extension of the Revised Multiplet Table? 
to cover the rocket ultraviolet region. Entries included 
laboratory intensities, assigned multiplet numbers, and, in 
many cases, her solar intensity predictions. 

At the available resolution of about 1A, nearly every 
observed line is a blend; as many as 10 possible contributors 
to a single line have been found. Likely contributors have 
been assigned to nearly every observed feature, and work 
is under way to estimate their relative importance. 

As in the previously known region, Fe J and Fe JJ are 
dominant and clearly contribute to a majority of the lines. 
In Fig. 1, the matching of many of the solar lines (indicated 
by dots) to the Fe arc spectrum is apparent. In the regions 
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just below 2750, 2630, 2550, and 2490A there is a piling up 
of intense iron lines. This causes the whole level of the 
solar radiation to fall off sharply as shown in Fig. 1. In 
addition, many strong single Fe lines are found throughout 
the spectrum. 

The great Mg JJ lines at 2803 and 2796A are of particu- 
lar interest. They appear as two bright emission lines in the 
center of a great absorption band running from 2775 to 
2825A. 

Dr. Menzel of the Harvard Observatory offered the 
above explanation of the observed spectra and suggested 
the following explanation: A strong eruption of hydrogen 
occurred about 1 hour before the rocket flight, and the 
Mg emission may have originated in the prominence. 
However, similar emission lines, unresolved, may be in- 
ferred from the October 10 spectra, at which time no im- 
portant prominences are known to have occurred. 

Twelve lines of Si J of great intensity were found. 
The one at 2882 and the group between 2507 and 2529A 
showed strong wings. The existence of wings on the group 
between 2208 and 2219A could not be proved. A strong 
line of C J was found at 2478A. 

In addition to the above elements one or more lines have 
been assigned as follows: definite, V J, V IJ, Cr IJ, Mn JI, 
Co I, and Al IJ; probable, Na J, Ni J, Ni JJ, Cr J, Co JI, 
Be J, and Al J; possible, P J, and Cu J. More definite as- 
signments will be made taking into account multiplet in- 
tensity and the relative abundances of the elements. 

There appear to be regions of general absorption between 
2886 and 2893A and between 2442 and 2472A. The finding 
list contains few lines in these regions; the absorption may, 
therefore, be molecular. 

A complete report of the analysis will be published later. 

1W. A. Baum, F. S. Johnson, J. J. Oberly, C. C. Rockwood, C. V . 
Strain, and R. Tousey, Phys. Rev. 70, 781 (1946) 
? Charlotte 


E. Moore, “A revised multiplet table of astrophysical 
interest,"’ Princeton Observatory. 





Lower Bounds for Eigenvalues 


R. J. Durrin 
Depariment of Mathematics, Carnegie Institute of Technology, 
Pittsburgh, Pennsyloania 
April 28, 1947 


ARTA! has given a method for finding lower bounds 

for the “ground state” eigenvalue of a vibrating 

membrane. He does not claim that his method is as nice 

as the Rayleigh-Ritz method for finding upper bounds, 

but in a numerical example he shows that fair accuracy is 
obtainable. 

The purpose of this note is threefold: (1) We extend 
Barta’s method to include the Schrédinger equation, (2) 
the boundary conditions of Barta are relaxed, for instead 
of «=0 on the boundary we have u>0, and this may be 
easier to apply in practice, (3) we give a very simple, but 
mathematically rigorous, proof. Barta makes use of the 
fact that the ground state eigenfunction is positive; the 
proof of this fact is very long and sophisticated. 

Consider a continuous real function, ¥, which vanishes 
on the secondary of a region, R, and satisfies the -dimen- 
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sional Schrédinger equation, AY+(A— V)y=0. For sim- 
plicity we assume that V is a finite function. Let u be any 
function with bounded second derivatives satisfying u>e 
for a positive constant, «. The continuous function r=y/u 
vanishes on the boundary and may be assumed to have a 
positive maximum in R (otherwise take —y). At a maxi- 
mum point the first derivatives of r vanish, and the second 
derivatives are negative or zero. Adding the second deriva- 
tives gives Ar=(uAy—yAu)/u*<0. Thus Au/u>Ay/y 
= V—Aso that \> V—Au/u. A fortiori, \>min.(V—Au/u) 
for all points of the region. This gives the theorem: The 
eigenvalues of Schrédinger’s equation satisfy the inequality 
A>min.(V—Au/u) where u is a positive function in the 
region. 

The rest of this note is mainly of mathematical interest 
and proves that on the boundary, S, it is only necessary 
to assume u>0. In case the region is finite, we suppose 
that S is sufficiently regular so that there is a function, g, 
which satisfies Ag+(h— V)g<0 on Sand g>1 in R, where 
h=min.(V—Au/u). If w=u+eg, then A\>min.(V—Aw/w) 
by what has just been shown. Clearly, V--(Au+edg)/ 
(ut+eg)>V+(hu— Vu—edg)/(ut+eg)=h+e(Vg—hg—Ag)/ 
(u-+-eg). In a neighborhood of the boundary the last term 
is positive, and in the remainder of the region it must 
vanish with ¢«, so that A>h, Q. E. D. 

In an infinite region we take g=1 in the above proof, 
that is, w=u-+e. Assume, moreover, that outside a suffi- 
ciently large sphere h— V <0. Thus V—Aw/w>h+e(V—h)/ 
(u-+-eg). Outside the sphere the last term is positive, and 
inside the sphere it vanishes with ¢, so that A>h. 

The problem discussed here could be generalized for an 
equation in which is multiplied by a positive function. 
On the other hand, the writer has been unable to treat by 
similar methods the eigenvalues of the acoustic equation, 
the equation of the clamped plate, or Maxwell's equations. 


1 J. Barta, Comptes rendus, 204, 472 (1937). 





On a Connection between the Fountain Effect, 
Second Sound, and Thermal Conductivity 
in Liquid Helium II 
LoTHAR MEYER 
Institute for the Study of Metals, University of Chicago, 
Chicago, Illinois 
AND 
WILLIAM BAND 
Institute for Nuclear Studies, University of Chicago, 
Chicago, Illinois 
April 25, 1947 


E have investigated the possible phenomenological 

relations between the observed properties of 

liquid helium II. We assume that in the absence of a 

constriction the internal forces responsible for the fountain 

effect are still present but produce an internal momentum 
density M such that 


uVT=Vp=—dM/dt, (1) 


where JT is the temperature and p is the fountain effect 
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pressure difference. We then apply the first law of thermo. 
dynamics to a heat transported with the momentum 


V-h=V-(«M/m)=— pCdT/dt. (2) 
We solve (1) and (2) simultaneously to give a wave equa- 


tion for momentum and temperature similar to second 
sound. The phase velocity is 


V, = (eu/mpC,)}, (3) 


where y» is the fountain effect coefficient,'? pC, jg the 
specific heat* per cm*, and e is the heat associated with a 
mass m transported with the momentum. Empirical 
values of u, V. (second sound velocity),*® and pC, below 
1.8°K yield values of ¢/m which fit a monotonically jn. 
creasing curve from zero at 0°K to the heat content at 
the lambda point (7)). 
From the equation 


d(epz/m)/dT = pC, (4) 


we calculated the concentration p;/p of the liquid taking 
part in the momentum. The curve connecting these values 
is approximately equal to (7'/7)),* going therefore through 
zero at 0°K and unity at 7), which it should do from general 
considerations. 

The observed finite heat conductivity in liquid helium 
requires a wave equation having a relaxation (damping) 
term added to (1), namely, 


@T  1dT 
vern 1 etal /Y # 


Combined with the data on thermal conductivity® this 
gives a time of relaxation + of the order 10-* to 10~ sec., 
going to zero at 7). Equations (5) and (1) with the damp- 
ing term lead to an alternative form of the heat current as 
a convection process: 


h= pQuz, (6) 


where pQ is the heat content per cm? of liquid helium, and 
uz is the velocity associated with the steady momentum 
density as limited by the relaxation term. The data on 
thermal flow give u, as a function of T and d7’/dx, leading 
to values practically identical with those of the super 
critical velocities through slits.? ’ 

The relaxation causes an absorption and dispersion of 
second sound near the lambda point, formally analogous to 
the transmission of electromagnetic waves in a conducting 
medium. This may be responsible for the rapid approach to 
zero of the observed phase velocity as 7 approaches 7}, 
and leads to the expectation that under ideal conditions 
(small amplitude and high frequency) the velocity of 
second sound should rise up to the immediate neighborhood 
of 7) just as the maximum in the curves of the heat con- 
ductivity and the fountain effect shift towards the lambda 
point under appropriate conditions.* 

We are studying the physical meaning of these relations. 

1H. London, Proc. Roy. Soc. London A171, 484 (1939). 
_2L, Meyer and J. H. Mellink, Physica, 's-Gravenhage (to be pub- 
WH. Keesom, Helium (Amsterdam, 1942). 

4V. P. Peshkov, Vestnik Akademii Nauk 4, 117 (1945); Nature 157, 
200 Ferra k, W. Fairbank, and C. T. Lane, Phys. Rev. 71, 477(A) 
(1947). airDank, . Fairbank, an - i. ne, ys. Nev. , 
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Errata: The Collision of Neutrons with 
Deuterons and the Reality of 
Exchange Forces 
(Phys. Rev. 71, 558 (1947)] 

H. S. W. Massey AND R. A. BUCKINGHAM 
University College, London, England 
AND 
Conductivity of Sodium-Ammonia Solutions 


(Phys. Rev. 71, 563 (1947)] 


D. K. C. MacDoNALD AND K. MENDELSSOHN 
Clarendon Laboratory, Oxford, England 
AND 
A. J. Bircn 
Dyson Perrins Laboratory, Oxford, England 


HE Editor regrets that Fig. 2 of the first of the above 
papers has been interchanged with Fig. 1 of the sec- 
ond paper. The captions, however, are right, as they stand. 





_Erratum: The Quadrupole Moments and 
Spins of Br, Cl, and N Nuclei 


[Phys. Rev. 71, 644 (1947)] 


C. H. Townes, A. N. HoLDEN, J. BARDEEN, AND F. R, MERRITT 
Bell Telephone Laboratories, Murray Hill, New Jersey 


HE Editor regrets that Figs. 1 and 2 of the above 
paper have been interchanged. The captions for the 
figures, however, are correct as they stand. 





Isotopic Changes in Cadmium 
by Neutron Absorption 
A. J. DEMPSTER 


Argonne National Laboratory, Chicago, Illinois 
May 1, 1947 


N a recent letter B. J. Moyer, B. Peters, and F. H. 

Schmidt! report that the absorption of thermal neu- 
trons in cadmium that had been enriched in the isotope at 
mass 113 was six times normal, and conclude that this 
isotope is mainly responsible for the large absorption in 
this element. 

The mass spectra reproduced in Fig. 1 illustrate another 
phenomenon that also indicates the absorbing isotope, 
namely, the gradual alteration of the isotopic constitution 
of an element by long exposure to neutrons. The upper 
mass spectrum shows the normal isotopes at masses 110, 
111, 112, 113, and 114, with Cd"*=12.3, and Cd™*=28 
percent. The faint isotopes at 106 and 108 do not show on 
the prints. The lower spectrum shows the changed iso- 
topic abundance in a sample scraped from the surface of 
a piece of metal that had been used to absorb intense 
neutrons from a pile over a long period of time. These 
neutrons include a considerable number with energies 
above the thermal range. 

Photometric measurements, made by using Nier's 
measurements of the normal cadmium abundances as 
standards, showed that the isotope at mass 113 was reduced 
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Fic. 1. Normal cadmium (above) and isotopes altered 
by neutron absorption (below). 


to 1.6+0.2 percent, about one-eighth of its normal abun- 
dance, and that the isotope at mass 114 was increased to 
39.5+1.5 percent. Within the limits of the accuracy at- 
tained, the increase in 114 is equal to the decrease in 113. 
Any changes in the abundances of the isotopes at masses 
110, 111, 112, and 116 were less than five percent. It was 
noticed, however, that in both samples the isotopes at 110 
and 111 had a greater difference in abundance than is 
indicated by Nier’s measurements (12.8 and 13.0 percent). 
We may conclude that the other main isotopes have ab- 
sorbing cross sections less than one-fortieth of that of 
Cd"*, and that if any (m, 2m) or other reactions occur, 
they are comparatively rare. 


eae” Moyer, B. Peters, and F. H. Schmidt, Phys. Rev. 69, 666 





Radiation from Electrons in a Synchrotron 


F. R. Evper, A. M. Gurewitsca, R. V. LANGMutR, 
AND H. C. PoLiock 
Research Laboratory, General Electric Company, 
Schenectady, New York 
May 7, 1947 


IGH energy electrons which are subjected to large 

accelerations normal to their velocity should radiate 
electromagnetic energy.’~* The radiation from electrons 
in a betatron or synchrotron should be emitted in a narrow 
cone tangent to the electron orbit, and its spectrum should 
extend into the visible region. This radiation has now been 
observed visually in the General Electric 70-Mev synchro- 
tron. This machine has an electron orbit radius of 29.3 
cm and a peak magnetic field of 8100 gausses. The radiation 
is seen as a small spot of brilliant white light by an ob- 
server looking into the vacuum tube tangent to the orbit 
and toward the approaching electrons. The light is quite 
bright when the x-ray output of the machine at 70 Mev 
is 50 roentgens per minute at one meter from the target 
and can still be observed in daylight at outputs as low as 
0.1 roentgen. 

The synchrotron x-ray beam is obtained by turning off 
the r-f accelerating resonator and permitting subsequent 
changes in the field of the magnet to change the electron 
orbit radius so as to contract or expand the beam to suit- 
able targets. If the electrons are contracted to a target at 
successively higher energies, the intensity of the light radia- 
tion is observed to increase rapidly with electron energy. 
If, however, the electrons are kept in the beam past the 
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peak of the magnetic field and then expanded to a target, 
the intensity of the radiated light appears to be indepen- 
dent of the energy at which the electrons are removed from 
the beam. This is to be expected, for in a given machine the 
radiation is proportional to the fourth power of the elec- 
tron energy. The light radiation is not observed if the beam 
is contracted before its energy is about 30 Mev. When the 
electron beam has been accelerated to the peak of the 
magnetic field and then decelerated to low energy, a 
rough measurement of the phase angle over which the light 
was visible gave a value of 90-100 degrees. The light was 
viewed through a slotted disk rotating at synchronous 
speed. 

If the r-f resonator is turned off a short time before the 
peak of the magnetic field, the electron beam slowly con- 
tracts to a radius just larger than that of the interior 
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target and then expands as the magnetic field decreases 
In this case, the observer no longer sees a single point of 
light but a short line with extension in the plane of the 
orbit. 

The light emitted from the beam is polarized with the 
electric vector parallel to the plane of the electron orbit 
It disappears as the observer rotates a piece of Polaroid 
before the eye through ninety degrees. An investigation of 
the spectral distribution of the energy is in progress and 
will be reported. 

This work has been supported by the Office of Naval 
Research under contract N5ori-178. 

1D. Iwanenko and I. Pomeranchuk, Phys. Rev. 65, 343 (1944), 

+L. 1. Schifl Rev. Sci. Inst. 17, 6 (1946), 


4J. S. Schwinger, Phys. Rev. 70, 798 (1946). 
5H. C. Pollock et al., Phys. Rev. 70, 798 (1946). 
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Proceedings of the Southeastern Section of The American Physical Society 


MINUTES OF THE SALISBURY, NoRTH CAROLINA, MEETING ON APRIL 4-5, 1947 


HE thirteenth annual meeting of the South- 

eastern Section of the American Physical 
Society was held at Catawba College, Salisbury, 
North Carolina, on Friday and Saturday, April 
4-5, 1947. Approximately 180 members and 
guests attended. Local arrangements were made 
by a committee headed by M. L. Braun. 

The regular program consisted of thirty-nine 
papers, abstracts for thirty-one of which are 
appended hereto. Abstracts of the other papers 
will be found in the July-August issue of the 
American Journal of Physics. One abstract did 
not reach the Program Committee in time to be 
included in the printed program, but the paper 
was presented at the meeting. The title was, 
“On the Division of Nuclear Charge in Fission,”’ 
by R. D. Present, of the University of Tennessee. 

As additional features the following invited 
papers were presented: “Some Aspects of the 


Theory of Neutron Diffraction by Crystals,” by 
Fredrick Seitz; ‘New Developments in Instru- 
ments for Counting and Detecting Nuclear 
Particles,” by Walter Jordan; ‘‘Microwave Spec- 
tra,’ by Walter Gordy; “Report on the Oak 
Ridge Institute of Nuclear Studies,” by W. G. 
Pollard; “Educational Programs in Atomic 
Energy for the Southeast,’”’ by Karl Morgan; 
and “The Training School Program at Clinton 
Laboratories,”’ by Fredrick Seitz. 

At the business meeting the election of the 
following officers for the Section was announced: 
Chairman, M. L. Braun; Vice-Chairman, E. H. 
Dixon; Secretary, Eric Rodgers; Treasurer, H. 
F. Henry; and Member of the Executive Com- 
mittee, C. B. Crawley. 

The 1948 meeting will be held at the Uni- 
versity of Tennessee in Knoxville. 

Eric RopGers, Secretary 


ABSTRACTS 


1. An Automatically Recording Grating Spectrograph 
for the Infra-Red. ALvin H. NIELSEN, The University of 
Tennessee, Knoxville.—A recording prism-grating infra-red 
spectrograph with a spectral range of 1-25y has been built 
in the University of Tennessee shop. Spherical mirrors are 
employed with the exception of the principal collimator, 
which is an off-axis paraboloid, and the final mirror which 
focusses a reduced slit image on the thermocouple. An 
airtight cover permits the removal of carbon dioxide and 
water vapor from the spectrograph. All controls emerge 
from the base so the instrument need be opened only when 
changing absorption cells or gratings. The grating is 
mounted on a ball-bearing shaft and is turned slowly by a 
motor and worm gear arrangement. A counter attached to 
the drive trips a special pen on the recorder and thus puts 
a wave-length scale on the chart. A vacuum thermocouple 
on which the final slit image is focussed causes a high- 
sensitivity Leeds & Northrup galvanometer to deflect a 
beam of light across two photo-cells. The differential output 
of these cells is amplified by a d.c. amplifier described by 
Bell, Noble, and Nielsen.! The final amplified signal is 
reproduced on charts by either a Bristol recording milliam- 
meter, or Leeds & Northrup Speedomax Type G. Slides 
illustrating the details of the instrument and its per- 
formance will be shown. 


1E. E. Bell, R. H. Noble, H. H. Nielsen, Rev. Sci. Inst. 18, 48 (1947). 
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2. Mechanism of Luminescence of Infra~Red-Sensitive 
Phosphors. Ferp E. WiLiiams, University of North 
Carolina.—The model of potential energy versus configura- 
tion coordinate proposed by Williams and Eyring' to 
correlate diverse properties of luminescent solids has been 
found applicable to infrared-sensitive phosphors. An 
explanation of the complex afterglow of the rare earth- 
activated strontium sulfoselenides is that electron traps of 
various depths are active. Calculations from the model 
indicate that at high temperatures only the deepest traps 
are active, therefore, assuming a predominantly mono- 
molecular rate-determining ‘process, the afterglow should 
be exponential. This is verified by an analysis of the data 
of Shrader* for SrS:Sm, Eu, and the thermal trap depths 
are found to be related to the optical depths as expected. 
As also predicted from the model, pronounced temperature- 
dependence of luminescent efficiency is encountered and 
can be precisely reproduced assuming an activated step 
in the luminescent process competing with the radiationless 
recombination of the excited electron and center. The fre- 
quency factors of the specific rate constants are small, 
indicating changes in multiplicity. The contributions of 
first- and second-order rate processes and the probable trap 
distribution are computed from buildup, afterglow,? and 
stimulation’ data. 


iF. E. Williams and H. Eyring, Phys. Rev. 69, 257 (1946). 
2 R. E. Shrader, Paper 15 at Am. Opt. Soc. Mtg. in Cleveland (March 


1946). 
+R. T. Ellickson and W. L. Parker, Phys. Rev. 70, 290 (1946). 








3. The Evaporation of Gallium. S. K. Haynes, Clinton 
Laboratories and Vanderbilt University—A thin B-ray 
source was made by evaporating gallium onto a thin sheet 
of polystyrene. Since gallium was not among the metals for 
which evaporation procedures were given by Olsen, Smith, 
and Crittenden!’ an investigation was made to find the best 
procedure. A modification of the procedure for aluminum? 
was finally used. The gallium was warmed just above its 
melting point and “buttered” onto a cooler clean open 
helix of 30 mil tungsten wire, drop by drop, until the 
tungsten was almost covered. When the tungsten was 
flashed to a dull red in vacuum, a fraction of the gallium 
reacted with it, while the remainder dropped off. The 
tungsten retained all the gallium ready for the final 
evaporation after a second “buttering” and flashing. 
Several other methods were tried with less success, includ- 
ing a close wound conical helix of tungsten and a tantalum 
crucible heated by a filament. 


1 j: App. Phys. 16, 425 (1945). 
2 J. Strong, Proc. in Exp. Phys. (Prentice Hall), p. 168. 


4. A Mechanism for Ductile Failure in a Notched Steel 
Bar.* Paut E. SHEARIN, University of North Carolina.— 
A mild-steel bend bar, 0.75’ 0.75" X4.6”", having four 
15 percent Izod notches undergoes a mode of failure in the 
slow-bend test in which the head of the advancing crack 
remains sensibly parallel to the rootline of the initiating 
notch, when failing in a ductile manner.' In such a test 
the work associated with the fracturing of a unit-cross- 
sectional area (unit work) is to a first-order approximation 
a linear function of the distance to the neutral axis, 
(a—x)/2, where x represents the distance the crack has 
traveled and a is the original retained depth of the specimen 
at the roots of the notches. This is so since, as the crack 
opens, a wedge of material between the crack front and the 
neutral axis is undergoing plastic deformation. As the 
crack progresses the height of this wedge decreases. In- 
tegration of f[(¢—x)/2] yields the result that unit work 
is proportional to the mean distance of the area from the 
final crack terminus. Experimental results indicate the 
essential correctness of this result. 


* Work done under contract ONR, NRL, N6ori-227. 
1C. Zener and J. H. Holloman, J. App. Phys. 17, 69 (1946). 


5. Voltaic Effects in Waterglass. E. Scott Barr, 
Tulane University.—Sodium silicate in the waterglass form 
shows two voltaic effects which make it a complicating 
factor when used as a binder in cements which come into 
contact with electrical measuring leads. It has been found 
that during the drying period of the material there are 
voltage differences between different portions of the 
material sufficient to produce appreciable indications on a 
sensitive galvanometer. Furthermore, during the drying 
period, and apparently to some extent after drying, water- 
glass polarizes very noticeably. This latter effect precludes 
accurate potentiometer measurements of developed e.m.f., 
and measurements of resistance with a Wheatstone bridge 
or ohmmeter. 

In the work now in progress, two identical wire windings, 
of copper or nickel, are put at the ends of a glass plate 
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which is then coated with a thin layer of waterglase, | 
has been found that the winding in the more moist sand 
of the film is positive. Voltages up to 50 millivolts and 
more have been observed between the windings, 

effects and the decay of voltage after the film has been 
polarized by an external voltage source are to be discussed. 


6. Conductivity of Sodium-Ammonia Solutions. MELyyy 
EISNER AND JOSEPH W. STRALEY, The University of Nor 
Carolina.—Ogg' and Hodgins* have reported the existence 
of persistent currents in dilute solutions of Na in ammonia 
at liquid air temperatures. Boorse* and others have been 
unable to confirm these observations. An attempt is being 
made to reproduce these experiments. The apparatus con. 
sists of a search coil of 1010 turns wound on a Bakelite 
spool whose inside diameter is 1 cm. This coil is capable 
of detecting the magnetic field due to a current change in 
a neighboring circuit of 0.4 ampere-turns. Two types of 
cells have been used, the first consisting of slide cover 
separated by a mica spacer as described by Ogg, the second 
being a glass toroid of 2 cm diameter. Of seventy-five tests 
forty have shown negative results. Although slight de. 
flections were thought to be observed in the remaining 
tests, they were not large enough to be distinguished with 
certainty from the random motions of the galvanometer, 
Work is being continued. 

1 Phys. Rev. 70, 93 (1946). 


? Phys. Rev. 70, 568 (1946). 
3 Phys. Rev. 70, 92 (1946). 


7. Diffuse Reflection from Dielectrics. Wa. A. Rensz, 
Louisiana State University—When light strikes an un- 
polished surface of a dielectric, some of it will be absorbed 
or transmitted and the rest will be sent back into the 
original medium in many or all directions. Diffuse reflection 
is said to occur. Previous theories of diffuse reflection either 
have been of a semi-empirical nature, or have limited the 
reflecting mechanism to that of elementary regular mirrors 
oriented in various ways on the illuminated surface. 

A new theory is here developed which treats the diffusely 
reflected light as being made up of three parts: (1) light 
regularly reflected from elementary mirrors, (2) light 
scattered from edges and points on the surface, (3) light 
which, after entering the surface and suffering absorption, 
has been deflected back into the incident medium in some 
manner. Fresnel’s equations from Maxwell's electromag- 
netic theory and classical equations for the scattering of 
light by an anisotropic molecule permit explicit expression 
of (1) and (2). Simplifying assumptions allow (3) to be 
estimated. The final expressions contain formulae for 
finding the total intensity and the intensities of the 
polarized components in terms of six parameters related 
directly to the roughness of the surface and to the re 
fractive index and absorption properties of the dielectric. 


8. Positive Grid Characteristics of a Triode. Gzorce W. 
Woop, Louisiana State University—Theory concerning 
the triode is meager. However, a recent treatment of the 
diode by Jaffé, wherein the electrons are assumed to have 
uniform initial velocities, provides a means of extending 
our knowledge of the triode. By operating the grid of a 
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triode at a positive potential, it becomes a virtual cathode 
discharging electrons with uniform initial velocity and thus 
the theory becomes applicable. The positive grid charac- 
teristics of a number of radio receiving triodes are deter- 
mined and matched with those determined from the theory. 
An agreement, satisfactory under given conditions, is 
found in the region of low plate voltages. 

of the empirical curves from the theoretical 
ones is appreciable because the practical tubes tested do 
not, by any means, correspond to the ideal conditions of 
the theory: the electrodes are not plane; the electrostatic 
field is not homogeneous; secondary emission from the grid 
is present; contact potentials exist; and there is always a 
velocity distribution among the electrons. While the inves- 
tigations described cannot be taken as conclusive, they 
point the way towards a better theory of the triode. 


9. Effect of Scattering on Experimental Determinations 
of XR for High Energy Particles. F. T. Rocers, Jr., The 
University of North Carolina.—In certain apparatuses for 
determining XR for a high energy particle (R is the radius 
of curvature which the path would have if the particle were 
to traverse a suitable electrostatic field X perpendicularly 
to the X-vector), X may be available in a convenient non- 
ideal configuration, R may be calculated from the par- 
ticle’s deflection (s) in this X, and uncertainties' may be 
introduced into the measured deflection by a thin scat- 
tering medium in the particle’s path before entering X. 
Thus if the proper X-value can be extracted from the 
available configuration, then the uncertainty in XR is 


ew =OavL/s+As/s, (1) 


where @a4y is a suitable measure of scattering,? L is a 
dimension, and As is a measure of errors of “triangulation,” 
both terms being essentially Gaussian. A discussion will be 
presented of the estimated lower limit to which ey might 
be reducible. 


1 The author is indebted to Dr. J. C. Street for having pointed out 
the im; ce of the @4y-contribution to ey. 
?E. J. Williams, Phys. Rev. 58, 292 (1940). 


10. Two Non-Linear Springs. G. W. Crawrorp, North 
Carolina State College of Agriculture and Engineering, AND 
F. T. RoGers, Jr., The University of North Carolina.—The 
equation 

T=t+46(L—1)+A(L—-l¥, (1) 


where ¢, 5, and A are constants and where (L—/) is the 
extension beyond a minimum length (/), has been proposed! 
for the tension of strings (about 0.0004-in. diameter) in 
certain commercial string-galvanometers. To exhibit this 
concept on a large scale, two plane springs of 0.13 cm 
diameter phosphor-bronze wire were made in grid form; 
nine grid elements were used, each 0.75’ X 1.5”; they were 
rectangular in one, not in the other. The constants ¢, 5, and 
4 were found to be, respectively, : 


0, 1.82105 dynes/em, —0.0586X105 dynes/cm*, (2a) 
0, 6.1105 dynes/cm, 1.68X105 dynes/cm*, (2b) 


and represented the non-linearities to within A.D.’s of 
2X10‘ and 6X 10* dynes, respectively; these representa- 
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tions are valid for (L—1) up to 10 cm and 4 cm, respéc- 
tively, which extensions approached the elastic limits of the 
two springs. 

iF, T. Rogers, Jr., Rev. Sci. Inst. 12, 351 (1941). 


11. A Redetermination of the Nuclear Magnetic 
Moment of Sodium by the Polarization of Resonance 
Radiation. MarGcuerireE M. Rocers, The University of 
North Carolina.—A redetermination of the nuclear mag- 
netic moment of sodium was made by measuring the degree 
of polarization in the resonance radiation from sodium 
vapor subjected to various magnetic fields. This effect 
is well known and has been previously measured by 
Larrick' and by Ellett and Heydenburg.? In the present 
experiment, great care was exercised to reduce the self- 
reversal in the sodium source by means of a specially 
designed discharge tube. 

The variation of the degree of polarization with field 
strength was compared with theoretical values from com- 
putation of the Paschen-Back intensities of the hyperfine 
structure based on a spin of §. The values of the Landé 
g-factor and nuclear magnetic moment resulting from 
fitting the observed curve to the theoretical one gave good 
agreement, within the accuracy of the experiment, to the 
accepted, precise values of Millman and Kusch of 1.4765 
and 2.215 


1 Larrick, Phys. Rev. 46, 581 (1934). 
? Ellett and Heydenburg, Phys. Rev. 46, 583 (1934). 


12. Quantum Mechanics of Vibrating String. NaTHAN 
RosEN, The University of North Carolina.—As a simple 
example of the quantization of a continuous medium, the 
transverse oscillations of a perfectly flexible string stretched 
under tension between two fixed supports are treated by 
means of wave mechanics. The energy levels and wave 
functions of the stationary states are determined. One 
obtains the well-known infinite zero-point energy, which 
would have to be discarded to give finite values for the 
energy of the string. The degree of degeneracy of the 
system is considered. 


13. Problem of Homogeneous Incompressible Fluid 
Sphere in General Relativity. BRowntze NEWMAN AND 
NATHAN RosEN, University of North Carolina.—The 
Schwarzschild and related solutions of the gravitational 
field equations for the problem of the homogeneous sphere 
impose the conditions 7'=77=7=—p and T=p 
=constant. Eddington! pointed out that for an incom- 
pressible homogeneous fluid the last condition should be 
replaced by 7.*=constant. An attempt has been made 
to solve the problem for this case. Starting with a line 
element of the Schwarzschild form, one obtains a set of 
field equations which can be reduced to a single equation 
with one unknown. This equation is, however, of a rather 
complicated nature, and thus far it has not been possible 
to find its general solution in closed form. Two particular 
solutions for the interior of the sphere can be obtained, but 
they are not satisfactory from the physical standpoint. It 
may be that the problem will have to be solved by nu- 
merical integration of the differential equation. 


1A. S. Eddington, Mathematical Theory of Relativity (Cambridge, 
1923), p. 170. 











14. The Limits of Visibility and the Maximum Bright- 
ness of Continuous Spectra at Various Color Temperatures. 
Otto STUHLMAN, JR., University of North Carolina.— 
Preliminary data on the visibility of radiation at low inten- 
sities and low color temperatures indicate that the sensi- 
bility to threshold brightness conforms to the results of 
Weaver and Hussing (J. O. S. A. 16, 29, 1939). At higher 
color temperatures the visibility departs from the above 
results. The maximum brightness wave-length parallels the 
blue end of the spectrum. 


15. Preliminary Experiment on Convection Currents in 
a Porous Medium. H. L. Morrison, North Carolina State 
College of Agriculture and Engineering, AND F. T. ROGERs, 
Jr., The University of North Carolina.—To assess a theory 
recently advanced! for the minimum temperature gradient 
necessary to cause convection of a liquid in a porous 
medium, an experiment has been initiated using glycerol 
in a sand having a permeability of about 95 darcys. The 
mixture was placed in a beaker to a depth of 10 cm; ther- 
mocouples were used to observe temperatures; heat was 
applied from below; and microscopic observations were 
made of suspended impurities. When the particles of 
impurity moved, the temperature gradient was about 80°C 
locally, but varied widely throughout the depth of the 
mixture. The theory clearly calls for a gradient of excess 
of 200°C/cm if uniform and if boundary conditions such 
as those in the experiment are negligible (there is ap- 
parently no information on the viscosity of glycerol above 
30°C, so that logarithmic extrapolation is being used to 
higher temperatures until suitable measurements can be 
completed). Weighing all these uncertain factors together, 
it can only be concluded at present that the theory predicts 
excessive minimum gradients, a conclusion which was 
anticipated in the presentation of the theory. 


1C. W. Horton and F. T. Rogers, Jr., J. App. Phys. 16, 367 
(1945). 


16. Alternative Basic Theories for the Carbon-Pile Type 
of Voltage-Regulator. CHARLES W. Duss anp F. T. 
Rocers, Jr., The University of North Carolina.—A second 
theory for the carbon-pile type of voltage-regulator,' 
developed on the assumption of coulombic friction, predicts 
one or more half-cycles of cosine waves in the output- 
voltage for a step-function change of input-voltage; this 
is in disagreement with the apparently exponential wave- 
form which is observed? and which the first theory 
predicts. A third theory, based on coulombic friction and 
neglecting armature inertia as compared with inductive 
reactance, leads to an exponential character in the output- 
voltage in which the time-constant suddenly decreases 
(recovers, sharply) when the armature begins to move. 
All three theories have been compared to a minimum of 
experimental data: for a 16 percent step applied to a 
certain model of regulator, they predict 1.9 percent, zero 
to 16 percent, and 1.9 percent to 16 percent changes 
(after transients) in output-voltage, while the measured 
change is 3.6 percent. All three theories appear in their 
present forms to predict response-times at least an order 
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of magnitude smaller than the 30 milliseconds observed for 
this model. 


1W. G. Neild, Trans. AIEE, 63, 839 (1944), 
2F. T. Rogers, Jr., Instruments 19, 256 and 381 (1946), 


17. Measurement of the Ionization Potential of Gases 
in Photoelectric Cells. ARTHUR WALTNER, The University 
of North Carolina.—When an alternating e.m.f. is applied 
across a gas-filled photoelectric cell the current wave form 
passes through three distinct phases as the Voltage js 
increased: (1) For very low voltages (2 or 3 volts) a half. 
wave rectified pattern is observed; (2) as saturation sets 
in a rather close approximation of a square wave appears; 
(3) when the voltage is sufficiently great to produce ions 
by collision an “ionization hump” appears at the top of the 
square wave. By observing the voltage at which this hump 
first appears, it is possible to approximate the ionization 
potential of the gas. A more accurate determination cap 
be made by photographing the wave form and measuring 
the phase angle at which the hump begins and ends, Cor. 
rections must be made for the finite velocity of the elec. 
trons as they leave the cathode. This correction is made by 
potentiometric measurement. Values which are correct to 
within 5 percent of the accepted values have been obtained, 


18. Magnetic Facilities Provided by the Kensington 
Magnetic Laboratory. B. P. Ramsay AnD F. L. Yost, 
Naval Ordnance Laboratory.—During the war, the Carnegie 
Institution of Washington operated for the Naval Ord- 
nance Laboratory special facilities which were known as 
the Kensington Magnetic Laboratory. The facilities 
included coil systems of unusual size, design, and con- 
struction. These were employed to generate in the labora- 
tory magnetic conditions and signals which are useful in 
the development of Naval Ordnance. Generically, these 
coil systems have been called magnetic generators. Sole- 
noids, Helmholtz coils, and gradhelms* were specific types 
of generators in use. The largest solenoid was 3 ft.X4 ft. 
X40 ft. Several helmholtzes were framed in 20-ft. cubes; 
and gradhelms were comparable in size. In various appli- 
cations, signals from several types of magnetic generators 
were superimposed. The signals were controlled in contour, 
amplitude, period, and phase by simulators of various 
types. The periods of the signals which were generally 
employed fell in the range from 1 sec. to 1000 sec. Matters 
of general technical interest regarding construction, 
accuracy, and interaction will be discussed. 


aod H. Shortley and A. May, J. May, J. App. Phys. 16, 841-843 


19. A Photoelectric Nephelometer. N. UNDERWOOD 
AnD A. H. DoeRMANN, Vanderbilt University.—A nephe 
lometer is described that does not require a precision 
optical system or an electronic amplifier. It is sensitive 
enough to measure the growth of B. coli from a starting 
concentration of approximately 2X 10* bacteria per ml. 
The readings are reliable from 107 to 10° bacteria per ml. 
A graph of the photo-current due to the scattered light 
versus the number of B. coli per ml is given. The B. coli 
culture is a saturated 24 hr. nutrient broth culture in the 
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lag phase. In this range of bacterial concentrations the 
scattered light is proportional to the number of scattering 
centers. The instrument is well adapted to the study of 
growth curves of bacterial cultures and the lysis phenomena 
that occur when the organisms are attacked by bac- 


a ae beam of light passes through the vertical 
test tube containing the specimen. The light scattered 
downward is measured with a barrier layer photo-cell. The 
horizontal incident light produces a photo-current of 
2.4X 10-8 amperes. The light scattered by a specimen whose 
titer is 1.6 10° B. coli per ml produces a photo-current of 
1.2X10-* amperes. In this case the scattered light is 1/200 


of the incident light. 


20. Mass Measurement with a Single Field Mass Spec- 
trometer. E. P. Ney AnD A. K. MANN, University of Vir- 
ginia.—An improved Nier type 60 degree mass spectrom- 
eter has been used in the measurement of the mass ratio 
H?*/H.**. Measurements of this type have previously 
been made with “‘double focussing” or “velocity selecting”’ 
mass spectrographs.' The value of the H,‘/H;' ratio was 
computed from the measured H;'*+*+/H,‘** ratio and com- 
pared with Aston’s value.? The agreement was well within 
the accuracy of the measurement. This is significant 
because it indicates that precision mass measurements as 
well as relative abundance measurements can be made with 
a single field mass spectrometer of sufficiently high resolv- 
ing power. The primary advantage of this method lies in 
the comparatively short time required for the experiment. 
The fifty-six observations of the double separation of 
H?*+—H** were made in approximately 15 hours. Data 
and a mass spectrum will be shown. 


1K. T. Bainbridge and E. B. Jordan, Phys. Rev. 50, 282 (1936). 
2F. W. Aston, Nature 137, 357 (1936). 


21. The Production of Constant High Rotationa! Speed. 
J. W. Beams AnD J. L. Youna, III, University of Virginia. 
—A steel rotor is supported in a vacuum (less than 2 x 10-* 
mm Hg) by an electro-magnetic suspension and is spun by 
a magnetic field rotating at a fixed constant speed. The 
rotating magnetic field is produced by two pairs of coils in 
which the phase of the current differs by 90 degrees. The 
two pairs of coils are placed symmetrically around the rotor 
and at right angles to each other. The frequency of the 
current in the coils is controlled by a piezoelectric crystal 
circuit. The induced currents in the rotor start it spinning. 
Because of the very small rotor drag, the rotor continues 
to accelerate up to the frequency of the rotating field where 
it “locks in” and spins at the same frequency as the field, 
ie., it operates as an armature of an induction motor until 
it reaches full speed and then operates as a synchronous 
motor. Rotor hunting is found to damp out quickly so that 
the rotor speed is as constant as the frequency of the 
piezoelectric crystal circuit. With a 1.59 mm. rotor spinning 
at 10° r.p.s. the temperature rise was less than 6°C. The 
arrangement makes possible a very constant high rota- 
tional speed mirror. 


22. Neutron-Proton Interaction. H. M. MosELEY AND 
NaTHAN RosEN, University of North Carolina.—The poten- 
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tial energy of interaction between a neutron and proton in 
the triplet state is taken to be of the form —B/(e"/*—1). 
The Schrédinger equation with this potential energy is 
solved exactly for the ground state of the deuteron, giving 
a relation between B and a to account for the observed 
binding energy. The problem of neutron-proton scattering 
is then treated on the basis of a combination of *S and 'S 
scattering. For the singlet state the interaction potential 
is taken to have the same form as for the triplet but with a 
different coefficient B’. The latter is fitted to the scattering 
data at low energies. It is found that taking a=1X10-" 
cm, B=60 Mev, and B’=38 Mev (corresponding to an 
energy of about 60,000 ev for the virtual singlet state) gives 
satisfactory agreement with experiment. Further work 
using spin-dependent forces is planned. 


23. Neutron-Proton Scattering. MELVIN EISNER AND 
NATHAN ROsEN, The University of North Carolina.—Since 
the exact law of interaction between a neutron and a 
proton is not yet known, it is of interest to consider various 
possible interactions and compare predictions based on 
them with experiment. In the present case it is assumed 
that the interaction potential energy is of the form 
V=— Vo sech*(r/a), where Vo has different values for the 
triplet and singlet states. The constants were chosen to 
give the correct binding energy for the ground state of the 
deuteron and to give a total neutron-proton scattering 
cross section of 19.7X10- cm* at E=0.01 Mev. It is 
found that for a=1X10-" cm agreement with scattering 
data is obtained. For the triplet state Vo is about 112 Mev, 
for the singlet, about 77 Mev, corresponding to a virtual 
singlet level of about 70,000 ev. 


24. Tolerance Concentrations of Radioactive Sub- 
stances. KARL Z. MorGAN, Clinton Laboratories.—A series 
of equations is developed in a form convenient for use in 
determining the tolerance amount of radioactive material 
that can be permitted safely in the human body. The body 
organ that sustains the maximum damage from a given 
radioactive material determines the body tolerance. It is 
assumed that this body organ rids itself of the radioactive 
isotope according to a simple exponential function such 
that the effective half-life, 7, is given by the equation: 
T=TeTr/Te+Tr, in which Te is the body elimination 
half-life and Tr is the radioactive decay half-life. When a 
parent substance is taken into the body and later disin- 
tegrates into a daughter radioactive substance, the disin- 
regration rate of the parent is equal to the disintegration 
rate of the daughter when the disintegration rate of the 
daughter is a maximum. Graphs are given comparing these 
particle disintegration rate curves with the energy dis- 
sipation rate curves. These equations are used in deter- 
mining the tolerance concentration in air and in water, for 
external and internal exposure to a number of radio- 
isotopes. The calculated concentrations check rather 
closely with the experimental values that have been 
established. The maximum permissible radiation dose is 
assumed to be 0.1 roentgen per day for gamma radiation 
and 0.01 equivalent physical roentgens per day for internal 
exposure to alpha rays. The general safe tolerance concen- 
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trations determined for beta or gamma radioisotopes is 
10-74 curies per cc in air or 5X 10~ yw curies per cc in water 
and for alpha emitters it is 10-'y curies per cc in air or 
10-4 curies per cc in water (1p curie taken as 3.7 X< 10* 
dis/sec.). 


25. A Portable Instrument for the Detection of Fast 
Neutrons in the Presence of Gamma Radiation. Pau W. 
REINHARDT, Clinton Laboratories.—The instrument con- 
sists of two ionization chambers mounted on a yoke as- 
sembly and connected so that the unbalance ionization 
current can be read on a Lindemann electrometer. One of 
the chambers is paraffin-coated inside and filled with 
butane at 30 pounds pressure. The other is filled with 
argon at about 35 pounds pressure. The pressures in the 
two chambers are so adjusted that their gamma ray 
ionization currents balance. Recoil protons ejected from 
the paraffin and butane by fast neutrons produce a dif- 
ferential ionization current which is read on the electrom- 
eter. A fast neutron flux of about 5 neutrons per cm? per 
second can easily be detected. 

The electrometer circuit provides for 5 fixed electrometer 
sensitivities, and an adjustable sensitivity control for 
maintaining constant sensitivity over a long period of time. 
The collecting voltage is supplied from a 450-volt battery 
pack. Except for the ionization chambers, the entire 
assembly is housed in an aluminum box 3} X4}" X8}”. 


26. Preliminary Evaluation of “Special Fine Grain Par- 
ticle Film” for Use in Neutron Monitoring. J. S. CHEKA, 
Clinton Laboratories.—Eastman’s Special Fine Grain Par- 
ticle Safety Film is used at Clinton Laboratories to monitor 
personnel for neutron exposure. A portion of the film is 
shielded by cadmium, which stops thermal neutrons, and a 
part of it is exposed without a shield, intercepting both 
thermal and fast neutrons. 

Calibrations were made to determine the sensitivity of 
the emulsion for both thermal and fast neutrons in terms 
of proton tracks observable in a dark field under a micro- 
scope of 970 magnification. The ratio 3.01 10* n,/ob- 
served track was noted for the N"(n,p)C™ reaction, and 
the ratio, 2.05 X 10‘ n;/observed track was noted for proton 
recoils. Evaluating expected track densities, based on 
N-atom density for m» and H-atom density for my, the 
first process was found to be about 53 percent effective and 
the latter 46 percent effective in producing recognizable 
tracks of three grains or more in this emulsion. 

Based on tolerance values of 4700 m,/(cm? sec.) and 
266 my/(cm? sec.), one week's tolerance doses are indicted 
by 20 tracks/50 fields for mi», and 17 tracks/50 fields for my. 
However, since different batches of this film were found to 
differ in sensitivity, a calibration must be made on each 
batch before use. 


To be Read by Title 


27. Integrating sin*x/x. S. M. CHRISTIAN, Agnes Scott 
College.— 


z sin*t 


: — a = $[7+log-(2x) —Ci(2x)]. 


y =0.577,215,7 and Ci(z) is the cosine integral. 
Brown, M. S. Corrington, Tolminson Fort, J. K, . P, 
and others suggest that this formula gives quicker and more 
reliable results than the series previously used. 











x Integral | x Integral | x Integral 
0.1  0.004,992| 1  0.423,691 | 11 1.833.309 
a .019,867 2 1.052,246 | 12 1.896, 801 
aa .044,330 3 1.218,516 | 13 1,903,509 
4 .077,897 4 1.267,117 | 14 1,949,275 
Be .119,906 5 1.462,629 | 15 2,005,723 
6 .169,539 6 1.555,951 | 16 2.013,281 
7 .225,841 7 1.573,438 | 17 2.043,656 
8 .287,743 8 1.682,002 | 18 2.094,072 
9 .354,096 9 1.755,531 | 19 2.103,836 
10 —1.764,264 | 20 2,123,538 


1S. M. Christian, Phys. Rev. 69, 546 (1946). 


28. Differential Equation for the Over-All Behavior of g 
Servo-Operated Recording Potentiometer. F. T. Roczrs, 
JR., AND MarGueritE M. Rocers, The University of 
North Carolina.—A type of recording potentiometer which 
is widely used industrially, contains a continually running 
motor which slowly actuates the recording pen at regular 
intervals (~2 sec) through a galvanometer-controlled 
clutch. At the instant of clutch-engagement, the gal- 
vanometer indicates the difference between the input- 
terminal emf (V) and the recorder emf (v), and so assures 
that a predetermined fraction (k) of V—v is transmitted 
by the clutch in the time interval (At) as a Av. That is, the 
over-all ‘“‘smoothed” behavior of the instrument follows the 
differential equations 

dv/dt=k(V—v), -M£ V—v< M,)\ 
dv/dt=kM sgn(V—v), | V—v| > M. | (1) 


The second of these equations indicates the predetermined 
limit of response of the galvanometer. For a typical instru- 
ment of this type, with 0< v¢ 10 mv, the quantities M and 
k have been measured so as to allow the use of Eq. (1) in 
correcting v-records for instrumental lags: M=3 mv, and 
k=0.0005 sec. for low circuit-resistances. 


29. The Rational Design of the Computor for a Certain 
Class of Mechanical System Containing Feedback- 
Control.—F. T. RoGers, Jr., The University of North 
Carolina.—Consider a mechanical system in which a motor 
(constant speed wy) operates a load at speed w, through 
first a speed-controller (delivering speed w= F(@)-wy, 
where F(@) is specified or known) and then a “take-off” 
unit; here @, an input to the controller, essentially governs 
the ratio w/wy. If some fraction (f) of w be withdrawn 
through the take-off unit from ¢=0, and ¢= /fudt be 
converted to @ through a computor, @=G(¢@), the problem 
is so to design this computor (or function G) that w behaves 
in the predetermined way, w = (¢t). Clearly, the computor- 
design equation results from the elimination of ¢ between 


o=f f° 2(x)dx (1) 


and 


F(8)=2(¢)/om. (2) 
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An interesting simple case is that in which Q(t) is to be 
ential in time; here Eqs. (1)-(2) yield G linear in ¢ 
for F linear in @. 


30. Statistical Significance of Some Measurements on 
the Permeabilities of Unconsolidated Sands. H. L. 
Morrison, North Carolina State College of Agriculture and 
Engineering, AND F. T. Rocers, Jr., The University of 
North Carolina.—The permeabilities (k) of some sands 
have been repeatedly measured by the method of Terzaghi,} 
using the vertical percolation of water through 90-cm 
columns of the sands; the pressure head was carefully 
measured as a function of time in each case, but no other 
particular precautions were taken. The k-values for a 
local red sand of 40-45 mesh had a mean value of 63 
darcys and a standard deviation of 7.5. Mean k-values for 
the clean natural sand from White Lake (Bladen Co.), 
N. C., were 90.3 darcys for the 40-50 mesh component, 
105.5 for the unscreened sand, and 125.1 for the (remaining) 
30-40 mesh component; standard deviations in each case 
were between 1 percent and 2 percent. All these measure- 
ments were made with columns of 1.34 cm diameter; 
doubling the diameter in the case of the 30-40 mesh white 
sand changed the mean k-value by about 2 percent. The 
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above information clearly allows estimates to be made of 
the significance of a mean k-value for such sands as found 
by a severely restricted series of measurements of low 
precision. 

1M. Muskat, Flow of Homogeneous Fluids T: Porous Media 
(McGraw-Hill, 1937), p. 84. — 


31. On the Reduction of Data Which Are Values of a 
Definite Integral. F. T. Rocers, Jr., The University of 
North Carolina.—Data which represent physical observa- 
tions may, for example in measuring the work done by a 
force of varying magnitude, of necessity be in the form! 


kW= J * e(x)b(x) dx. (1) 


Here e(x) would be a sought-for “‘unit-work function” ; the 
constant & and the function b(x) would be known; and y 
and z would be available as variables. Clearly the function 
e(x) can be obtained from the data by differentiation with 
respect to the limit x or y, 


e(y)= —k(OW/dy)/b(y), e(s)=k(@W/dz)/b(z), (2) 
if the experiment is properly designed. 
1 See, for example, Phys. Rev. 61, 728 (1942). 


